The  Pennsylvania  State  University 

APPLIED  RESEARCH  LABORATORY 

P.O.  Box  30 

State  College,  PA  16804 


ACOUSTIC  EMISSIONS  FROM  UNSTEADY 
TRANSITIONAL  BOUNDARY  LAYER 
FLOW  STRUCTURES 

By 


R.  C.  Marboe 
G.  C.  Lauchle 


Technical  Report  No.  TR  00-002 
October  2000 


Supported  by:  L.  R.  Hettche,  Director 

Applied  Research  Laboratory  Applied  Research  Laboratory 


Approved  for  public  release;  distribution  unlimited 

Ultw  QTiiiiiTi'i  IS> ZQiZQ'jSgg)  4 

20001116  001 


•  '3->' 

-  :X 

"J 

-  V3& 


'  %$/: 


REPORT  DOCUMENTATION  PAOi 


form  Ags&F&tfi&d 

0*m  MOL  Q7Q4-0m 


PvSstte  rqv&Tlszq  to  cn*  ccMgeaow  orf  «  ea&notsd  tc  aw@ra§9  t  397  raqsstfwa.  m 

artfeavra  trt  fisaemsaw^  to  data  mrarfssL  an*  coottotift?  c*ti  to  cassettes*  <rf  «KtaraaaSo». 

ceiteaw  o<  *fer"w«a*.  «**wa«g  wmbc^^^usb^  to  torsten.  to  Wtetompm  *se4&ert an  * 
Qgra  Kffiogy,  S«l»  Aftte^ewTvA^ 31383-4383.  ^  »  to  Oflto  of  swg  ayggK.  Pag 


1.  Answer  USE  OS&Y  rt^s*®  1 2.  REPORT  OATS 

I  October  2000 


V*Q  tzas®  fey  rsssweg™  wsracsore!* 


3-  &&PQW  vm  WeSWcSw 

Ph.D.  Dissertation  /  Acoustics 


Acoustic  Emissions  from  Unsteady  Transitional  Boundary 
Layer  Flow  Structures 


6.  AUTHORS) 

R.  C.  Marboe,  G.  C.  Lauchle 


1.  PfSPOmmia  ORGANIZATION  NA&lg<S)  AN©  A0®8ES$<i§) 

Applied  Research  Laboratory 

The  Pennsylvania  State  University 

P.0.  Box  30 

State  College,  PA  16804-0030 


9.  SPONSORING  /  MONITORING  AGENCY  NAMI(S)  AND  AOOflgS^SS) 

Applied  Research  Laboratory 

The  Pennsylvania  State  University 

P.0.  Box  30 

State  College,  PA  16804-0030 


ii.  sumiMEimav  notes 


1 U.  ^$TmOTION/AVASiAi8U13f  STATEMENT 


Distribution  Unlimited,  Approved  for  Public  Release 


i®.  spommsm  i  mossst  o®m 

agency  report  mtm a 


13.  ABSTRACT  (Mdmmum  2CO  w&m) 

The  research  described  critically  examines  the  popular  approaches  to  modeling  the 
radiation  mechanisms  and  attempts  to  bring  some  degree  of  closure  to  the  physical  and  practical 
significance  of  noise  and  pseudo-noise  originating  in  the  laminar-to-turbulent  transition  zone 
within  a  natural  boundary  layer.  Several  refinements  were  made  to  theoretical  models  for 
directly  radiated  noise  based  upon  the  Liepmann  analogy  for  fluctuating  displacement  thickness 
including  the  incorporation  of  a  semi-empirically  derived  space-time  correlation  function  for  the 
intermittency  indicator. 

A  quiet  airflow  facility  was  developed  to  measure  the  direct  acoustic  radiation  from  a 
naturally  transitioning  boundary  layer.  Real-time  acoustic  intensity  measurement 
instrumentation  was  developed  if  measurements  of  isolated  spots  in  otherwise  laminar  flow  had 
been  necessary.  This  technique  uses  a  hot  film  signal  from  the  transition  structure  to  obtain  the 
coherent  output  intensity  (COI).  Model  predictions  are  compared  to  the  measured  acoustic 
radiation  from  a  naturally  transitioning  boundary  layer. 
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ABSTRACT 


The  true  acoustic  radiation  contribution  of  boundary  layer  flow  structures  has 
been  the  subject  of  debate  for  many  years.  The  research  described  here  critically 
examines  the  popular  approaches  to  modeling  the  radiation  mechanisms  and  attempts  to 
bring  some  degree  of  closure  to  the  physical  and  practical  significance  of  noise  and 
pseudo-noise  originating  in  the  laminar-to-turbulent  transition  zone  within  a  natural 
boundary  layer.  This  includes  improving  models  to  include  recent  computational  and 
experimental  statistics,  an  evaluation  of  model  sensitivities  to  input  parameters,  and  the 
applicability  significance  for  situations  of  engineering  relevance. 

There  have  been  significant  efforts  over  the  last  twenty  years  to  model  the 
statistics  of  the  wall  pressure  fluctuations  resulting  from  flow  structures  within  the 
laminar-to-turbulent  boundary  layer  transition  zone.  Models  for  the  directly  radiated 
noise  from  these  structures  have,  in  turn,  depended  upon  these  statistics.  Measurements 
of  the  space-time  statistics  for  these  pressures  allow  further  development  of  prediction 
codes.  Several  refinements  have  been  made  to  a  theoretical  model  for  the  directly 
radiated  noise  based  upon  the  Liepmann  analogy  for  fluctuating  displacement  thickness 
including  the  incorporation  of  a  semi-empirically  derived  space-time  correlation  function 
for  the  intermittency  indicator  and  a  dipole  source  description  for  the  induced  wall  normal 
velocity.  A  similar  model  exists  consisting  of  a  Lighthill  acoustic  analogy  using  a  two- 
fluids  model.  An  additional  model  for  radiation  by  vortex  structures  within  the  transition 
zone  is  qualitatively  reviewed.  The  predominant  effort  in  the  last  eight  years  has  been  to 
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apply  direct  numerical  simulation  methods  to  solving  for  acoustic  radiation.  These  efforts 
are  reviewed  to  help  define  their  useful  role  in  predicting  sound  radiation  from  transition. 
The  role  of  pressure  gradient  in  axisymmetric  body  flows,  flat  plate  flows,  and  over  air- 
and  hydrofoils  is  investigated  for  incorporation  into  the  model  through  a  sensitivity 
analysis  on  affected  flow  parameters. 

A  quiet  airflow  facility  was  developed  and  refined  to  measure  the  direct  acoustic 
radiation  from  a  naturally  transitioning  boundary  layer  using  time-averaged  acoustic 
intensity.  Real-time  acoustic  intensity  measurement  instrumentation  was  developed  in 
case  measurements  of  isolated  turbulent  spots  in  otherwise  laminar  flow  had  been 
necessary.  This  technique  uses  a  hot  film  signal  from  the  transition  structure  to  obtain  the 
coherent  output  intensity  (COI). 

Comparisons  are  made  of  model  predictions  to  the  measured  acoustic  radiation 
from  a  naturally  transitioning  boundary  layer.  Radiated  noise  measurements  isolating  the 
direct  radiation  from  a  transition  zone  demonstrated  similar  dependence  with  axial 
location  within  the  transition  zone  as  previous  wall  pressure  measurements.  This 
included  radiated  noise  levels  near  mid  transition  zone  that  were  higher  than  the  noise 
level  produced  by  a  fully  turbulent  boundary  layer  by  2-3  dB.  The  levels  measured 
suggest  that  the  radiation  from  transition  flow  structures  is  multipolar  and  has  a  higher 
radiation  efficiency  than  the  quadrupole  TBL  sources.  Transition  noise  per  unit  area  is 
greater  than  TBL  noise  per  unit  area.  Thus,  the  contribution  to  overall  directly  radiated 
flow  noise  from  the  transition  zone  in  typical  engineering  applications  is  negligible 
compared  to  the  radiation  from  the  much  larger  area  of  fully  turbulent  flow. 
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CHAPTER  1 


INTRODUCTION 


1.1  Background 

Much  debate  exists  over  the  true  acoustic  radiation  contribution  of  boundary  layer 
flow  structures.  Is  the  near-  and  far-field  noise  significant?  Current  knowledge  is  based 
on  analytical  theories  and  a  very  limited  database  of  measured  data  where  it  is  believed 
that  transition  noise  exists.  Unfortunately,  much  of  the  data  from  these  measurements  are 
questionable  due  to  possibility  of  contamination  by  other  noise  sources.  It  seems  that 
sonar  self-noise  is  affected  by  transition  and  that  a  radiation  noise  field  exists,  but  a 
refined  theoretical  or  computational  model  and  subsequent  experimental  verification  is 
required. 

The  research  described  in  this  dissertation  critically  examines  the  popular 
approaches  to  modeling  the  radiation  mechanisms  and  brings  a  degree  of  closure  to  the 
physical  and  practical  significance  of  noise  and  pseudo-noise  originating  in  the  laminar- 
to-turbulent  transition  zone  within  a  natural  boundary  layer.  This  effort  includes  updating 
of  models  to  include  computational  and  experimental  statistics  plus  an  evaluation  of 
model  sensitivities,  applicability,  and  significance  for  situations  of  engineering  relevance. 

There  have  been  significant  efforts  over  the  last  twenty  years  to  model  the  wall 
pressures  resulting  from  flow  structures  within  the  laminar-to-turbulent  boundary  layer 
transition  zone.  Knowledge  of  the  wavevector-frequency  spectrum  of  the  wall  pressures 
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can  be  used  to  predict  indirect  radiation  from  excitation  of  the  structure  over  which 
transition  occurs. 

Josserand  and  Lauchle  (1990)  showed  that  the  transition  zone  wall  pressure 
fluctuation  wavevector-frequency  spectrum,  near  the  convective  wavenumber  domain,  is 
equal  to  the  wavevector-frequency  spectrum  for  wall  pressure  fluctuations  under  a  fully 
turbulent  boundary  layer  (Figure  1.1)  multiplied  by  an  intermittency  factor,  y.  This  factor 
varies  monotonically  from  zero  at  the  beginning  of  transition  to  unity  for  the  fully 
turbulent  boundary  layer.  The  transition  zone  wavevector-frequency  spectrum  is: 

^riwuA ~  Y(V()  (1-1) 

Here,  Xx-(x-x0)/Ax  and  the  function  <PTBL  used  by  Josserand  and  Lauchle  (1990)  was 
developed  by  Chase  (1980).  However,  this  equation  is  not  valid  over  the  low 
wavenumber  or  sonic  region.  It  is  thus  the  intent  of  this  dissertation  to  refine  equation 
(1.1)  with  emphasis  on  the  prediction  of  the  direct  radiation  from  transition  flow 
structures  at  sonic  and  supersonic  wavenumbers. 

1.2  Boundary  Layer  Transitional  Flows 

Since  the  scope  of  this  dissertation  is  limited  to  the  modeling  of  direct  radiation 
from  the  transition  zone,  it  is  best  to  refrain  from  detailing  all  of  the  fundamental  work 
that  has  been  done  to  determine  the  fluid  mechanics  of  transition  mechanisms  and  the 
measurements  of  turbulent  flow  structures  which  form  the  basis  for  the  described  models. 
The  reader  is  directed  to  the  review  articles  of  Tani  (1969),  Gad-el-Hak  (1989),  and 
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Lauchle  (1991).  Some  of  the  more  important  prior  developments  relevant  to  the  present 
research  will  be  discussed  briefly  in  the  next  two  sections. 

1.2.1  Boundary  Layer  Flow  Structures 

Much  of  the  modeling  of  flow  structures  in  the  transitional  boundary  layer  has 
been  based  on  the  modeling  of  flow  structures  found  in  the  fully  developed  turbulent 
boundary  layer  (TBL).  The  characterization  of  these  structures  has  come  predominately 
through  experimental  measurements,  particularly  from  those  utilizing  laser  doppler 
velocimetry.  The  models  are  highly  empirical. 

The  description  of  flow  structures  and  characteristics  within  the  TBL  are  found  in 
Bershadskii  and  Tsinobar  (1994),  Honkan  and  Andreopoulos  (1999),  Klebanoff  et  al. 
(1992),  Koh  (1993),  Krogstad  and  Antonia  (1994),  Schlichting  (1979),  and  White  (1974). 
Wall  pressure  models  for  the  TBL  are  given  by  Corcos  (1963),  Chase  (1991),  Capone  and 
Lauchle  (1995),  Graham  (1997),  Keith  and  Abraham  (1997),  and  Singer  (1997).  Radiated 
noise  models  are  found  in  Powell  (1960),  Chase  (1987),  Barker  (1974),  Bergeron  (1973), 
Howe  (1998),  Laufferet  al.  (1964),  Pierce  (1981),  and  Tam  (1975). 

1.2.2  Transitional  Flows 

As  laminar  flow  over  a  flat  or  curved  surface  approaches  a  critical  Reynolds 
number,  hydrodynamic  instabilities  occur.  Examples  are  Tollmien-Schlichting,  Rayleigh, 
and  Kelvin-Helmholtz  for  different  flow  regimes.  [See  Huerre  and  Monkewitz  (1990), 
Panton,  (1984),  Schlichting  (1979),  or  White  (1974)  for  more  details.]  For  the  flat  plate, 
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two-dimensional  Tollmien-Schlichting  waves  form  in  the  boundary  layer.  These 
instabilities  become  non-linear  and  three-dimensional.  Bursting  occurs,  followed  by  the 
formation  of  Emmons  spots  (Emmons,  1951).  These  spots  are  turbulent  regions  within 
an  otherwise  laminar  flow  (Figures  1.2  and  1.3)  which  grow  and  coalesce  to  become  fully 
turbulent  flow.  The  transition  region,  of  length  Ax  beginning  at  x0,  is  usually  defined  as 
the  zone  in  which  distinct  Emmons  spots  exist. 

Therefore,  within  this  transition  region,  the  flow  can  be  thought  of  as  a  classical 
mixture  of  turbulent  and  laminar  flow  zones  with  concentrations  determined  by  the 
intermittency  function  y(tjx)  and  1  -y(rjJ,  respectively.  Here,  y(rjx)  is  the  time  average  value 
of  I(x,t),  the  intermittency  indicator  function  which  has  the  value  of  zero  for  a  laminar 
state  and  unity  for  a  turbulent  state.  The  system  of  coordinates  and  relevant  flow 
characteristics  are  shown  in  Figure  1.4.  This  function  is  typically  used  for  y=0.  For 
infinite  span,  we  expect  I(r{ to  be  statistically  homogeneous  in  the  rjz  direction  but 
non-homogeneous  in  the  rjx  direction.  So,  the  time  average  value  will  be  independent  of 
t]t  and  vary  from  0.0  at  t\x- 0  to  1.0  at  rjx=Ax.  The  burst  frequency,  N(rjJ,  describes  the 
mean  number  of  Emmons  spots  which  pass  a  point  per  unit  time.  The  general  behavior  of 
these  functions  through  the  transition  zone  can  be  seen  in  Figure  1.5.  While  the  ideal 
indicator  function  is  constructed  with  step  functions,  a  physically  realizable  indicator 
includes  rise  times  as  shown  in  Figure  1.6.  This  rise  time  accounts  for  the  time  for  the 
flow  to  change  from  a  laminar  flow  value  to  a  turbulent  flow  one  and  conversely 
(Lauchle,  1989).  As  shown  in  Figure  1.6,  the  rate  of  change  of  the  indicator  function  is 
described  by  e</,(.  If  this  change  in  flow  state,  with  inherent  increase  in  displacement 
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thickness  from  laminar  to  turbulent,  were  considered  as  a  piston  movement,  then  a  shorter 
time ,  t;,  would  result  in  higher  source  surface  velocity  and  a  higher  radiated  noise.  This 
time  should  be  a  function  of  the  flow  velocity  and  other  characteristics  of  transition  zone 
which  are  fluid  dependent. 

The  transition  zone  length  can  be  predicted  from  empirical  formulae  of  the  form 

=  C  Rex.  (1.2) 

where  the  Reynolds  numbers  are  based  on  the  subscripted  length  and  freestream  velocity, 
Chen  and  Thyson  (1971)  found  C=60  and  n=2/3  for  blunt  nosed  bodies  which 
accounted  for  positive  pressure  gradients  and  Mach  number.  However,  values  have  been 
derived  from  the  flat  plate  experimental  data  of  Dhawan  and  Narasimha  (1958):  C=15 
and  n= 0.8,  and  from  Josserand  and  Lauchle  (1990):  C=17  and  n=0.8.  Narasimha  et  al. 
(1984)  argued  that,  because  the  spots  do  not  propagate  at  consistent  speeds  proportional 
to  freestream  velocity,  every  pressure  gradient  results  in  a  unique  intermittency 
distribution,  and  Equation  (1.2)  must  be  adjusted  accordingly. 

Once  again,  there  have  been  many  investigations  covering  transition  mechanisms 
and  flow  structures  described  in  the  literature.  The  following  are  some  representative 
references  that  were  found  to  be  helpful.  Modeling  of  transition  mechanisms  is  described 
by  Barenblatt  (1990),  Gogineni.  and  Shih  (1997),  and  Stewart  and  Smith  (1992). 
Definitions  of  the  characteristics  of  turbulent  spots  themselves  can  be  found  in 
Barrow  et  al.  (1984),  Gad-el-Hak  and  Fazle  Hussain  (1986),  Henningson  and  Alfredsson 
(1987),  Jahanmari  et  al.  (1996),  and  Van  Atta  et  al.  (1982).  A  description  of  vortices  in 


5 


spots  is  available  in  Haidari  and  Smith  (1994).  Details  on  spot  intermittency  through  the 
transition  zone  is  provided  by  DeMetz  and  Casarella  (1973),  Farabee  et  al.  (1974),  Krane 
(1992),  Krane  and  Pauley  (1992),  Narasimha  (1957),  Schneider  (1995),  and  Zhang  et  al. 
(1995).  The  wall  pressures  resulting  from  turbulent  spot  formation,  growth,  and 
convection  are  well  documented  in  Bowles  and  Smith  (1995),  Ching  and  LaGraff  (1995), 
Gedney  (1979),  Gedney  and  Leehey  (1989),  Josserand  and  Lauchle  (1989),  and  Mautner 
(1982).  Studies  on  the  interaction  of  turbulent  spots  are  given  by  Cherny  and  Pauley 
(1994),  and  Coles  and  Savas  (1979).  The  growth  and  convection  of  turbulent  spots 
varies  with  the  pressure  gradient  over  the  transition  zone.  Several  investigations  have 
addressed  the  variations  with  both  favorable  and  adverse  pressure  gradients.  These 
include  Bandyopadhyay  and  Ahmed  (1993),  Blair  (1992  a,b),  Clark  et  al.  (1996),  Dey 
and  Narasimha  (1991),  Jeon  and  Kang  (1995),  Katz  et  al.  (1990a,b),  Ramesh  et  al. 

(1996),  Seifert  and  Wygnanski  (1995),  and  Walker  and  Gostelow  (1990). 

1.3  Goals  and  Scope  of  the  Present  Investigation 

This  investigation  was  begun  under  a  three  year  Office  of  Naval  Research 
Fundamental  Research  Initiatives  Program:  Hydroacoustics  of  Boundary  Layer 
Transition.  The  overall  objective  was  to  identify  the  physical  mechanisms  that  cause 
sound  radiation  from  intermittent  (particularly  transitional)  boundary  layer  flows  over 
rigid  surfaces.  Specifically  this  requires: 

1.  evaluation  of  existing  theoretical  and  semi-empirical  models  and  new 

developments  or  modifications  based  on  new  or  refined  space-time  statistics; 
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2.  experimental  determination  of  the  statistics  of  the  unsteady  fluid  mechanics  of 
these  flow  structures,  in  particular: 

-  develop  scalar  descriptors  of  turbulent  spot  growth, 

-  measure  space-time  development  of  transition  zone  fluctuating 
displacement  thickness  (for  individual  and  multiple  interacting  spots), 

-  relate  spot  growth  descriptors  to  sound  production  including  the 
characteristic  time  constants  associated  with  the  acoustic  radiation 
model; 

3.  obtain  unambiguous  measurements  of  the  acoustic  emissions  from  isolated 
boundary  layer  flow  structures  including  natural  transition  and/or  individual 
turbulent  spots;  and 

4.  use  the  existing  experimental  database  to  validate  the  theoretical  and 
semi-empirical  models. 

The  second  goal  was  the  basis  for  an  associated  investigation  which  is  reported  in 
the  Ph.D.  dissertation  by  Krane  (1992)  and  complements  the  measurements  by  Josserand 
(1986).  The  first,  third,  and  fourth  items  were  the  basis  for  the  present  investigation.  If 
this  research  were  viewed  from  a  verification  and  validation  process  (Coleman  and 
Steele,  1999)  perspective,  the  models  have  undergone  verification  to  assure  correct 
implementations  of  equations  and  address  computational  limitations  but  the  validation  is 
really  more  of  an  early  evaluation  of  the  applicability  of  the  may  models.  We  are  still  in 
the  discovery  phase  of  the  physical  mechanisms  which  are  contributing  to  the  transition 
zone  radiated  noise  which  has  been  experimentally  determined  for  the  first  time. 
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The  structure  of  this  dissertation  is  as  follows.  Important  boundary  layer 
parameters  and  definitions  are  covered  in  Chapter  1.  Five  approaches  to  modeling  the 
radiated  noise  from  a  transition  zone  are  covered  in  Chapter  2..  These  include:  a 
fluctuating  shear  stress  model  [Lauchle  1980],  a  fluctuating  displacement  thickness 
model  using  the  Liepmann  analogy  [Lauchle  1981,1989],  a  two-fluids  model  using  the 
Lighthill  analogy  [Somette  and  Lagier  1984a, b,  Lagier  and  Somette  1986,  and  Audet,  et 
al.  1989a,b],  direct  numerical  simulation  [Wang,  Lele,  and  Moin,  1996a],  and  a  vortex 
motions  model  [Kambe  1986  and  Hardin  1991].  Recent  activity  has  focused  on  the  use 
of  direct  numerical  simulation  (DNS)  using  full  Navier-Stokes  equations  solvers  to 
predict  statistical  quantities  relevant  to  a  Lighthill  approach.  The  problem  statement  then 
evolves  from  examining  some  predictive  comparisons  for  three  test  cases  for  which 
experimental  data  are  believed  to  be  influenced  by  transition.  These  include  a  flat  plate  in 
air,  an  axisymmetric  body  in  air,  and  an  axisymmetric  body  in  seawater.  It  is  apparent 
that  pressure  gradient  plays  an  important  role  in  the  pressure  fluctuations  and  its  effect 
will  also  be  discussed. 

The  computational  model  refinements  are  detailed  in  Chapter  3.  Based  on  the 
results  of  Josserand  and  Lauchle  (1990),  it  is  possible  to  make  refinements  to  Lauchle’s 
(1981)  model  by  incorporating  the  empirically  derived  space-time  correlation  function. 

In  addition,  several  other  refinements  can  be  made  using  measurements  by  Krane  and 
Pauley  (1993),  correlation  simplifications,  plus  accounting  for  retarded  time  differences 
and  vortical  motions.  Predictions  for  several  relevant  flow  conditions  are  presented  in 
Chapter  4. 
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It  is  important  for  any  model  to  have  experimental  validation.  The  direct  radiation 
from  transition  flow  structures  (distinct  from  measured  wall  pressures)  has  not  previously 
been  measured  in  isolation.  Chapter  5  contains  a  description  of  the  quiet  airflow  facility 
modified  for  these  measurements  and  the  flow-through  anechoic  chamber  of  the  Applied 
Research  Laboratory  at  The  Pennsylvania  State  University.  Acoustic  intensity 
instrumentation  and  signal  processing  was  employed  to  overcome  facility  noise.  Time- 
averaged  acoustic  intensity  measurement  techniques  were  used  successfully  for  the 
measurements  of  directly  radiated  noise  from  a  natural  transition  zone.  This 
measurement  includes  the  interactions  of  multiple  spots.  These  measurements  are 
described  in  Chapter  6.  In  anticipation  of  a  possible  need  for  additional  noise  rejection, 
real-time  acoustic  intensity  measurement  instrumentation  was  developed.  A  process  was 
developed  to  use  a  hot  film  signal  from  the  transition  flow  structure  to  obtain  coherent 
output  intensity  (COI).  This  technique  would  have  been  necessary  for  measurements  of 
the  radiated  noise  from  individual  spots  generated  and  convecting  in  otherwise  laminar 
flow.  Conclusions  and  recommendations  for  further  investigation  are  provided  in 
Chapter  7. 

Mathcad  ®  software  from  Mathsoft,  Inc  was  used  for  the  computational  modeling. 
Printouts  of  the  worksheets  are  provided  in  Appendix  A.  A  chronology  of  the 
development  of  the  Quiet  Airflow  Facility  is  detailed  in  Appendix  B.  Finally,  details  of 
the  electronics  developed  for  the  real-time  acoustic  intensity  measurement  are  given  in 
Appendix  C. 
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Figure  1 . 1  General  features  of  the  wavevector-frequency  spectrum  of  the  wall 

pressure  fluctuations  under  a  turbulent  boundary  layer  where  kz=0  and  kx  is 

the  independent  variable  for  constant  frequency. 


Figure  1.2  Photograph  of  flow  on  a  heated  body  showing  transition  flow  structures 
[from  Lauchle  1991]. 


11 


Plan  View 


Elevation  View  on  Centerline 

y 


Figure  1.3  Schematic  views  of  a  turbulent  spot  and  its  growth  with  downstream 
distance  [from  several  sources]. 
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Laminar  Transition  Turbulent 

Boundary  Layer  Boundary  Layer 


ZX  =  T]X/Ax 


Figure  1.4  System  of  coordinates. 
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(a) 


Figure  1.6  The  modified  indicator  function  7(0  is  shown  in  (a)  and  its  first  time 
derivative  in  (b)  [from  Lauchle  1981]. 
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CHAPTER  2 


MODELS  FOR  BOUNDARY  LAYER  TRANSITION 
ZONE  ACOUSTIC  RADIATION 

The  modeling  of  flow-induced  noise  from  boundary  layer  flows  has  been  well 
researched  for  many  years.  Several  papers  are  significant  and  form  the  basis  for 
investigation  of  transition  noise.  The  works  of  Lighthill  (1952, 1954),  Ffowcs  Williams 
(1967, 1969, 1977, 1982),  Curie  (1955),  Crighton  (1975),  and  Blake  (1986),  and  Howe 
(1998)  provide  a  good  general  background. 

Several  approaches  have  been  taken  to  model  the  unsteady  pressures  due  to 
boundary  layer  transition  processes.  These  have  predominately  used  the  Lighthill  analogy 
as  a  basis.  Lauchle  (1981)  took  a  more  phenomenological  approach  by  modeling  the 
displacement  thickness  change  due  to  transition  as  a  piston  motion  with  non-uniform 
velocity  distribution.  This  was  later  supported  through  a  Lighthill  approach  by  Somette 
and  Lagier  (1984a).  Audet  et  al.  (1989a)  then  claimed  the  experimental  appearance  of 
monopole-like  acoustic  sources.  While  the  gross  fluctuations  of  the  growth  in 
displacement  thickness  seem  to  be  modeled  well  by  the  piston-like  sources,  there  is  also 
much  interest  in  modeling  the  acoustic  sources  due  to  the  vortical  structures  of  the  spots. 
In  this  chapter,  four  of  these  models  will  be  presented.  Direct  numerical  simulation 
efforts  will  be  discussed  briefly. 
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2.1 


Fluctuating  Shear  Stress  Using  a  Ligh thill  Analogy 


Lauchle  (1980)  used  the  Lighthill  acoustic  analogy  to  find  dipole,  quadrupole,  and 
octupole  components  of  radiated  noise  from  transition.  He  assumed  radiation  efficiencies 
scale  on  Mach  number  according  to  M3  and  Ms  for  the  dipole  and  quadrupole  sources 
respectively.  Since  a  low  Mach  number  flow  was  also  assumed,  he  neglected  the 
quadruples  and  octupoles  due  to  the  fluctuating  Reynolds  stresses  and  their  images. 
However,  their  large  amplitudes,  relative  to  those  of  the  wall  shear  stresses  that  form  the 
dipole  contribution  may  compensate  for  the  lower  radiation  efficiency  and  make  this  a 
poor  assumption. 

The  Lighthill  equation  for  sound  radiation  from  fluid  dynamic  sources  is: 


=  J!Ll 

dt 2  a*,2  dxfy 


(2.1) 


where  Ti}  is  Lighthill’s  stress  tensor  given  by: 

Tu  =  p¥ro^(pVc\  (2.2) 

It  should  be  noted  that  Fraudelplonika  (2000)  suggests  that  this  stress  has  an  additional 
contribution,  but  that  was  not  considered  in  the  context  of  this  development.  The 

tensor  of  Equation  (2.2)  consists  of  the  fluctuating  Reynolds  stress  tensor,  puuj ,  the 
fluctuating  viscous  shear  stress  tensor,  and  a  “correction”  term  which  is  essentially 
zero  for  hydrodynamic  flows. 
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Employing  Curie's  equation  and  Powell's  (1960)  extension  for  image  flow  with  a 
rigid,  planar  surface,  Lauchle  derived  an  estimate  of  the  acoustic  pressure  due  to  shear 
stress  fluctuations  only: 


P(r,t)  * 


COS0 


2  7t  r  co  dt 


T  ff  xo  (V'“)  dS(\)>  M  «  1 

at  J  J  c 


(2.3) 


where  0  is  the  angle  from  the  x  axis  to  the  observation  point.  This,  of  course,  assumes 
that  the  contribution  of  the  longitudinal  and  lateral  quadrupoles  is  insignificant  based  on  a 
M  ~2  reduction  in  radiation  efficiency  for  low  Mach  number  flows. 

Using  intermittency  to  model  the  wall  shear  stress  in  a  mix  of  laminar  and 
turbulent  states,  dzjdt  reduces  to  (xT-xL)  dl/dt  =  o(x;)  i(x.,t)  .  This  results  in  an 

equation  for  the  power  spectral  density  of  the  sound  radiated  per  unit  spanwise  width  of 
transition: 


dG(r,  co)  (o2cos2da2(x0) 

3*  '  4  rt2  rJ  c*  JIJ  * 


*X4’t) 


eim  dr]x  d$x  dx  . 


(2.4) 


This  assumes  that  the  space-time  correlation  function  for  the  indicator  can  be  separated, 
the  same  as  Corcos  (1963)  did  for  the  homogeneous,  fully  developed  turbulent  boundary 
layer.  This  requires  that: 

(7/0  =  Rx(r\x£x,x)Rz(\£z)  .  (2.5) 


18 


It  is  assumed  that  the  indicator  function,  i(xyt),  is  stationary  in  time  and  homogeneous  in 
x3.  A  longitudinal  space-time  correlation  function  was  then  formulated  based  on  physical 
notions  of  the  transition  process.  The  details  will  not  be  presented  here. 

Incorporation  of  this  correlation,  which  is  itself  a  function  of  the  intermittency, 
spot  burst  frequency  rate,  and  a  flare  parameter  ( a  )  (defining  the  decay  envelope),  yields 
a  final  form  of  Equation  (2.4)  reduced  to  a  single  integration: 


dGjrfi o)  =  cos29  g\xo)  uo  (Ax)2 
2  7T2  r2  cl 


F(kc Ax,  a* Ax, 


(2.6) 


with 


u  '  1  272  7  2  n  -  e *  4185  z\  e  ~  4185  z* 

F(k  Ax,  a’ Ax,  =S.)  =  (kc Ax)2  f  T - 1  272  z  il— £ IS - - 

uo  {  [  6.472  (ujuf  Z2  e ' 837  z'  +  (kcAx)2  j  [  (a  Ax)2  +  (kcAx)2  J72 


2  a  Ax 


[  [  (a ’Ax)2  +  (kcAx)2  I72 


-  e'a  iuZ  sin(it  Ajc  Z  -  <p)  +  e o-z>  sjn[£ Ax  (i-Z)  -  <p] 


dZ 


(2.7) 


where  ^  =  tan1  (a  Ax/ kgAx),  kc=co/uc  and  Z  =  rj/Ax .  Based  on  limited  experimental 
data,  this  model  appears  to  significantly  under-predict  the  transition  zone  radiated  noise. 
This  led  to  a  search  for  a  more  efficient  mechanism  which  would  produce  the  higher 
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measured  levels.  It  does  not  appear  to  be  advantageous  to  pursue  any  further 
modifications  to  this  model  that  uses  only  the  wall  shear  stresses  as  the  source  term. 

2.2  Fluctuating  Boundary  Layer  Displacement  Thickness  Using  a 

Liepmann  Analogy 

In  order  to  account  for  the  previously  neglected  terms  for  the  fluctuating 
velocities,  Lauchle  (1981,1989)  used  the  Liepmann  (1954)  analogy  rather  than  calculating 
the  Reynolds  stresses  in  transition  as  required  by  the  Lighthill  equation. 

The  Liepmann  analogy  uses  the  fluctuations  in  the  boundary  layer  displacement 
thickness  to  describe  the  normal  velocity  of  a  vibrating  surface.  This  is  the  “surface”  that 
the  potential  outer  flow  feels  as  it  is  displaced  by  the  boundary  layer.  A  wave  equation  is 
then  solved  for  the  acoustic  radiation  of  this  vibrating  surface.  The  results  of  this 
approach  differ  from  Lighthill’ s  analogy  only  in  that  the  unknown  source  is  described  in 
the  boundary  condition  and  not  in  the  source  terms  of  the  Lighthill  inhomogeneous  wave 
Equation  (2.1).  The  use  of  displacement  thickness  fluctuations  in  acoustic  modeling  was 
also  recognized  by  Lighthill  (1958)  and  Howe  (1981). 

2.2.1  Formulation  Based  on  Lauchle  Theoretical  Space-Time  Correlation 

Function 

A  solution  for  the  Helmholtz  equation  is  obtained  using  the  Green’s  function,  G0 
for  an  infinite  half  space  with  a  Neumann  boundary  condition  (essentially  a  baffled 
"piston"  with  velocity 
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(2.8) 


G0(xy,t,x )  dS 


with  boundary  conditions  of 


dp, 

dn 


=  P 


dt  ’ 


(2.9) 


/  .x  ,3<Dm 
VhM  =  (-^)l,=0 


dt 


7-  A8’ 

dt  z 


(2.10) 


where  n  represents  the  direction  normal  to  the  surface,  v„  is  the  fluctuating  normal 
velocity  and 


8‘Ovv)  =  [i  ~/(iyn »/)]  b[(x)  +  i(\,\,t)  bfa). 


(2.11) 


A  power  spectral  density  formula  was  derived  as 


/  7  /_! 

8  it2  [l+(a>(,.)2]  J.  {  i  U- 


7 1  l*-0’+4l 


(2.12) 


With  an  estimated  space-time  correlation  function  for  the  intermittency  indicator,  the  far- 
field  solution  for  the  power  spectrum  of  the  radiated  pressure  per  spanwise  width 


becomes 


**•»>  -  P  <A8->  “-A  F-a  tx.  a'&x,  i) 


dz 


8  7C2  r2 


Ax  M, 


(2.13) 
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where 


F*  = 


(k  Ax)2 


[1  +  (k  Ax)2  (-^i)2] 
Ax 


F(k  Ax,  a* Ax,  — ) 

W„ 


(2.14) 


with  F(klAx,  a  Ax,  u/u0)  given  by  Equation  (2.7). 

More  generally,  the  normal  velocity  in  Equation  (2.10)  is  created  both  by  flow 
over  the  spot  and  the  spot’s  growth  itself.  Wygnanski,  et  al.  (1976)  found  that  the  normal 
velocity  is  directed  toward  the  plate  at  the  spot  leading  interface  and  outward  over  the  rest 
of  the  spot.  This  is  shown  in  Figure  2.1. 

The  opportunity  exists  to  make  several  modifications  to  this  model  based  on 
recent  experimental  results.  These  are  described  in  general  in  Lauchle  (1991),  Marboe 
and  Lauchle  (1992, 1994)  and  in  detail  in  Chapter  3. 


2.2.2  Formulation  Based  on  Josserand-Lauchle  Measured  Space-Time 

Correlation  Function 

Theoretical  forms  for  the  space-time  correlation  function  were  the  only  functions 
that  could  be  employed  prior  to  the  actual  measurements  of  these  statistics  in  the 
transition  zone  by  Josserand  (1986).  This  is  now  suggested  for  substitution  into  Equation 
(2.12). 

The  development  of  a  semi-empirical  space-time  correlation  function  was 
performed  by  Josserand  (1986)  using  a  zero  pressure  gradient,  natural  transition  zone  on  a 
flat  plate  in  the  ARL  Penn  State  1.22  m  wind  tunnel.  From  Josserand  and  Lauchle 
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(1990): 


(l-r.)  expj-  5  (4  ♦  200  |t-±|)  |t-^|  }  •  expj^J^L 


|i4i 


•  exrt  - 


1260  |t-— | 


15.75A*  - 


1  +  71  |  A 


> 

IM 

•  exp 

*  IM  1 

1  ♦  14.2  |Aj  J 

0.0014  +  |  AJ  ,  +  1300k-— I2 

uc 

(2.15) 


where  A  =  -ln[l-exp(-4.27/Ax)],  (units  of  Ax  are  in  meters), 

yu  =  y(rjx)  and  yd  =  y(rjx+£J  when  the  reference  location  is  upstream  of  the  second 
location, 

yu  =  y(r}x+£J  and  yd  =  yfrjJ  when  the  reference  location  is  downstream  of  the 
second  location. 


2.2.3  Chase  Statistical  Modification 

Chase  (1989)  has  suggested  that  the  first  exponent  term  in  Equation  (2.15)  can  be 
replaced  by  a  more  general  term  based  on  the  probability  of  state  transition  to  turbulence 
of  form: 


exp  -  — 4— |t-^|  •  (2.16) 

However,  Josserand  found  that  both  terms  predict  the  measured  correlation  values  well 
for  0.3 <;y5;0.7;  but,  that  outside  of  these  limits,  the  Chase  theoretical  term  tended  to  over¬ 
predict  the  measured  correlation.  Both  terms  are  available  as  options  in  the 
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computational  codes  developed  for  this  dissertation. 


2.3  Two-Fluids  Model  Using  a  Lighthill  Analogy 

The  two-fluids  model  was  first  proposed  by  Somette  and  Lagier  (1984a,b),  further 
developed  with  Lagier  and  Somette  (1986),  and  applied  by  Audet  et  al.  (1989a, b).  A 
Lighthill  approach  was  used  to  validate  Lauchle’s  phenomenological  approach.  The 
acoustic  radiation  described  by  the  two-fluids  (states)  model  is  also  analogous  to  a 
fluctuating  piston  powered  by:  (1)  fluctuations  due  to  turbulence  (vorticity)  below  the 
interface,  and  (2)  "growth"  processes  due  to  the  laminar-turbulent  transition  (accounted 
for  by  Lauchle  (1981)).  The  normal  quadrupole  terms  are  found  plus  one  monopole  due 
to  entropy  production  by  viscous  dissipation  (dependent  upon  the  time  derivatives  of  the 
intermittency  indicator  function)  and  another  monopole  due  to  kinetic  energy  dissipation 
(dependent  upon  the  spatial  derivatives  of  the  indicator  function). 

Following  the  Lighthill  analogy,  the  fluctuating  density  is  found  to  be  p '  =ypT  +  <P 
with“acoustic  potential”,  <P 


a2o 

dt2 


-  cl  V2  <D  =  s0  +  5 


and  the  two  source  terms,  S0  and  S,  given  by 


(2.17) 


so  =  co  (2  VrVpr  +  pr  v2y)  - 


Ar  fPr) 

Vt  dt 2  dt  dt 


(2.18) 
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S(r,t)  =  p. 


d2VVj 


dxi  dXj 


d\W  d2vLvL 

■  '  ‘  -  (1-7)  1 


dxj  dXj 


dxt  dXj 


-  2  i3 L  -L(xt-tl)  -  dy 

dxt  dXj,J  ,J  dx,  dXj  V  11  ,j) 


(2.19) 


a2©,. 


+  2  Vy*V(p  r-p  L)  +  ip  T-p  L)  V2y) 


ax; 


where  v(>  is  the  fluid  particle  velocity  component,  Ttj  is  the  fluctuating  viscous  stress 
tensor,  and  with  Gy  defined  by  (using  shear  viscosity,  t],  and  volume  viscosity,  0 


dy 


•  i  .,',1 


*  (v  r-v  L)  ■  Vy 


(2.20) 


(In  this  notation,  y(x,t)  =  l(x,t) ).  The  S0  term  comes  from  the  wall  pressure  fluctuations 
that  are  due  to  a  fully  turbulent  flow  and  whose  amplitude  is  weighted  by  the 
intermittency  indicator.  It  contains  monopole,  dipole,  and  quadrupole  terms.  The  second 
term,  S,  represents  the  incompressible  part  of  the  velocity  fluctuations  in  the  transitioning 
boundary  layer.  These  are  monopole  sources.  The  ratio  of  SJS  is  proportional  to  Mach 
number,  so  S0  could  be  neglected  in  hydrodynamic  cases. 

Audet  et  al.  (1989a, b)  developed  a  model  for  the  wall  pressure  spectrum  using  the 
two  fluids  theory  to  get  a  fluctuating  pressure  field  of  ypT  +  /V™,,,.  accounting  for  the 
contribution  of  the  turbulent  wall  pressure  modified  by  the  intermittency  and  the  radiated 
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acoustic  field  (calculated  using  the  Lauchle  (1981)  theoretical  intermittency  indicator 
correlation).  This  acoustic  pressure  is  given  by 


P 


acoustic 


2-  f  — —  A5*  - - -  dS. 

*  {  \r-r0 \  dt2 


(2.21) 


An  example  of  this  calculation  is  shown  in  Figure  2.2.  Equation  (2.21)  has  the  advantage 
over  Equation  (2.13)  (Lauchle  1981)  of  predicting  radiation  as  a  function  of  rj/dx  because 
of  accounting  for  retarded  time  differences  in  the  surface  integral  for  vn(x,t).  However, 
Audet,  et  al.  (1989a, b)  used  simplifying  assumptions  for  their  calculations  such  as 
assuming  a  far  field  receiver  and  separation  of  spatial  dependence  from  temporal 
dependence  which  limits  the  usefulness  of  this  approach. 

Certainly  the  two-fluids  model  pioneered  by  Lagier  and  Somette  provides  a  more 
rigorous  model  of  possible  transition-induced  pressure  fluctuations,  which  contribute  to 
both  the  nearfield  and  the  farfield.  However,  using  this  method  to  get  an  estimate  of 
acoustic  radiated  pressures  becomes  an  almost  intractable  problem  for  the  engineer.  This 
would  seem  to  be  better  suited  to  use  with  direct  numerical  simulation  and  staggering 
amounts  of  computer  time  and  memory.  When  the  sufficient  simplifications  are  made  to 
allow  a  calculable  estimate,  it  provides  little  benefit  when  compared  to  the  fluctuating 
displacement  thickness  approach. 
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2.4 


Vortex  Motions 


Kambe  and  Minota  (1981)  and  later  Kambe  (1984,  1986)  used  the  viscosity  terms 
for  vortices  to  predict  quadrupole  and  low  amplitude  monopole-like  radiation  from  vortex 
motion  particularly  vortex -edge  interaction.  According  to  Kambe  (1986),  the  transition 
process  evolves  from  Tollmien-Schlichting  waves  which  are  amplified  and  become 
associated  with  the  concentration  of  vorticity  along  discrete  lines.  These  lines  distort  into 
vortex  loops  which  continue  to  distort  and  grow  as  turbulent  spots  until  coalescence.  The 
stresses  at  the  vortex  interface  become  a  significant  sound  source. 

Hardin  (1991)  argues  that,  by  a  process  of  elimination,  the  major  source  of  sound 
in  the  turbulent  boundary  layer  is  the  formation  of  horseshoe  vortex  structures  by  the 
interaction  of  the  streamwise  vorticity  with  the  primary  spanwise  vorticity.  This  may  be 
an  important  source  of  sound,  in  addition  to  the  normal  velocity  of  the  fluctuating 
displacement  thickness. 

There  have  been  many  advocates  of  vortex  noise  generation  modeling.  The 
leaders  have  been  Powell  (1964, 1995),  Obermeier  (1977, 1979, 1980, 1985),  Mohring 
(1978),  Howe  (1998a, b,  1999a, b),  and  Mohring  et  al.  (1983).  The  dynamic  response  of  a 
vortex  to  interaction  with  a  sound  field  [Lund,  1989],  a  surface  [Kambe,  Minota,  and 
Ikushima,  1985],  and  another  vortex  [Kambe  and  Minota,  1983]  have  all  been  measured 
experimentally.  Two  additional  experiments  have  yielded  good  qualitative  data  for 
vortex  formation  behind  a  vortex  generator  [Littell  and  Eaton,  1991]  and  in  boundary 
layer  transition  [Williams,  Fasel,  and  Hama,  1984] . 

While  some  aspects  of  vortex  motions  are  incorporated  in  the  fluctuating 
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displacement  thickness  and  two-fluids  models,  they  probably  have  not  been  fully 
accounted  for  in  the  computational  approximations  for  the  normal  velocity  and  therefore 
warrant  being  highlighted  as  an  additional  modeling  possibility. 

The  contribution  of  vortex  sound  could  also  play  more  of  a  role  in  flows  with  non¬ 
zero  pressure  gradients  due  to  growth  and  stretching. 

2.5  Direct  Numerical  Simulation 

The  state  of  the  art  in  computational  fluid  dynamics  application  to  solution  of  the 
acoustic  problem  in  fluid  flow  can  be  best  summarized  with  a  quote  from  Crighton 
(1975): 

“Lack  of  knowledge  of  features  of  the  flow  does  not  entitle  us  to  make 
assumptions  and  approximations  except  with  the  greatest  care.  We  are 
dealing  with  a  small  by-product,  and  approximations  which  may  be  quite 
safe  as  far  as  the  flow  dynamics  are  concerned  may  be  fatal  for  the  by¬ 
product.  Vast  errors  may  easily  be  incurred  in  this  way  through 
apparently  harmless  approximations,  especially  in  underwater  flows, 
where  the  Mach  numbers  are  always  extremely  small,  and  the  acoustic 
energy  is  a  correspondingly  minute  fraction  of  the  total  energy....  Highly 
accurate  numerical  predictions  are  simply  not  to  be  expected  -  a  reflection 
of  the  fact  that  noise  levels  usually  vary  far  more  widely  than  do  the 
parameters  which  describe  the  noise-producing  flow.  ” 

Having  read  that,  it  is  important  to  state  that  significant  progress  has  been  made  in 
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the  last  eight  years.  Some  of  the  best  work  has  been  done  by  Wang,  Lele,  and  Moin 
(1994, 1996a,  1996b).  Their  modeling  focuses  on  the  quadrupole  source  functions  in  the 
solution  of  the  full  incompressible  Navier-Stokes  equations.  They  find  that  negligible 
sound  is  created  during  primary  and  secondary  instability  stages  but  that  disintegration  of 
the  detached  high  shear  layer  and  vortex  shedding  near  the  boundary  layer  edge  results  in 
amplification  of  the  Reynolds  stress  quadrupoles,  particularly  at  the  frequency  of  vortex 
shedding.  Their  conclusion  is  that  the  acoustic  quadrupoles  and  images  dominate  up 
through  the  formation  of  spots  and  then  in-plane  dipoles,  due  to  surface  shear  stress, 
dominate  in  late  transition.  [Note  that  this  is  the  modeling  approach  of  Lauchle  (1980) 
discussed  in  section  2.1].  The  small  scale  features  and  large  convected  eddies  are 
considered  to  be  inefficient  acoustic  sources  despite  larger  amounts  of  energy.  They  also 
show  that  small  structures  traveled  faster  than  the  lambda  vortex  which  elongates  the 
disturbance  region  yet  neighboring  spot  interaction  is  negligible.  This  group  has  also 
developed  a  hybrid  technique  for  low  Mach  numbers  that  uses  a  Navier-Stokes  solver  for 
incompressible  hydrodynamic  terms  and  a  Lighthill  analogy  for  the  acoustic  contribution. 
Notably,  an  exit  boundary  correction  is  also  been  developed  for  the  quadrupole  source 
calculated  in  the  truncated  domain,  V0 

e,/»  -  -f2  fT,/ySd\  ♦  |  /  VJ,fyJ)d*y  (2.22) 

or  w  c 

Vo  ^ O 

given  by  the  second  term  which  is  equal  to  the  time  derivative  of  the  Lighthill  stress 
fluxes  carried  by  convecting  eddies  across  the  exit  surface.  In  contrast,  Kloker, 
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Konzelmann,  and  Fasel  (1993)  use  a  “re-laminarization”  zone  downstream  of  transition 
before  passing  through  the  outflow  boundary  to  suppress  reflected  flow  disturbances  in 
direct  numerical  simulation  (DNS). 

A  very  important  part  of  any  DNS  modeling  effort  is  the  introduction  of  a 
disturbance  to  cause  breakdown  to  turbulence.  An  excellent  overview  is  provided  in  Rist 
and  Fasel  (1995).  They  discuss  the  use  of  temporal  vs.  spatial  models.  They  recommend 
that  most  “real  life”  cases  should  only  use  a  spatial  model,  but  the  penalty  is  in  processing 
resources.  They  also  discuss  a  mean  flow  which  is  obtained  from  a  time  averaged 
disturbance  flow.  This  quantity  is  like  what  Krane  (1992)  measured.  This  is  contrasted 
with  a  steady  base  flow  that  is  obtained  from  numerical  solution  of  the  Navier-Stokes 
equations  before  the  unsteady  flow  field  is  calculated.  They  also  found  a  high  shear  layer 
structure  that  sat  on  top  of  the  lambda  vortex  structure  and  instantaneous  maps  of  normal 
velocity  (see  Figure  2.3)  that  show  a  pushing  of  the  flow  structures  away  from  the  wall. 

A  number  of  computational  studies  have  achieved  results  which  may  contribute  to 
our  qualitative  understanding  of  transitional  flow  structures.  Coherent  structures,  similar 
to  bursts  in  the  fully  developed  TBL,  were  calculated  by  Rempfer  and  Fasel  (1994). 

There  is  a  spike  stage  in  their  result  which  yields  fluctuating  velocities  that  are  significant 
for  acoustic  modeling.  Ducros  et  al.  (1996)  performed  a  large  eddy  simulation  (LES)  of 
boundary  layer  transition.  This  took  80  hours  (CRAY  YMP)  using  LES  (compared  to 
800  hours  using  DNS,  Rai  and  Moin  (1993))  and  a  filtered  structure  function  on  the 
velocity  field  to  eliminate  large  scale  oblique  perturbations.  Their  results  show 
qualitative  progress  and  could  aid  in  interpretation  of  the  role  of  vorticity  in  the  acoustic 
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signature.  But  obviously,  this  is  not  a  time  effective  means  for  prediction  or  design. 
Yang,  et  al.  (1992)  performed  computational  analyses  for  natural  transition  in  an  adverse 
pressure  gradient.  Similarly,  the  works  of  Goldstein  et  al.  (1992),  Goldstein  and  Leib 
(1993),  plus  Henningson  and  Kim  (1991)  show  results  which  provide  some  insight  into 
the  roles  of  vorticity,  local  separations,  and  spot  growth.  Wilde  and  Rose  (1997) 
concluded  from  their  use  of  a  temporal  model  of  K-type  transition  and  spatial  model  of 
bypass  transition  that  the  generated  sound  comes  from  the  streamwise  dipole  of  the  wall 
shear  stress. 

The  work  of  Joslin  (1991)  on  the  modeling  of  the  effects  of  compliant  walls  on 
primary  and  secondary  instabilities  in  boundary  layer  transition  warrants  further 
investigation  relative  to  efficient  numerical  simulation  compared  to  DNS.  This  may  also 
have  some  applicability  in  further  refining  the  fluctuating  displacement  thickness  model. 

Abraham  and  Keith  (1997)  have  shown  that  DNS  simulations  for  the  fully 
developed  TBL  can  produce  a  reasonable  estimate  for  wall  pressure  events  associated 
with  sweeps  and  bursts.  They  also  show  that  the  number  of  positive  and  negative  events 
occur  in  comparable  densities  and  that  their  average  spectra  are  lower  in  level  by  3  to  7 
dB  than  their  peak  spectral  amplitudes.  This  is  significant  in  light  of  the  positive  and 
negative  normal  velocity  distributions  sensed  for  the  passage  of  for  turbulent  spots 
described  in  Krane  and  Pauley  (1995). 
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2.6  Direct  and  Indirect  Radiation  From  Complex  Geometries 

The  research  described  here  is  intended  to  critically  examine  the  popular 
approaches  to  modeling  the  radiation  mechanisms  and  bring  some  degree  of  closure  to  the 
physical  and  practical  significance  of  noise  and  pseudo-noise  originating  in  the  laminar- 
to-  turbulent  transition  zone  within  a  natural  boundary  layer.  This  effort  includes 
updating  of  models  to  include  recent  computational  and  experimental  statistics  plus  an 
evaluation  of  model  sensitivities,  applicability,  and  significance  for  situations  of 
engineering  relevance  such  as  for  the  noses  of  bodies  of  revolution  and  for  airfoils. 

Although  little  has  yet  been  done  to  model  transition  noise  for  transition  zones 
specifically  on  curved  surfaces,  Howe  (1981)  has  made  some  generalizations.  He 
contends  that  displacement  thickness  fluctuations  will  give  efficient  sound  radiation  in 
the  regions  of  a  surface  or  body  where  the  local  radius  of  curvature  is  smaller  than  the 
acoustic  wavelength.  For  an  airfoil,  this  would  be  in  the  leading  and  trailing  edge  regions 
under  compact  conditions. 

2.6.1  Vicinity  of  a  Half-Plane  Leading  Edge 

The  effect  of  scattering  of  near-field  energy  due  to  the  surface  geometry  near  the 
leading  edge  or  trailing  edge  of  a  plane  is  to  enhance  the  shear  stress  dipole  sources  in  the 
transition  region  based  on  formulations  by  Howe  (1999a,c).  These  sources  may  be  less 
important  for  flat  walls  compared  to  the  quadrupole  sources  in  the  turbulent  spots  but,  for 
a  curved  leading  edge,  tend  to  contribute  to  the  low  frequency  signature  following  a  U6 
velocity  scaling.  Two  factors  play  into  this.  First,  due  to  the  curvature  of  the  leading 
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edge,  the  transition  tends  to  occur  in  a  region  of  adverse  pressure  gradient  which  will 
accelerate  the  transition.  Second,  the  nearfield  of  the  transition  region  hydrodynamic 
unsteadiness  will  extend  forward  to  the  leading  edge  and  diffract  around  the  curved 
surface. 

Howe  models  the  transition  region  with  an  ‘upwash’  velocity  on  the  body  surface 
using  the  Biot-Savart  formula  giving  the  (time  dependent)  distributed  vorticity  induced 
instantaneous  velocity  field  as  though  the  surface  was  not  present.  He  uses  the 
proportional  local  blocked  surface  pressure  to  represent  the  spectrum  of  the  upwash 
velocity.  The  fully  developed  TBL  spectrum  is  used  with  the  addition  of  an  intermittency 
weighting  factor. 

The  assumptions  for  this  model  include: 

•  very  low  Mach  number,  2-D  flow; 

•  symmetric,  finite  thickness  plane  (better  termed  a  hydrofoil); 

•  hydrofoil  chord  is  acoustically  non-compact; 

•  hydrofoil  span  is  much  greater  than  the  boundary  layer  thickness; 

•  hydrofoil  thickness  is  acoustically  compact,  and  that  the  leading  edge  is  within  the 
nearfield  of  the  transition  zone; 

•  zero  angle  of  attack  and  boundary  layer  fluctuations  on  opposite  sides  of  the  foil 
are  equal  and  independent  so  an  estimate  of  one  side  can  just  be  increased  by  3 
dB; 

•  the  integration  of  the  vorticity  is  restricted  to  the  non-linear  region  of  the 
boundary  layer  above  the  viscous  sublayer,  so  bound  vorticity  on  the  hydrofoil  is 
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ignored;  and 

•  high  Reynolds  number,  the  contribution  of  tangential  shear  stresses  is  ignored. 

Using  the  Green’s  function  for  a  rounded  leading  edge,  Howe  approximates  the  farfield 
acoustic  pressure  as 

p  to  JIF0  sin1/2\|/  sin(— )  e'*01*1 

p(x, co)  =  - - -  <f  0*(y)  v  (y,co)  dS(y)  ,  |x|=»°°  (2*23> 

%s[2k\x\  Js 


where 

0,  \|/  define  the  observer  orientation  direction  x  /  |x| 

0*(y)  is  the  velocity  potential  describing  the  flow  of  unit  speed  in  a 
clockwise  sense  over  the  profile  S  of  the  hydrofoil 
vto  is  the  normal  component  of  the  upwash  velocity. 


When  a  conformal  mapping  is  done  for  the  hydrofoil  surface,  equation  (2.23)  is 
approximated  as: 


UJc 

<D(x,oo)  *  0.35 — - 


n 


-  (— — r)( 1  )  sin  v|/  sin2(^)  |F(co,/:  )|2)  ,  |x| 

\r\2  1+a  2 


(2.24) 


where 
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(2.25) 


F(a>,kc)  =  f  [1  +  e'^h'(s)]  e'V  s=s(Q,  s0=s&0) 

a,  P  define  the  nose  shape  in  the  conformal  mapping  of  the  profile  S 

onto  the  real  C,  axis; 

h  is  the  constant  thickness  of  the  hydrofoil; 

L  is  the  span  of  the  hydrofoil;  and 

d>pp(s,cfl)  is  the  point  pressure  spectrum  of  TBL  wall  pressure  fluctuations 
(To  account  for  the  transition  zone,  Opp(s,to)  includes  a 
weighting  with  y(s).) 

Equation  (2.24)  represents  the  far  field  sound  due  to  turbulent  hydrodynamic  fluctuations 
on  the  upper  surface  of  the  hydrofoil. 

2.6.2  Bodies  of  Revolution 

We  know  that  transition  in  an  adverse  pressure  gradient  will  result  in  higher 
radiated  and  structural  excitation  levels  due  to  the  possibility  of  separation  and  re¬ 
attachment  of  the  flow.  The  Haddle  and  Skudrzyk  (1969)  data  and  Perraud  (1989)  data 
are  representative  of  the  case  of  separation-induced  transition  on  the  nose  of  a  body  of 
revolution.  These  are  presented  in  section  2.7.1. 

The  rms  pressure  level  measured  under  the  transition  region  on  a  1.5  caliber  ogive 
nose  was  reported  by  Arakeri  (1975)  to  be  some  2  to  2.5  times  higher  than  the  rms 
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pressure  measured  further  back  under  the  fully  developed  TBL.  The  hemispherical  nose 
shape  does  promote  laminar  separation  and  the  rms  pressure  levels  measured  in  and 
around  the  separation  reattachment  point  are  10  times  higher  than  the  TBL  rms  pressure 
levels.  The  separation-induced  wall  pressure  fluctuations  are  4  to  5  times  higher  than 
those  induced  by  transition.  Similar  measurements  were  performed  by  Huang  and  Shen 
(1989).  Several  axisymmetric  headforms  were  tested  in  the  Anechoic  Flow  Facility  at 
David  Taylor  Model  Basin  where  wall  pressure  fluctuations  induced  by  either  natural 
transition  or  laminar  separation  were  determined.  They  arrived  at  the  same  conclusions 
as  Arakeri  regarding  the  magnitude  of  the  pressure  fluctuations. 

Arakeri,  et  al.  (1991)  found  a  scaling  equation  for  the  flow  noise  received  at  the 
stagnation  point  of  axisymmetric  bodies  of  the  form: 

S„(f)  -  S.W  +  50  log(-^E.)  *  20  logA)  +  A Sd  *  A Sr  (2.26) 

U  m 

where  m,p  refer  to  model  and  prototype  conditions, 

D  is  body  diameter, 

ASd  refers  to  losses  due  to  surface  diffraction,  and 
AST  refers  to  a  reduction  in  levels  due  to  transducer  configuration  and 
directivity. 

However,  the  measured  data,  shown  here  as  Figure  2.4,  appear  to  yield  a  Ujs 
dependence.  This  agrees  with  the  theoretical  prediction  of  Lauchle  (1989).  This  also 
results  in  better  agreement  with  the  Haddle  and  Skudrzyk  data  presented  as  Figure  1 1  in 
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Arakeri  et  al.  (1991),  reproduced  with  additions  as  Figure  2.5.  Therefore,  it  is  suggested 
that  equation  (2.26)  be  modified  as  (see  also  Lauchle,  1992): 

s„(f)  -  sm(f)  *  75  log(^)  *  20log(i)  ♦  AS,  *  ASr  .  (2.27) 

The  models  described  in  this  dissertation  all  suggest  that  higher  radiated  pressure 
levels  will  result  from  (1)  a  shorter  streamwise  extent  of  the  transition  zone,  (2)  an 
increase  in  the  normal  velocity  due  to  an  increase  in  Ad* /At,  or  (3)  a  decrease  in  tt.  Each 
of  these  situations  can  be  associated  with  the  noses  of  axisymmetric  bodies  and  adverse 
pressure  gradients.  An  example  of  the  effect  of  a  change  in  Ax  is  shown  in  Figure  2.6.  It 
compares  predicted  radiated  noise  using  Equation  (2.13)  for  Ax=0.06l  m  as  appropriate 
for  the  hemispherical  nose,  and  z!jc=0.15  m  as  would  be  calculated  by  Equation  (1.1)  for  a 
flat  plate  at  the  same  Reynolds  number.  A  shorter  transition  zone  results  in  a  slight 
increase  in  the  peak  amplitude  of  the  power  spectrum  and  shifts  the  peak  to  a  higher 
frequency. 

Howe  (1999a)  has  also  extended  his  model  to  axisymmetric  bodies  with  similar 
assumptions  as  described  in  section  2.6.1. 

2.6.3  Airfoils  and  Hydrofoils 

Just  as  trailing  edge  flows  are  critical  to  the  direct  radiation,  excitation,  re¬ 
radiation  and  scattering  of  noise  for  airfoil  trailing  edges  (  see  Brooks  and  Hodgson, 

1981,  Gershfeld,  1997,  Howe,  1998b,  1999b, c),  so  may  transition  near  the  leading  edge  of 
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an  airfoil  act  as  a  noise  source.  Of  course,  the  characteristics  of  the  transition  location 


and  extent  is  quite  dependent  upon  the  chordwise  distribution  of  thickness  for  the  airfoil. 
This  is  well  discussed  in  the  papers  by  Halstead  et  al.  (1995a,b,c,d)  and  Lou  and 
Hourmouziadis  (1999).  These  investigations  encompassed  operation  with  periodic  flows 
due  to  upstream  blade  wakes  in  the  former  and  a  downstream  rotating  flapper  valve  in  the 
latter  reference.  Transition  as  well  as  laminar  separation  and  turbulent  reattachment  were 
observed.  However,  in  that  case,  the  growth  of  the  instability  did  not  match  that  for 
natural  transition  on  a  flat  plate.  It  was  suggested  by  Wilde  and  Rose  (1997)  that 
transition  noise  may  be  most  important  in  cases  where  the  surface  area  covered  by 
transitioning  flow  is  comparable  to  fully  turbulent  flow  area.  This  is  potentially  the  case 
with  airfoils  used  in  turbines  which  encounter  bypass  and  separation-induced  transition 
due  to  inflow  wakes  (Addison  and  Hodson,  1991  and  Halstead  et  al.,  1995a,b,c,d).  The 
applicability  of  any  of  the  reviewed  models  for  transitional  flow  on  blades  would  most 
likely  be  the  substitution  of  the  appropriate  leading  edge  Green’s  function  into  the 
fluctuating  displacement  thickness  model  following  the  work  of  Howe  (1998, 1999a). 

2.7  Need  for  Refining  the  Fluctuating  Boundary  Layer  Displacement  Thickness 

Model  in  a  Natural  Transition  Region 

It  is  now  possible  to  make  some  refinements  to  the  Lauchle  (1981)  model  with  the 
availability  of  the  empirically  derived  Josserand  space-time  correlation  function 
(Josserand  and  Lauchle,  1990).  It  should  be  noted  that  a  similar  treatment  can  also  be 
incorporated  into  the  two-fluids  model. 
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Audet  et  al.  (1989b)  point  out  that  a  significant  difference  between  their 
prediction  from  Equation  (2.21)  and  Lauchle’s  (2.13)  is  that  they  leave  the  factor  1/  lr-r0  f2 
inside  the  integral.  This  is  important  when  calculating  the  directly  radiating  near-field 
contribution  to  measured  wall  pressure  fluctuations  whereas  Lauchle  made  a  far-field 
(and  compact  source)  approximation  of  1/  lr-r0 12  -  1/r2  which  becomes  a  constant. 
Therefore,  this  term  should  be  brought  back  inside  the  integral  for  making  predictions 
from  a  transition  zone  which  may  not  appear  compact  or  necessarily  in  the  far-field  to  a 
transducer.  This  is  often  the  case  for  buoyant  vehicle  and  sonar  self-noise  tests. 

Using  the  chosen  intermittency  space-time  correlation  function  in  Lauchle’s 
Equation  (2.12)  provides  a  refined  equation  for  the  radiated  pressure  power  spectrum: 


G(rj)  = 


p2  (A8*)2  o>4 

°°  Ax  W/2  '-n.r  1  Strand 

r  r  r  r  r  1  1 

4  7I2  [1  +(cot(.)2] 

J  J  J  J  J  Ix-yl  |jc-(y+^)| 

-»  0  W/2  -ii,  -1.5n,tana  1  1  1  v  ’  1 

(2.28) 

•  WtA-O  c‘"  ^  K  *1*  *1,  *  • 

The  limits  of  integration  for  d£t  have  been  decreased  from  infinity  to  a  correlation  length 
limit  approximated  by  the  spot  transverse  extent.  Chen  and  Thyson  (197 1)  found  this  to 
be  ±  1.5  tjx  tan  a  where  a  is  the  spot  boundary  angle  to  the  flow  axis  (a  =  10  degrees). 
Also,  assuming  homogeneity  in  the  rjz  direction,  Equation  (2.25)  becomes  a  quadruple 
integral  with  a  transition  zone  width  factor,  W,  in  the  numerator  of  the  leading  coefficient. 

The  approach  taken  by  Lauchle  using  a  Liepmann  analogy  has  the  advantage  of 
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being  able  to  incorporate  both  the  gross  piston  like  fluctuation  of  the  boundary  layer 
displacement  thickness  and  the  acoustic  sources  associated  with  the  growth  of  the  spot 
domain  involving  the  horseshoe  vortex.  The  gross  piston  fluctuation  is  described  by  AS’, 
Ax,  uc,  tj  and  N.  The  spot  size  limitation  is  accounted  for  by  1.5  rjx  tan  a.  Whatever  spot 
dynamics  (and  resulting  unsteady  pressure  sources)  haven’t  been  accounted  for  properly 
in  an  empirically  derived  space-time  correlation  at  the  wall  will  become  the  necessary 
future  refinements. 

In  trying  to  identify  these,  some  useful  information  has  been  found  on  spot 
statistics  in  Wygnanski  et  al.  (1976).  Experimental  evidence  was  cited  to  establish  that 
two  artificially  generated  spots  displaced  laterally  did  not  affect  each  others’  growth  rate. 
They  also  provide  detailed  descriptions  of  spot  dynamics  as  shown  in  Figures  2.7  and  2.8. 
The  spot  shape  as  a  function  of  normalized  surface  point  passage  time  in  Figure  2.7 
shows  the  variation  of  the  extent  of  the  spot  with  height  off  the  surface.  This  variation  is 
important  when  trying  to  model  the  spatial  distribution  of  normal  velocity  to  define  the 
motion  of  an  equivalent  "piston".  The  variation  of  vn  with  normalized  time  at  two 
elevations  (Figure  2.8)  shows  how  difficult  it  is  to  define  The  celerity  of  the 

spot's  leading  and  trailing  edges  vary  with  y.  Gutmark  and  Blackwelder  (1987)  reported 
measuring  leading  and  trailing  edge  celerities  which  increased  with  y  from  0.88 1/„  and 
0.5817.  at  y^O.lS^  to  0.93 £7.  and  0.6077.  at  y=Slam.  Similar  results  were  found  from  data 
of  other  investigators. 


40 


2.7.1  Example  Applications  of  Fluctuating  Displacement  Thickness  Model 

Measurements  of  wall  pressures  and  an  apparent  upstream  acoustic  propagation  of 
transition  zone  emissions  have  been  made  on  an  axisymmetric  0.32  m  diameter  body  in 
the  CEPRA  19  anechoic  wind  tunnel  by  Perraud  (1989).  Comparisons  are  made  to 
pressure  spectra  measured  using  large  flush-mounted  transducers.  In  this  experiment, 
r  =  0.3  m,  Ax  =  0.20  m,  u/U„  =  0.8,  and  AS'  was  estimated  to  be  0.25  mm.  The  spectral 
shape  and  amplitude  were  predicted  reasonably  well  with  0.05  ^  uj/Ax  <,  0.15. 

Computational  results  are  compared  to  the  measured  data  in  Figure  2.9.  Low- 
frequency  energy  appears  centered  around  the  most  unstable  Tollmein-Schlichting  wave 
frequencies;  this  contribution  is  not  predicted  by  the  current  theoretical  model. 

A  comparison  is  also  made  of  predicted  radiated  noise  to  the  buoyant  body  data  of 
Haddle  and  Skudrzyk  (1969).  While  there  are  also  contributions  from  the  fully  developed 
TBL  and  noise  from  the  stabilizing  fins,  which  were  found  (Lauchle  and  Brungart  1988), 
using  TBL  noise  estimates  from  Blake  (1986)  and  Lauchle  (1977),  and  fin  trailing  edge 
noise  estimates  from  Brooks  and  Hodgson  (1981),  to  be  much  lower  than  the  measured 
noise.  Results  for  a  typical  operational  condition  of  U„  -  15.44  m/s,  Ax  =  0.061  m,  r  - 
1.0  m,  u/U„  =  0.8,  body  diameter  of  0.48  m  and  AS'  estimated  as  0.43  mm  have  been 
made.  These  are  all  shown  in  Figure  2.10.  The  predicted  spectral  levels  for  the  directly 
radiated  noise  from  the  fully-developed  TBL  seems  inconsistently  low  relative  to  the 
predicted  transition  noise  when  compared  with  the  airflow  measurements  described  in 
Chapter  6.  It  is  believed  that  several  factors  may  contribute  to  this.  The  low  wavenumber 
data  used  to  estimate  the  TBL  noise  may  have  been  underestimated  by  the  wavevector 
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filtering  used  in  the  experiment.  The  scaling  relationship  used  by  Blake  assumes  a  rigid 
flat  wall  and  also  neglects  the  effect  of  proximity  of  the  source  to  a  nose  of  the  body  as 
examined  by  Howe  (1999a).  There  could  also  have  been  some  separation  and 
reattachment  of  flow  near  the  noise  due  to  slight  variation  in  angle  of  attack  as  the 
buoyant  body  ascended. 

Neither  application  of  the  model  is  particularly  satisfying.  Likewise,  each  of  the 
examples  presented  by  Lauchle  (1981)  seem  to  have  the  best  fit  of  predicted  versus 
measured  data  for  non-dimensional  rise  time,  0.45  <  Uj/Jx,<  1.0.  However,  this  tends 
to  describe  only  the  shape  of  the  higher  frequency  portion  of  the  spectrum.  Clearly, 
improvements  to  the  model  are  in  order. 

It  is  important  to  point  out  that  in  underwater  flow  situations,  the  common 
presence  of  inhomogeneities,  particularly  bubbles,  can  account  for  a  considerable  increase 
in  radiated  noise.  For  example,  Crighton  (1975)  shows  that  air  bubble  concentrations 
with  a  void  fraction  of  less  than  1%  in  turbulent  flow  can  increase  the  radiated  noise 
power  by  up  to  50  dB  due  to  excitation  by  passage  of  the  bubbles  through  high  pressure 
gradient  flows. 

2.7.2  Study  Objectives:  General  Approach  Towards  Characterizing  the  Acoustic 

Emissions  from  Isolated  Transition  Flow  Structures  and  a  Natural 

Transition  Zone 

Given  the  possibility  of  high  radiating  energy  from  the  transition  zone, 
understanding  the  boundary  layer  transition  process  is  perhaps  as  important  from  the 
point  of  view  as  a  structural  excitation  mechanism  as  it  is  as  a  direct  radiator.  It  is  beyond 
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the  scope  of  this  dissertation  to  pursue  all  of  the  fluid-structure  coupling  implications  and 
local  excitation  of  bubbles  which  must  be  dealt  with  to  truly  predict  the  transition  related 
signature  in  an  underwater  situation.  However,  this  may  be  the  role  in  which  transition  is 
most  important  for  problems  of  practical  concern.  It  certainly  does  not  appear  that  direct 
radiation  from  such  a  small  area  as  the  transition  region  (and  probably  only  half  of  it  in 
terms  of  spot  growth)  is  significant  when  compared  to  fully  developed  TBL  noise  over 
the  surface  area  of  the  rest  of  the  body  or  plane  in  zero  or  variable  pressure  gradients. 
Nevertheless,  we  are  interested  in  pursuing  improvements  in  both  aerodynamic  and 
hydrodynamic  modeling  of  boundary  layer  transition  noise  mechanisms  and  validation 
through  quality  experiments  because  it  is  reasoned  that  transition  flow  structures  form  the 
basic  building  blocks  of  all  subsequent  downstream  flow  behavior  and  its  acoustics. 

The  problem  then  is  that  two  basic  modeling  approaches  exist  for  the  direct 
radiation  from  a  naturally  transitioning  boundary  layer.  While  the  two-fluids  model  may 
be  the  more  inclusive,  it  is  also  the  more  difficult  to  apply  in  practice.  The 
phenomenological  Liepmann  approach  is  the  most  amenable  to  incorporation  of  recent 
experimental  descriptions  of  the  normal  velocity  fluctuations.  This  study  has 
concentrated  on  Liepmann  model  modifications  and  further  qualitative  assessments.  A 
series  of  experiments  have  been  performed  for  evaluation  of  the  models. 

At  this  point,  the  value  of  direct  numerical  simulation  may  be  to  provide  insight 
on  the  necessary  modeling  of  substructure  dynamics  that  influence  the  spatial  definition 
of  the  normal  velocity  used  in  the  fluctuating  displacement  thickness  model,  and  hence 
contribute,  in  their  own  right,  to  the  acoustic  pressure,  particularly  that  in  the  nearfield. 
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This  bridges  the  gap,  to  some  extent,  between  the  two-fluids  model  using  the  Lighthill 
analogy,  and  the  Liepmann  approach  used  for  the  fluctuating  displacement  thickness 
model. 


Descriptions  of  the  Fluid  Motions 
associated  with  an  isolated 
Turbulent  Spot  in  a  Laminar 
Boundary  Layer. 


Figure  2.1 


Vortex  view  of  a  turbulent  Spot  [from  Audet  et  al.  1989a]. 


Figure  2.2  The  normalized  rms  wall  pressure  fluctuations  in  the  transitional  boundary 
layer  measured  and  predicted  for  a  flat  plate.  The  curves  are  predictions 
which  account  for  spatial  and  frequency  filtering  of  the  sensors.  Curve  1 
is  ypT,  Curve  2  is  the  radiated  acoustic  source  term,  and  Curve  3  is  their 
sum.  The  symbols  are  measured  data  for  3  velocities,  [from  Audet  et  al. 
1989a,b]. 
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Figure  2.5  Verification  of  the  proposed  scaling  law  with  the  Haddle  and  Skudrzyk 

(1969)  results. - ,  OSc  model  results  with  U=  12.2  m/s  and  diameter 

D— 0.075  m,  o,  model  results  extrapolated  to  Haddle  and  Skudrzyk 
conditions  of  £4=18.5  m/s  and  D=0.42  m  without  corrections;  -  - 
extrapolated  results  corrected  for  diffraction  loss  scaled  on  two  different 
power  laws;  A,  the  Haddle  and  Skudrzyk  experimental  results  [the  5th 
power  law  is  from  Arakeri  et  al.  1991;  the  7.5  power  law  is  from  Lauchle 
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Figure  2.6  Comparison  of  radiated  1/3  octave  levels  predicted  by  Equation  (2.13) 
for  the  Haddle  and  Skudrzyk  (1969)  test  conditions  using  a  flat  plate 
(0.15  m)  and  a  hemispherical  noise  prediction  (0.061  m)  for  transition 
zone  length  Ax. 
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Figure  2.7  Two  cross-sections  through  a  spot:  a  comparison  of  criteria,  U  =9.4  m/s 
(a)  elevation,  (b)  plan,  -x-.  Coles  and  Barker  (1974);  — ,  Wygnanski  et  al. 
(1976).  Shape  based  on  2%  elevation  in  U:  ••••,  LE  ensemble  average; 
— ,  TE  ensemble  averaged  [from  Wygnanski  et  al.  19761. 


Figure  2.9  Pressure  spectra!  densities  measured  under  the  laminar  boundary  layer  on 
an  axisymmetric  body  in  the  CEPRA  19  wind  tunnel  and  comparison  with 
computational  model  results  [from  Perraud  1989], 
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Figure  2. 10  Measured  and  predicted  radiated  noise  spectra  for  the  Haddle  and 

Skudrzyk  (1969)  buoyant  bodies  operating  at  15.44  m/s  [from  Lauchle 
and  Brungart,  1988]. 


CHAPTER  3 


MODEL  REFINEMENTS 


A  major  objective  of  this  research  program  is  to  continue  to  develop  refinements 
to  the  radiated  noise  models  for  transition  noise.  These  theoretical  extensions  proposed 
are  necessary  to  complete  the  area  of  investigation  of  transition  zone  acoustics  in  a  zero 
pressure  gradient.  As  discussed  previously,  there  are  two  major  approaches  which  can  be 
summarized  in  Figure  3.1.  Lauchle  began  with  a  Lighthill  acoustic  analogy  and  moved  to 
a  Liepmann  analogy  while  Somette  and  Lagier  pursued  a  more  rigorous  Lighthill 
approach.  These  are  indicated  by  the  top  and  bottom  boxes  in  the  figure.  The  right  hand 
box  represents  the  use  of  the  theoretical  space-time  correlation  function  whose  form  was 
estimated  using  phenomenological  assumptions.  Lauchle  (1981),  Lagier  and  Somette 
(1986),  and  Audet,  et  al.  (1989)  used  this  form  to  make  estimates  of  radiated  noise 
spectra.  Similarly,  the  left  hand  box  represents  the  semi-empirical  form  developed  by 
Josserand  and  Lauchle  (1990).  The  dashed  line  between  the  two  correlation  functions 
indicates  an  attempt  to  derive  the  constants  of  the  theoretical  function  from  the  empirical 
function.  The  lines  and  dates  linking  the  models  to  the  correlation  function  form  boxes 
indicating  when  certain  modifications  were  made  to  the  models  as  they  were  carried 
toward  validation  with  measured  spectra.  The  following  sections  describe  some  of  the 
further  refinements  and  analysis  that  are  being  performed.  In  order  to  satisfy  the  input 
requirements  for  both  models,  it  has  been  necessary  to  review  previous  fluid  dynamic 
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experiments  which  have  identified  structures  within  turbulent  spots  and  developed 
empirical  descriptions  for  spot  growth. 


3.1  Fluctuating  Boundary  Layer  Displacement  Thickness  Using  a  Liepmann 
Analogy  •  Theoretical  Space-Time  Correlation  Function 

For  turbulent  spot  substructures,  one  can  describe  the  resultant  normal  velocity 

distribution  with  two  terms,  using  the  substantial  derivative: 


D5* 

Dt 

...  TI  56*  x  58* 

v„(jc,0  =  £/ -  +  - 

"  dx  dt 


(3.1) 


The  first  component  describes  the  growth  of  a  spot  in  a  fixed  frame  of  reference.  It 
should  have  strong  statistical  descriptors  including  spatial  correlations.  It  has  been  found 
by  Krane  (1992)  that  the  normal  velocity,  described  by  the  first  component,  is  directed 
toward  the  wall  at  the  spot  leading  interface  and  outward  over  the  rest  of  the  spot  which 
requires  a  higher  normal  velocity  for  the  potential  flow.  A  streamwise  dependence  was 
measured  by  Krane  (1992).  A  spanwise  dependence  has  not  been  measured.  The 
adjacent  inward  and  outward  velocities  produce  a  weak  dipole  at  low  Mach  number.  The 
second  term  is  due  to  the  growth  of  the  spot  from  the  creation  of  internal  vorticity.  An 
instantaneous  vorticity  field  was  measured  by  Sankaran  et  al.  (1988, 1991).  This  is 
shown  in  Figures  3.2  and  3.3.  The  vortex  axes  are  observed  to  be  inclined  to  the  plane 
wall.  As  sources  of  sound,  these  vortices  are  likely  to  be  quadrupole  in  nature.  It  has  been 
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argued  by  Somette  and  Lagier  (1984a)  that  weak  monopoles  are  created  due  to:  (1) 
entropy  production  via  viscous  dissipation,  and  (2)  kinetic  energy  dissipation.  Spatial 
growth  may  occur  due  to  the  streamwise  increase  in  the  number  of  vortical  substructures. 
However,  when  Equation  (3.1)  is  handled  with  retarded  time  variables,  spatial  derivatives 
become  time  derivatives  according  to: 


d_  ^  _a 

dx  c  dt 


Then,  it  reduces  to: 


^  ,  U„  as*  58* 

v"fc')  =  '  T  17  *  IT 


=  (1  -M.)  — 

dt 


38^ 

dt 


for  «  1 


(3.2) 


(3.3) 


So  we  really  want  to  concentrate  on  those  sources  which  describe  the  temporal 
fluctuations  of  the  displacement  thickness,  which  are  obviously  due  to  the  creation  and 
interactions  of  substructure  vortices. 

3.1.1  Computational  Model  Refinements 

Several  refinements  have  been  performed  for  the  computer  program  GLPROG81 
which  implements  Equation  (2.13).  These  are  described  in  the  rest  of  this  section. 
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Similar  refinements  were  also  made,  where  appropriate,  to  GLPROG80  which 
implements  Equation  (2.6)  for  the  fluctuating  shear  stress  model..  Examples  of  all  of 
these  Mathcad  files  are  provided  in  Appendix  A. 

3.1.1. 1  Conversion  of  Software  Environment 

The  entire  program  was  converted  from  FORTRAN  with  its  many  specific 
software  package  mathematical  and  graphical  function  calls  to  the  Mathcad  2000 
Professional  ®  (Mathsoft,  Inc.)  software  package.  This  has  several  advantages.  First,  it 
is  a  Microsoft  Windows  supported  personal  computer  product  rather  than  written  for  use 
on  a  Digital  Vaxstation  computer.  Second,  the  math  is  presented  in  true  mathematical 
symbolic  form  rather  than  a  FORTRAN  breakdown  into  variable  names  and  code  lines. 
This  makes  the  presentation  much  easier  to  follow  and  understand. 

3.1.1.2  Intermittency  Function  Models 

Josserand  derived  a  new  intermittency  factor  from  his  experimental  data: 

y(ti  )  =  l-e  'd+3.4Z)z2  where  Z=—  (3.4) 

1  U  Ay 


instead  of 


yOljc)  =  i-« 


-4.185Z2 


(3.5) 


which  was  developed  by  Farabee,  et  al.  (1974).  The  first  function  follows  from  a 
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formulation  developed  by  Gedney  (1979)  (and  Gedney  and  Leehey,  1989)  for  a  combined 
Dirac  line  source  and  constant  source  hypothesis.  The  comparison  of  these  two  relations 
is  shown  in  Figure  3.4.  The  Farabee  relation  is  seen  to  have  its  highest  slope  earlier  in  the 
transition  zone.  This  carries  through  to  modest  increases  in  the  predicted  amplitudes  as 
well.  The  Josserand  model  predicts  intermittencies  that  agree  with  experimental  data 
better  than  the  Farabee  model.. 


3.1. 1.3  Burst  Frequency  Models 

A  new  spot  burst  frequency  description  was  also  derived  by  Josserand  (1986)  based  on 
the  Gedney  (1979)  hypothesis  mentioned  in  Section  3. 1.1.2: 


M 


N(ZX)  =  2.38  (-2-) 
Ax  M 


[1-Y (ZJ]  In[ 
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i-y  CZJ 
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(3.6) 


instead  of 


N(Z  )  =  1.272  (— )  Z  e-4185(Zi)2  (3.7) 

Ax 

due  to  Farabee,  et  al.  (1974).  Josserand  modified  the  constant  in  Equation  (3.6)  to  yield  a 
better  match  with  his  experimental  data.  Note  that  Equation  (3.6)  can  be  evaluated  with 
the  intermittency  function  described  by  either  Equation  (3.4)  or  (3.5).  A  comparison  of 
these  three  formulations  is  shown  as  a  normalized  burst  rate,  B(Z)  (where  N(ZX)  is 
multiplied  by  Ug/Ax),  in  Figure  3.5.  It  can  be  seen  that  the  magnitudes  of  the  peak  burst 
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rate  for  the  Josserand  function  variants  are  about  4  time  larger  than  that  for  the  Farabee 
model.  The  non-zero  values  of  B(Z)  predicted  by  Equation  (3.6)  beyond  the  end  of  the 
transition  zone  seem  physically  unrealistic. 

Of  course,  our  real  interest  is  in  how  the  use  of  these  different  relations  will  affect 
the  calculation  of  the  radiated  noise.  The  results  of  these  combinations  of  substitutions  of 
intermittency  and  burst  rate  are  shown  for  one  case  in  Figure  3.6.  The  measured  data 
come  from  a  pressure  sensor  flush  mounted  in  the  laminar  flow  region  on  a  body  of 
revolution  in  the  CEPRA  19  wind  tunnel  (Perraud  1989).  The  parameters  for  that  test 
affecting  the  computational  model  inputs  for  GLPROG81  are  listed  in  the  caption.  The 
normalized  rise  time  of  0.15  was  suggested  by  Perraud.  This  new  spot  burst  frequency 
description  results  in  a  6  to  8  dB  increase  in  predicted  radiated  noise  over  much  of  the 
frequency  range  of  interest.  The  impact  of  a  change  of  the  intermittency  model  is 
negligible  compared  to  that  of  the  burst  rate.  The  Farabee  form  (Equation  3.7)  of  the 
burst  frequency  function  results  in  a  slightly  lower  peak  frequency  as  well  and  a  larger 
sound  pressure  magnitude  below  50  Hz. 

However,  because  Josserand  (1986)  reported  a  variance  of  20%  on  his  calculated 
burst  rate,  which  makes  his  empirical  fit  somewhat  questionable,  it  seems  likely  that  a 
portion  of  the  predicted  high  levels  of  sound  are  due  to  this  over-prediction  of  the  burst 
rate. 

3.1.1.4  Application  of  Strip  Analysis 

Multiplication  of  the  result  of  Equation  (2.13)  by  the  transverse  width  of  the 
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transition  zone,  W,  yields  the  total  radiated  spectral  pressure  from  the  transition  zone. 

This  assumes  limited  spanwise  spot  interaction  which  leads  to  homogeneity. 

3.1.2  Incorporation  of  Measured  Normal  Velocity  Fluctuations 

Lauchle’s  assumed  monopole  description  of  the  normal  velocity  in  GLPROG81 
has  a  provision  for  substitution  of  the  normal  velocity  resulting  from  the  water  channel 
measurements  of  Krane  (1992).  He  measured  the  normal  velocity  using  a  mass  flux 
deficit  to  describe  the  displacement  thickness  (Figure  3.7).  The  velocity  profile 
measurement  technique  was  developed  by  Bruneau  and  Pauley  (1992, 1995).  This 
technique  tends  to  average  out  the  time  dependent  fluctuations  of  38* ,  but  it  does  reveal  a 
negative  mass  flow  deficit  combined  with  the  previously  assumed  (by  Lauchle,  1981) 
positive  one.  Krane  identified  the  streamwise  dependence  with  a  weak  dipole  which 
results  from  the  negative  normal  velocity  near  the  leading  edge  (overhang)  of  the 
turbulent  spot  and  a  positive  velocity  associated  with  the  passing  of  the  rest  of  the  spot. 
The  multipole  source  distribution  proposed  by  Krane  is  shown  in  Figure  3.8. 

Krane  and  Pauley  (1993)  also  found  that  the  interaction  of  two  turbulent  spots  had 
a  minimal  effect  on  the  mass  flux  deficit  (displacement  thickness)  and ,  therefore,  the 
normal  velocity  descriptor.  This  would  imply  that  their  normal  velocity  model  should 
apply  to  even  a  natural  transition  case.  For  both  the  single  and  the  interacting  spots,  a 
dependence  on  Zx  was  measured.  A  function  was  fit  to  this  data, 

TK(Zx)  =  0.30 -exp[- 1.45  -  3(Z/]  .  (3.8) 
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This  function  is  shown  compared  to  Krane’s  measurements  in  Figure  3.9. 

The  dipole  source  shown  in  Krane  (1992,  Figure  5.1)  was  used  to  replace  the 
monopole  source  in  GLPROG81.  The  proposed  quadrupole  sources  were  not  included 
because  of  their  reduced  radiation  efficiency  compared  to  the  dipole  and  they  do  not 
appear  prominent  for  all  velocities  shown  by  Krane.  However,  the  amplitude  of  the 
normal  velocity  appears  to  be  equal  for  the  positive  and  negative  events  and  was  assumed 
to  be  the  same  as  the  positive  event  modeled  by  Lauchle.  Therefore,  the  dipole  was 
modeled  (as  shown  in  Pierce  (1981,  pp.170)  as  modifying  the  monopole  pressure  by  the 
factor  of  2  sin[(k  d  cos0)/2]  where  d  is  the  separation  distance  of  the  two  sources  that 
make  up  the  dipole  and  the  observation  distance,  r»d.  This  can  be  further  simplified 
for  kd«l,  to  k  d  cos0.  This  would  be  appropriate  for  the  low  frequencies  for  which 
transition  noise  seems  important.  The  limitation  in  using  this  in  the  model  is  that  it 
would  result  in  and  amplitude  independent  of  frequency  at  higher  frequencies. 

The  separation  distance  can  be  modeled  as  2  uc  tj  which  is  then  equivalent  to 
d(Z)  =  2  A*  TKiZ)  ~  d  =  0.4  Ax.  So  GLPROG81  could  be  modified  by  multiplying 
the  previous  result  by  4  sin2[(0.2  k  Ax  cos0)]  .  However,  if  the  separation  distance 
between  the  two  dipole  sources  is  taken  to  be  a  random  variable,  then  the  mean  value  of 
the  pressure  spectrum  is  given  by 

<m  -  M  --  Md) ]  (3.9) 

where  Gmonopole  is  given  by  Equation  2.13  with  the  assumption  of  spanwise  homogeneity. 
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G(r,  ©)  = 


Wp2  (AS*)2  u0  u2c 


(3.10) 


8  712  r2 


F*(kc Ax,  a* A*, 


V/ 

Ax’ 


The  expected  value  is  calculated  for  the  following  dipole  source  multiplier 

fid')  =  4  sin2(-U  d'  cosG)  (3.11) 

This  expected  value  is  determined  with  the  assumption  of  a  uniform  probability 
distribution  of  d’  and  limits  of  0.1  Ax  <  d’<0.2Ax. 

-  /  p(d')Ad)  dd'  =  jj-L-  f  Ad1)  dd'  (3.12) 


This  then  becomes 


Elfid)]  =  2+[—  sm(0A.k  A*  cos9)]  [1-2  cos(0.1  k  Ax  cosG)] 
0.1  it  Ax  cos0 


(3.13) 


3.2  Fluctuating  Boundary  Layer  Displacement  Thickness  Using  a  Liepmann 
Analogy  -  Empirical  Space-Time  Correlation  Function 

3.2.1  Integration  of  Josserand-Lauchle  Correlation  Function 

As  discussed  previously,  the  substitution  of  the  semi-empirical  space-time 
correlation  function  developed  by  Josserand  (1986)  into  Equation  (2.12)  was  coded  for 
integration.  The  function,  from  Josserand  and  Lauchle  (1990),  is  repeated  here: 
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=  y-  ^  +  (1^>  exp(-  5  (4  +  200  |T'^)  |T'^  } '  exp(  o.ou  T^yJ  |t"^ 

1260  |t-^|  |  a, |  J  A  |  A  |  1 

- —  - - -  •  exp - 

1  +  71  | A,| |  1  +  14.2  |AjJ  0.0014  ♦  |AJ  ,  +  130o|t-i|* 

Ue 

where  A  =  -ln[l-exp( -4.27/Ax)],  (units  of  Ax  are  in  meters), 

A,  =  \  /  Ax  , 

yu  =  y(t]J  and  yd  -  y(tjx+£x)  when  the  reference  location  is  upstream  of  the  second 
location, 

yu  =  y(tjx+£J  and  yd  =  y(rjx)  when  the  reference  location  is  downstream  of  the 
second  location. 

Integration  of  this  function  has  proven  difficult  for  the  Mathcad  package.  The 
built-in  integration  was  bypassed  with  a  trapezoidal  integration.  Tthe  use  of  a  high 
number  of  integration  points  smooths  this  to  a  large  degree.  Several  approaches  have 
been  tried  to  ease  computational  burden.  One  is  to  set  physically  reasonable  limits  of 
integration.  This  resulted  in  limits  for  £  from  -rjx  to  1  -rjx  and  £  from  -1.5  rjx  tan  a  to  1.5 
tjx  tan  a.  In  an  attempt  to  get  better  use  of  the  computational  dynamic  range,  the  mean  was 
removed  from  the  cross-correlation  as  indicated  in  equation  (3.15)  and  the  last  term  of  the 
following  equation  was  ignored. 

■y®)  =  fR*y(T)e  ^dx  =  +  MA5(w)  (3.15) 
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Because  of  the  difficulty  in  integrating  the  Josserand  (1986)  form  of  the  space- 
time  correlation  function,  an  effort  was  made  to  look  for  simplification  of  the  function  to 
make  integration  easier.  However,  it  does  not  appear  that  this  function  has  a  form  that  is 
direction  separable  for  x  and  z. 

3.2.2  Chase  Statistical  Modification 

Chase  (1989)  has  suggested  that  the  first  exponent  term  in  Equation  (3.14)  can  be 
replaced  by  a  more  general  term  based  on  the  probability  of  state  transition  to  turbulence 
of  form: 


exp< 


yd  (  l^d) 


(3.16) 


However,  Josserand  found  that  both  terms  predict  the  measured  correlation  values  well 
for  0.3<^0.7;  but,  that  outside  of  these  limits,  the  Chase  theoretical  term  tended  to  over¬ 
predict  the  correlation.  Both  terms  are  available  as  options  in  the  computational  software. 
Limit  evaluation  has  determined  that  the  Chase  modification  to  the  Josserand-Lauchle 
correlation  function  does  not  improve  the  results  in  comparison  to  measured  data  for 
natural  transition.  A  comparison  of  the  Josserand  and  the  Chase  modified  correlation 
surfaces  over  the  transition  zone  is  shown  in  Figure  3.10.  The  differences  are  subtle  with 
the  Chase  term  having  deeper  troughs  in  the  amplitude  surface. 
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3.3  Two  Fluids  Model  with  Refined  Normal  Velocity 

In  tying  together  all  the  approaches,  the  Audet  and  Somette-Lagier  approaches  did 
not  provide  any  advances  because  of  simplifying  assumptions  to  allow  integration  and 
heavy  dependence  on  Lauchle’s  (1981)  functions.  The  refined  normal  velocity  description 
could  be  substituted  into  the  two  fluids  model  but  it  would  not  provide  an  advantage  over 
the  refined  fluctuating  displacement  thickness  model.  Knowledge  gained  of  the  radiated 
noise  component  magnitude  can  be  used  to  further  refine  the  two-fluids  model  of  Audet, 
et  al.  (1989b)  for  wall  pressure  fluctuations  as  input  to  the  fluid-structure  interaction 
problems  of  interest. 

Their  predictions  used  simplifying  assumptions  for  their  more  general  equation. 
The  retarded  time  in  the  equation  does  not  take  into  account  the  separation  distance  in  the 
2  point  correlation.  It  accounts  only  for  the  reference  point.  Global  stationarity  is 
assumed.  No  spot  interaction  is  assumed.  All  values  are  taken  for  tj= 0.  Spatial  ' 
dependence  is  separated  from  temporal  dependence. 

The  good  comparison  of  the  Lauchle  (1981)  model  with  CEPRA  wind  tunnel 
data  when  using  normalized  rise  times  of  0.05  to  0.15  must  be  viewed  with  reserve  since 
the  authors  indicated  that  some  data  were  contaminated  by  comer  vortex  growth 
overlapping  the  sensor  locations. 

3.4  Accounting  for  Vortical  Motions 

Another  objective  of  this  study  is  to  investigate  the  feasibility  of  incorporation  of 
effects  due  to  vorticity  into  either  of  Lauchle’s  models.  This  will  mainly  draw  upon  the 
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prior  work  of  Powell  (1995),  Hardin  (1991),  Lund  (1989),  and  the  many  papers  of  Kambe 
and  Obermeier.  Williams,  Fasel,  and  Hama  (1984)  also  made  a  great  experimental 
contribution  with  their  mapping  of  the  three-dimensional  vorticity  field  in  boundary  layer 
transition.  It  is  anticipated  that  this  would  form  a  second  source  term  in  the  fluctuating 
displacement  thickness  model  in  a  manner  similar  to  the  S  and  S0  functions  of  Lagier  and 
Somette  (1986).  Lagier  and  Somette  concluded  that  the  contribution  of  a  shear  interface 
at  the  edge  of  a  turbulent  spot  (due  in  large  measure  to  the  presence  of  these  vortices)  was 
insignificant  compared  to  the  overall  growth  of  the  spot. 

Sankaran,  Sokolov,  and  Antonia.  (1988)  and  Sankaran,  Antonia,  Bisset,  and 
Sokolov  (1991)  used  instantaneous  velocity  field  measurements  to  identify  a 
superposition  of  vortex  sources  within  a  turbulent  spot.  The  existence  of  vorticity 
formations  may  help  explain  the  Uj*  noise  dependence  by  looking  at  the  sources  as 
primarily  quadrupoles  with  some  monopole  like  bursting.  The  measurements  indicate 
that  the  spatial  growth  of  the  spot  is  directly  related  to  a  streamwise  increase  in  the 
number  of  substructures.  Substructures  are  generally  created  near  the  trailing  edge  of  a 
spot,  and  all  convect  at  the  same  velocity  0.8  Um.  This  means  that  the  cumulative 
contribution  from  these  substructures  will  be  greatest  near  the  midpoint  of  the  transition 
zone  because  the  number  of  spots  is  greatest  there.  Would  incorporation  of  a  distribution 
of  these  vortices  help  account  for  the  highly  spatially  dependent  normal  velocity 
component  due  to  587 dr?  Many  would  argue  that  the  hydrodynamic  fluctuations  due  to 
these  substructures  is  statistically  the  same  as  what  would  be  found  with  a  fully  turbulent 
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boundary  layer  however,  the  spatial  and  temporal  extent  is  limited  by  the  axial 
dependence  of  intermittency  and  the  burst  rate.  Hence,  the  vortical  contribution  is 
captured  in  the  use  of  TBL  statistics  in  models  such  as  Howe’s.  Therefore,  no  additional 
changes  are  made  to  any  of  the  models  since  they  either  utilize  the  modified  TBL 
statistics  or  not. 

3.5  Incorporation  of  Green’s  Function  for  Transition  Zone  Near  Half-Plane 

Leading  Edge 

Howe’s  model  as  described  in  section  2.6.1  has  been  implemented  in  a  Mathcad 
worksheet,  MSHTRANS99,  that  is  provided  in  Appendix  A.4. 

The  trick  on  the  integration  of  Equation  (2.25)  is  to  apply  a  window  function  that 
varies  with  frequency.  A  longer  window  is  required  for  low  frequencies  and  a  short 
window  is  required  for  high  frequencies.  Also,  to  avoid  aliasing,  a  much  higher  sampling 
rate  is  necessary  for  high  frequencies  than  for  the  low  frequencies.  Trapezoidal 
integration  is  also  used  rather  than  the  Mathcad  built-in  function. 

Howe’s  original  formulation  uses  the  intermittency  function  shown  in  Equation 
(3.5).  An  alternative  Mathcad  worksheet  uses  the  Josserand  intermittency  function 
(Equation  3.4).  The  impact  is  about  a  4  dB  decrease  in  the  radiated  pressure  levels  with 
the  Josserand  intermittency  function. 

An  alternative  use  of  the  formulation  provided  by  Howe  would  be  to  substitute  the 
pressure  spectrum  for  the  wall  pressures  as  determined  by  Josserand  (1986)  rather  than 
intermittency  modified  Chase  (1980)  description  of  TBL  wall  pressures  for  the  section  of 
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the  finite  thickness  half-plane  (hydrofoil)  surface  S  covering  the  transition  zone.  This  has 
not  yet  been  pursued. 
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1 


yFsnsb^Z)  :=  1  -  exp(-4. 185Z2)  Empirical  rolaftcns  from  Farabs®,  @t  al.  (1974). 

yJL(Z)  :=  I  -  exp[-(l  +  3.4Z)-Z2]  Empirical  r©la$on©  from  Josserand  and  lauchl© 

(1 980)  which  was  based  on  th®  work  of  G@dn©y 
(1979)  asessing  both  a  Dime  tin®  source  hypothesis 
and  a  constant  source  hypothesis  for  int©rmft©ncy. 


gure  3.4  Comparison  of  empirically  derived  intermittency  functions. 
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Figure  3.6 


Effect  of  substitution  of  refined  empirically  derived  equations  for  the 
intermittency  function,  y(qx) ,  and  turbulent  spot  burst  frequency,  N(r\x), 
for  CEPRA  19  axisymmetric  body  in  air,  Uj=  50  m/s,  uct;/dx=0.15, 

Ad*- 0.25  mm  and  Ax=0.2  m  (data  from  Perraud  1989).. 
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Figure  3.9  Comparison  of  Krane  (1992)  measured  non-dimensional  rise  time,  \\QtJAx, 
and  function  TK(ZX)  fit  to  the  data. 


Variation  with  x  and  with  i;x=1.0. 


Variation  with  i;x  and  with  x=0. 


Figure  3.10  Comparison  of  correlation  surface,  Rj(r|x,^x,0,t),  for  empirical 

space-time  correlation  function  for  original  Josserand  form  and 
Chasefirst  exponential  term  statistical  modification  for 
U=11.8  m/s,  Ax=1.83m,  with  r|x/Ax=0.50. 
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CHAPTER 4 


RADIATED  SOUND  POWER  PREDICTIONS 


4.1  Model  Comparisons 

The  three  transition  noise  models  based  on  Lauchle’s  theoretical  space-time 
correlation  function  were  each  exercised  to  predict  far  field  radiated  noise  from  natural 
transition  a  flat  plate  with  zero  pressure  gradient  with  seven  velocities  for  air  flows  and 
seven  velocities  for  water  flows  in  appropriate  ranges.  The  model  which  required  the 
integration  of  the  Josserand  empirical  space-time  correlation  function  was  run  for  selected 
conditions  for  comparison  to  the  other  models.  The  higher  frequency  results  for  that 
model  are  not  considered  to  be  appropriate  yet.  The  Mathcad  worksheets  for  these 
models  are  all  provided  in  Appendix  A: 

GLPROG80.MCD  (implementing  Equation  2.6,  Appendix  A.l), 

(G.  Lauchle  PROGram  1980) 

GLPROG8 1  .MCD  (implementing  Equation  2.13,  Appendix  A.2), 

(G.  Lauchle  PROGram  1981) 

KMTRANS00.MCD  (implementing  section  3.1.2,  Appendix  A.3), 

(Krane,  Marboe  TRAnsition  NoiSe  2000) 

JLTRANS99.MCD  (implementing  Equations  2.12  with  2.15,  Appendix  A.4),  and 
(Josserand,  Lauchle  TRAnsition  NoiSe  1999) 
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MSHTRANS99.MCD  (implementing  Equation  2.24,  Appendix  A.5) 

(M.  S.  Howe  TRAnsition  NoiSe  1999). 

The  models  have  been  run  using  either  or  both  of  the  Farabee  (Equations  3.5  and  3.7)  or 
the  Josserand  (Equations  3.4  and  3.6)  intermittency  and  spot  burst  frequency  functions. 

The  common  constants  for  all  runs  in  air  and  water  respectively  are  provided  in 
Table  4.1.  The  necessary  model  inputs  which  are  freestream  velocity  dependent  were 
calculated  using  TRANSITION  AIR.MCD  (Appendix  A.6)  and  TRANSITION 
WATER.MCD  (Appendix  A.7).  These  are  listed  in  Table  4.2. 


Table  4. 1  Common  Constants  for  Model  Exercising 


Variable 

Air 

Dimension 

Fresh  Water 

Uc/Uo 

0.7 

0.7 

W 

0.305 

m 

0.305 

r 

1 

m 

1 

Q 

1.2 

kg/m3 

999 

CO 

343 

m/s 

1525 

l) 

1.51  x  lO'5 

m2/s 

1.01  x  lO"6 

0 

0 

degree 

0 

Ret 

1500000 

1500000 

Pref 

20  x  10"6 

Pa 

1  x  lO"6 
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Table  4.2  Calculated  Model  Inputs 


Air 

FreshWater 

Uo 

Ax 

Ad* 

Uo 

A6* 

m/s 

cm 

cm 

m/s 

WM5EMM 

cm 

9 

86.57 

0.649 

3 

0.148 

12 

64.93 

0.487 

6 

lllK*YnG!l 

0.074 

15 

51.94 

0.389 

9 

0.049 

18 

43.28 

0.325 

12 

4.93 

0.037 

37.10 

15 

3.95 

18 

3.29 

27 

28.86 

0.216 

21 

2.82 

0.021 

4.1.1  Comparisons  of  Models  for  Air  Flows 

Comparisons  of  radiated  noise  from  the  shear  stress  dipole,  Liepmann  analogy 
monopole,  and  dipole  models  over  a  range  of  freestream  velocities  from  9  to  27  m/sec  are 
provided  in  Figures  4.1, 4.2,  and  4.3  respectively.  The  model  with  the  Josserand 
empirical  correlation  function  was  only  run  for  24  m/sec  (Figure  4.4).  The  Howe  model 
for  enhancement  of  a  natural  transition  near  the  leading  edge  of  a  half-plane  was  also  run 
for  the  24  m/sec  case  as  a  comparative  case  to  demonstrate  the  effect  in  the  change  to  the 
Green’s  function  (Figure  4.5).  All  of  the  models  are  compared  in  Figure  4.6  for  the  24 
m/sec  case. 

It  is  apparent  that  the  shear  stress  dipole  predictions  are  much  too  low  by  30  dB. 
The  Josserand  functions  tend  to  raise  the  levels  by  12  dB  compared  to  the  Farabee 
functions.  The  peak  frequency  in  the  spectra  (around  100  Hz)  increases  with  velocity  as 
expected.  The  Krane  dipole  results  in  lower  levels  for  frequencies  below  100  Hz  but  3 
dB  higher  levels  than  the  monopole  model  at  higher  frequencies.  This  is  because  of  the 
assumption  that  the  source  amplitude  (proportional  to  normal  velocity)  of  each  dipole 
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source  is  the  same  as  for  the  monopole,  compactness,  and  an  assumed  uniform 
distribution  of  dipole  separation  distances  [d  =  2  uct  =  2  Ax  (uct/Ax)].  This  tends  to 
smooth  out  the  nulls  associated  with  the  sine  function. 

4.1.2  Comparisons  of  Models  for  Water  Flows 

Comparisons  of  radiated  noise  from  the  shear  stress  dipole,  Liepmann  analogy 
monopole,  and  dipole  models  over  a  range  of  freestream  velocities  from  9  to  27  m/sec  are 
provided  in  Figures  4.7, 4.8,  and  4.9  respectively.  The  model  with  the  Josserand 
empirical  correlation  function  was  only  run  for  24  m/sec  (Figure  4.10).  The  Howe  model 
for  enhancement  of  a  natural  transition  near  the  leading  edge  of  a  half-plane  was  also  run 
for  the  24  m/sec  case  as  a  comparative  case  to  demonstrate  the  effect  in  the  change  to  the 
Green;s  function  (Figure  4.1 1).  All  of  the  models  are  compared  in  Figure  4.12  for  the  24 
m/sec  case. 

The  comparative  results  are  similar  to  those  seen  in  for  the  airflow  cases  but  are 
45  to  55  dB  higher  (after  accounting  for  26  dB  difference  in  reference  pressures)  for  the 
same  freestream  velocity. 

4.2  Role  of  Rise  Time,  t, 

With  a  factor  of  1  It?  in  the  amplitude  coefficient  of  the  power  spectrum,  the 
nondimensional  rise  time,  uj/dx,  (where  t,  is  defined  in  Figure  1.5)  has  usually  been  an 
independent  parameter  whose  value  is  chosen  to  get  the  model  results  to  match  the 
experimental  data.  It  has  taken  on  many  values  from  0.05  to  the  physically  implausible 
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1.36.  This  parameter  is  difficult  to  measure  and  would  seem  to  be  quite  dependent  on  the 
pressure  gradient.  However,  it  is  useful  to  look  at  some  values  which  have  been  extracted 
from  various  experimental  investigations  in  order  to  make  a  reasonable  first  estimate  for 
use  in  the  computational  code.  A  delta  wing  generated  vortex  as  seen  in  Emmons  spots 
produced  by  Littell  and  Eaton  (1991)  yielded  a  nondimensional  rise  time  of  0.05.  Use  of 
data  from  Huang  and  Shen  (1989)  gave  uj/Ax=Q.\%l.  Other  recent  investigations  on  a 
flat  plate  showed  a  value  of  0.30.  As  mentioned  previously,  Perraud  (1989)  found  best 
agreement  for  0.05<n,J/dx<0. 15.  Audet  et  al.  (1989)  used  a  value  of  1.36  which 
indicates  normal  velocity  fluctuations  continues  even  after  the  flow  is  fully  turbulent. 
Lauchle  (1981)  found  that  0.45  provided  the  best  agreement  with  the  Haddle  and 
Skudrzyk  (1969)  data.  It  appears  that  a  band  of  values  for  uj/Ax  around  0.20  will 
provide  a  reasonable  first  estimate.  Data  from  the  experiments  performed  by  Krane 
(1992)  in  the  fluid  mechanics  portion  of  the  ONR  project  indicate  that  for  an  isolated 
turbulent  spot  the  rise  time  is  O.ObcH^,  nJAx< 0.30  and  non-dimensional  fall  time  is 
0.13<wcf/ fal/Ax<OJ4.  Therefore,  uj/Ax  =  0.20 ,  was  used  for  each  of  the  models  for  the 
comparison  shown  in  Figure  4.13  for  the  airflow  condition  measured  in  the  experiment 
described  in  Chapters  5  and  6. 

4.3  Effects  of  Pressure  Gradient 

It  has  become  evident  that  transition  zone  noise  production  is  quite  dependent 
upon  the  pressure  gradient  over  the  transition  surface.  Narasimha  et  al.  (1984)  and 
Mautner  (1985)  found  that  a  positive  pressure  gradient  results  in  wall  pressure  fluctuation 
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spectral  shape  and  amplitude  very  close  to  those  for  a  zero  pressure  gradient  case. 
Mautner  found  that  for  an  adverse  pressure  gradient  that  the  nondimensional  wall 
pressure  power  spectrum  stayed  constant  to  a  nondimensional  frequency  of  f8*/U, -l 
rather  than  rolling  off  as  1//  after/<5*/l/j=0.1.  Mautner  concluded  that  this  broader 
frequency  forcing  function  would  result  in  a  greater  wall  pressure  field  coupling  with 
both  the  structure  and  the  acoustic  field  to  produce  more  intense  near  and  far  field  noise 
levels. 

Huang  and  Shen  (1989)  found  wall  pressure  levels  2  to  3  times  higher  for 
transition  in  adverse  pressure  gradients  compared  to  TBL  wall  pressure  levels.  Mautner 
and  Van  Atta  (1982)  found  rms  wall  pressures  for  individual  turbulent  spots  to  be  0.010- 
0.015  of  dynamic  pressure  (Vi  qU02)  for  zero  and  favorable  pressure  gradients  and  0.015- 
0.030  of  dynamic  pressure  in  adverse  gradients.  More  importantly,  in  the  adverse 
gradient,  the  wall  pressure  spectrum  was  nearly  constant  magnitude  with  increasing 
frequency  rather  than  the  typical  decreasing  slope.  This  is  enhanced  further  if  transition  is 
induced  by  laminar  boundary  layer  separation  and  re-attachment.  In  particular,  Huang  and 
Shen  found  that  wall  pressures  in  reattachment  zones  for  laminar  separation  were  10  to  15 
times  higher. 

The  models  all  have  a  sensitivity  to  the  length  of  the  transition  zone  and  the 
consequent  decrease  in  the  rise  and  fall  time  of  the  intermittency  indicator  function. 
Favorable  pressure  gradients  demonstrate  similar  or  slower  transition  than  a  zero  pressure 
gradient.  An  adverse  gradient  will  tend  to  accelerate  transition  and  shorten  the  transition 
zone  and  spot  growth.  Transition  in  an  adverse  pressure  gradient  will  result  in  higher 
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radiated  and  structural  excitation  levels  in  comparison  to  transition  in  a  zero  pressure 
gradient.  It  is  also  important  to  note  that  for  adverse  pressure  gradients,  it  is  very  possible 
that  some  amount  of  flow  separation  and  re-attachment  is  contributing  to  the  radiated 
noise. 


4.4  Effects  of  Scattering  by  a  Plane  Leading  Edge 

Howe’s  model  shows  the  impact  of  the  transition  zone  to  be  only  at  very  low 
frequencies  (below  about  200  Hz).  The  quadrupole  noise  normally  associated  with  the 
fully  turbulent  boundary  dominates  for  higher  frequencies  in  this  model  since  it  is  not 
restricted  to  the  transition  zone  alone.  The  effect  of  using  the  Josserand  intermittency 
function  is  to  decrease  the  radiated  noise  by  about  4  dB  (as  shown  in  Appendix  A.5). 

This  is  because  the  intermittency  is  used  to  multiply  the  point  pressure  spectrum  and  the 
Josserand  levels  are  slightly  lower  for  a  given  Z,  than  the  Farabee  derived  function.  In 
the  other  models,  it  is  the  Josserand  burst  rate  that  accounts  for  the  increased  levels 
calculated  when  the  Josserand  functions  are  used. 

4.5  General  Assessment  of  Models 

The  general  applicability  of  the  models  described  or  treated  in  this  dissertation  are 
assessed  as  follows: 

Fluctuating  Shear  Stress  Using  a  Li  gh  thill  Analogy  (GLPROG8Q)  -  Predicted 
noise  levels  are  much  too  low  (35  dB  below  other  models).  Even  the  use  of  the  Josserand 
intermittency  and  burst  frequency  functions  only  increases  the  results  by  12  dB. 
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Fluctuating  Displacement  Thickness  Using  a  Liepmann  Analogy.  Theoretical 
Space-Time  Correlation  Function  (GLPROG8D  -  The  premise  for  the  model  assumes  the 
growth  of  a  turbulent  spot  from  a  laminar  displacement  thickness  height  to  a  turbulent 
height  over  a  convection  time  in  the  transition  zone.  Based  on  the  highest  measured  burst 
rates  of  spots  being  between  25  %  and  60%  of  the  transition  zone  (Josserand  1986  and 
Farabee  et  al.,  1974),  predictions  have  appeared  to  be  best  using  a  normalized  rise  time, 
uct//Jx,  of  0.5  to  1.0.  Audet  et  al.  claim  a  best  fit  with  their  measurements  using 
normalized  rise  time  of  1.35  which  is  physically  implausible.  This  is  because  the 
monopole  model  of  the  sound  source  is  not  correct  and  a  higher  radiation  efficiency  is 
modeled  than  appropriate.  Thus,  slower  growth  times  must  be  used  to  compensate.  This 
is  where  the  substitution  of  the  dipole  description  of  vn(t)  yields  a  better  model  in  that,  at 
low  frequencies,  it  is  kd  cos  6  lower  than  the  levels  calculated  in  GLPROG81.  The  lower 
efficiency  dipole  is  necessary  to  allow  the  use  of  normalized  rise  times  of  0.20  (or  the 
variation  from  0.15  to  0.30  over  the  transition  region)  for  individual  spots  measured  by 
Krane  (1992)  in  an  open  water  channel.  This  dipole  is  implemented  in  a  variation  of 
GLPROG81  called  KMTRANS00.  For  k  d  >1,  the  sound  pressure  level  amplitudes 
predicted  by  GLPROG81  and  KMTRANS00  are  within  a  few  dB  of  each  other  because 
of  what  appears  to  be  two  sources  with  constructive  and  destructive  interference  and  a 
smearing  due  to  the  range  of  d  for  the  dipole  description  of  all  the  spots  in  the  transition 
zone.  The  prediction  by  KMTRANS00  matched  the  measured  spectral  shape  better  for 
frequencies  greater  than  200  Hz.  The  uncertainty  in  the  angle  of  the  dipole  axis  also  has  a 
great  effect  on  the  amplitude  of  the  prediction.  The  uncertainty  in  the  measurement  for 
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lower  frequencies  does  not  allow  a  meaningful  comparison. 

Fluctuating  Displacement  Thickness  Using  a  Liepmann  Analogy.  Empirical 
Space-Time  Correlation  Function  (JLTRANS89')  -  This  produces  prediction  levels 
slightly  higher  in  low  frequencies  than  GLPROG81  but  still  does  not  compute  well  in 
Mathcad.  While  it  might  provide  a  better  empirical  basis  for  modeling,  it  does  not  appear 
to  provide  enough  of  an  advantage  over  the  simpler  models  for  design. 

Two-Fluids  Model  -  Lagier  and  Somette  Formulation  -  This  model  has  two  source 
terms.  The  first  is  modification  of  TBL  source  levels  by  multiplication  with  the 
intermittency  factor.  In  effect,  this  makes  it  a  lesser  source  at  higher  frequencies  than  the 
nearby  (and  covering  a  much  greater  spatial  extent)  fully  developed  turbulent  boundary 
layer.  The  second  term  is  a  radiated  pressure  from  a  monopole  source  based  on  a 
Lighthill  analogy  that  is  effectively  the  same  as  the  fluctuating  baffled  piston-like  source 
due  to  the  Liepmann  analogy  used  by  Lauchle.  Their  formulation  provides  no  real 
improvement  over  GLPROG81,  and  is  computationally  more  complex. 

Boundary  leaver  Transition  on  Symmetric  Hydrofoils  -  Howe  Formulation 
(MSHTRANS99)  -  This  model  highlights  the  effect  of  the  scattering  and  diffraction  from 
a  nearby  leading  edge  on  the  low-frequency,  flow-induced  noise  spectrum  due  to  the 
transition  zone.  No  comparison  can  be  made  yet  against  wind  tunnel  or  water  tunnel 
measurements.  The  spectrum  derived  from  Josserand’s  empirical  space-time  correlation 
function  model  (Lauchle  and  Park,  2000)  could  be  substituted  for  the  Chase  spectrum 
model  for  the  calculations  covering  the  transition  zone.  The  impact  of  this  substitution  is 
not  known  but  is  not  expected  to  be  large.  The  low  frequency  amplitudes  may  be  low  due 
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to  only  using  the  intermittency  weighted  Chase  point  pressure  spectrum  and  not  including 
the  acoustic  pressures  measured  in  this  study  and  in  wall  pressures  in  a  transition  zone  in 
air  measured  by  Dufourq  (1984). 
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Figure  4.1  GLPROG80  predicted  radiated  pressure  levels  for  a  zero  pressure  gradient 
flat  plate  in  air  for  9  m/s  <  Um  <  27  m/s. 


90 


Sound  Pressure  Level  (dB//(20  uPaf/Hz/1m  Velocity  m/sec 


GLPR0G61  AIR  JJ 


GLPROG01  AIR  FF 


55 
50 
45 
40 
35 
30 
25 
20 
15 
10 

500  1000  1500  2000  2500  3000  3500  4000  4500  5000 
Frequency,  Hz 


FARABEE  Functions  (uct;  /Ax=0. 15)  JOSSERAND  Functions 


GLPROG81  AIR  FF  GLPROG81  AIR  JJ 

1 - — — ’ - . 80, - - — — -i . . 


-20 1 - - - * 

.  . 

_ _ 1 

1 

io1 

102  10* 
Frequency,  Hz 

104 

-20 1  * - - 

io1 

io2  io3 

Frequency.Hz 

io4 

Figure  4.2  GLPROG81  predicted  radiated  pressure  levels  for  a  zero  pressure  gradient 
flat  plate  in  air  for  9  m/s  <U„<21  m/s. 
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Figure  4.3  KMTRANSOO  predicted  radiated  pressure  levels  for  a  zero  pressure 
gradient  flat  plate  in  air  for  9  m/s  <U„<21  m/s. 


92 


dB  //  1  uPaA2/Hz  @rra 


Farabee  Intermittency  Function 


BL  Transition  Noise  For  Symmetric  Foil 


1  10  100  l  io3  i  io4 


Frequency  (Hz) 

“ “  Transition  Noise 
“  •  “  *  Quadrupole  Noise 
•••••  Total  Noise 


Figure  4.5  MSHTRANS99  predicted  radiated  pressure  levels  for  a  zero  pressure 
gradient  flat  plate  in  air  for  U„  =  24  m/s. 
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Figure  4.6  Comparison  of  models  for  radiated  pressure  levels  for  a  zero  pressure 
gradient  flat  plate  in  airflow,  U  =  24  m/s,zl<5*=0.243  cm,  Ax= 0.32  m, 
uJU m  =0.7,  r  =1.0 m. 
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Figure  4.7  GLPROG80  predicted  radiated  pressure  levels  for  a  zero  pressure  gradient 
flat  plate  in  water  for  3  m/s  <  U„  <  21  m/s. 
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Figure  4.8  GLPROG81  predicted  radiated  pressure  levels  for  a  zero  pressure  gradient 
flat  plate  in  water  for  3  m/s  <U„<  21  m/s. 


97 


Sound  Pressure  Level  (dB//(20  uPaf/Hz/l 


500  1000  1500  2000  2500  3000  3500  4000  4500  5000 
Frequency,  Hz 


500  1000  1500  2000  2500  3000  3500  4000  4500  5000 
Frequency,  Hz 


FARABEE  Functions  JOSSERAND  Functions 


KMTRANSOO  WATER  FF  KMTRANSOO  WATER  JJ 


Figure  4.9  KMTRANSOO  predicted  radiated  pressure  levels  for  a  zero  pressure 
gradient  flat  plate  in  water  for  3  m/s  <U„<  21  m/s. 
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Figure  4.11  MSHTRANS99  predicted  radiated  pressure  levels  for  a  zero  pressure 
gradient  flat  plate  in  water  for  U„  -  21  m/s. 
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Figure  4. 12  Comparison  of  models  for  radiated  pressure  levels  for  a  zero  pressure 

gradient  flat  plate  in  water,  Uj=2l  m/s,4<5*=0.021  cm,  4x=0.028  m, 
u/U„=0.7,  r  =1.0  m. 
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Figure  4. 13  Comparison  of  models  for  radiated  pressure  levels  for  a  zero  pressure 
gradient  flat  plate  in  airflow,  U=  21  m/s/lrf =0.278  cm,  Ax=0.9  m, 
uJU„  =0. 7,  r  =1.0  m. 
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CHAPTER  5 


EXPERIMENTAL  APPARATUS 

A  significant  effort  was  made  during  the  ONR  initiative  phase  of  this  investigation  to 
develop  a  facility  with  sufficiently  low  background  noise  which  maintained  laminar  flow 
over  at  least  the  first  third  of  a  test  plate.  The  components  of  the  entire  experimental 
apparatus  and  environment  include  the  ARL  flow-through  anechoic  chamber,  the  quiet 
flow  source,  the  flow  management  section,  nozzle,  diffuser,  and  the  various  test  plates 
that  have  been  created.  Further  information  on  the  facility  and  original  planned  research 
use  can  be  found  in  Marboe,  et  al.  (1991).  Subsequent  development  of  a  test  section 
located  on  the  floor  of  the  anechoic  chamber  by  Brungart  (1997)  utilizing  most  of  the 
Quiet  Airflow  Facility  components  appeared  to  be  much  closer  to  the  required 
background  noise  levels  to  conduct  the  desired  measurements.  This  facility  was  further 
enhanced  to  finally  create  a  quiet  enough  environment  for  the  measurements.  While  not 
all  of  the  Marboe,  et  al.  (1991)  apparatus  was  used  in  the  present  measurements,  it  is 
relevant  to  discuss  the  effort  that  went  into  trying  to  develop  a  facility  quiet  enough  to 
make  validation  measurements.  This  is  described  in  more  detail  in  Appendix  B.  The 
difficulties  encountered  do  emphasis  perhaps  how  insignificant  zero  pressure  gradient 
boundary  layer  transition  noise  is  when  compared  to  other  potential  noise  sources. 
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5.1  Quiet  Airflow  Facility 

This  facility  is  basically  a  quieted  air  source  connected  to  a  flow  management 
nozzle.  The  nozzle  can  form  an  open  jet  or  a  closed  entrance  into  a  test  section.  Several 
different  nozzles  exist.  One  has  a  semi-circular  ( 15.2  cm  radius)  orifice  situated  normal 
to  a  planar  baffle.  A  second  nozzle  has  a  square  orifice  18  inches  on  a  side.  The  current 
nozzle  has  a  15.2  cm  x  30.5  cm  rectangular  orifice.  The  entire  Quiet  Airflow  Facility  as 
finally  configured  for  this  investigation  is  shown  in  Figure  5.1. 

5.1.1  Flow-Through  Anechoic  Chamber 

A  low-frequency  acoustic  test  facility  to  study  the  acoustics  of  air  flow  and  air- 
moving  machinery  was  constructed  in  1986.  The  ARL  Penn  State  Flow-Through 
Anechoic  Chamber  was  designed  to  have  a  minimum  cut-off  frequency  of  70  Hz  (per  ISO 
standard),  using  fiberglass  wedges.  To  handle  the  air  flow  requirements  of  an 
experimental  Axial  Flow  Research  Fan,  the  air  inlet  wall  provides  100  square  feet  of  open 
area  to  insure  quiet  flow.  The  openings  are  uniformly  distributed  over  the  entire  wall 
opposite  the  test  facility  location  with  openings  omitted  in  the  wedge  blocks  adjacent  to 
the  side  walls  (Figure  5.2).  Two  layers  of  screen  and  perforated  sheet  metal  help  to  insure 
uniform  distribution  of  flow  with  minimum  turbulence.  The  exit  angle  between  the 
wedges  on  the  air  inlet  wall  was  designed  for  non-turbulent  flow.  This  final  wall  design 
was  established  following  a  scale  model  study  by  Quartararo  and  Lauchle  (1985). 

The  test  facility  opening  to  the  upper  control  room  is  placed  symmetrically  in  the 
chamber  so  that  flow  distortion  is  at  a  minimum.  Air  can  either  be  discharged  into  the 
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control  room  or  anechoic  chamber  with  flow  return  via  the  ceiling  plenum.  The  floor 
wedges  are  removable  to  gain  access  to  the  interior  of  the  test  volume  which  also  allows 
the  facility  to  be  used  as  a  hemi-anechoic  chamber.  A  modular  wire  mesh  platform  is 
available  for  mounting  experiments  within  the  volume. 

The  performance  characteristics  of  the  chamber  were  reported  in  Prout  and 
Marboe  (1990).  The  chamber  meets  the  requirements  of  both  ISO  3745  and  IEC  268 
(International  Standards  1972, 1977)  from  90  Hz  to  12.5  kHz.  Sound  pressure 
measurements  may  be  made  to  ISO  tolerances  within  1.0  to  4.0  meters  of  the  source  at 
frequencies  as  low  as  70  Hz.  Figure  5.3  shows  sound  level  meter  levels  plotted  on  a  set 
of  NCB  Balanced  Noise  Criterion  curves  (Beranek,  1988).  (These  are  revised  from  the 
PNC  curves  used  prior  to  1988).  The  empty  shell  noise  levels  were  approximately  NCB- 
26.  The  octave  band  background  noise  levels  in  the  anechoic  chamber  meet  a  balanced 
noise  criterion  of  NCB- 13  at  the  highest  octave  band  measured  (8000  Hz)  and  is  lower 
than  NCB- 10  in  frequency  bands  below  that.  Accordingly,  the  chamber  would  have  a 
rating  of  NCB- 13.  An  NCB- 15  level  had  been  predicted  for  the  chamber  when  the  wedge 
system  had  been  installed.  The  A-weighted  overall  level  of  the  ambient  noise  inside  the 
chamber  is  17  dB. 

5.1.2  Air  Source 

A  20  HP  centrifugal  blower  is  acoustically  and  mechanically  isolated  from  the 
wall  jet  apparatus  by  placing  it  in  a  sound  and  vibration  isolation  box  (Figures  5.4  and 
5.5)  that  is  located  in  the  control  room  adjacent  to  the  chamber.  The  inlet  to  the  blower  is 
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muffled  and  the  outlet  air  passes  through  a  treated  labyrinth  and  flexible  hose  to  the 
settling  chamber  of  the  wall  jet  facility  where  turbulence  management  screens  are  located. 
It  was  necessary  to  design  some  additional  sound  absorption  into  the  plenum  side  of  the 
sound  isolation  box  to  minimize  the  noise  traveling  down  the  flex  hose.  A  set  of  baffles 
was  designed  to  create  four  90  degree  turns  for  the  flow  so  that  higher  frequency  noise 
would  be  attenuated  by  the  Soundcoat  ®  2.5  cm  thick  foam-lead-foam  material  which 
was  used  to  line  the  box  and  the  baffles  as  shown  in  Figure  5.6.  This  treatment  plus  the 
intake  muffler  significantly  reduces  the  airborne  noise  from  the  blower. 

A  Buffalo  Forge  centrifugal  blower,  size  4  with  wheel  diameter  of  42  cm,  is  used 
as  the  prime  air  mover  (Figure  5.7).  It  is  powered  by  a  20  HP  AC  motor  driven  by  an 
ABB  inverter.  The  motor  is  water  cooled  to  decrease  the  build  up  of  heat  inside  the  foam 
lined  box.  Pillow  block  bearings  reduce  the  transmission  of  motor  induced  vibrations  to 
the  impeller.  A  rubber  bellows  attaches  the  blower  discharge  to  the  attenuation  plenum. 

A  Conic-Flow  Silencer  muffler  has  been  attached  to  the  inlet  to  the  blower  (Figure  5.5). 

A  capability  for  suction  for  flow  management  in  the  test  section  is  provided  by  the  blower 
intaking  air  from  a  10  cm  hose  attached  to  the  duct  between  the  muffler  and  the  rubber 
isolation  hose  into  the  blower.  A  throttle  can  be  mounted  on  the  muffler  intake  to  provide 
the  coarse  adjustment  of  the  suction  strength.  However,  neither  the  suction  nor  the 
throttle  capabilities  were  used  either  during  the  measurements  for  this  investigation. 

The  air  is  discharged  from  the  source  plenum  vertically  through  a  40  cm  diameter 
flexible  vinyl  hose  and  turned  with  a  120  degree  damped  sheet-metal  elbow.  A  section  of 
sheet-metal  40  cm  diameter  tubing  passed  through  the  portal  into  the  anechoic  chamber 
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(Figure  5.8)  before  being  turned  by  another  elbow  itself  into  another  Conic-Flow 
Silencer  muffler.  A  piece  of  flexible  hose  attached  the  silencer  outlet  to  a  diffuser  section 
for  entry  in  the  acoustic  settling  plenum  (Figure  5.9).  The  diffuser  is  lined  with  2.5  cm 
thick  Soundcoat.  The  plenum  is  lined  with  foam  acoustic  wedges.  The  wedges  are 
Illbruck  SONEXsuper  ™  15.24  cm  high  foam  with  3  wedges  in  a  30.5  cm  by  30.5  cm 
square.  Their  sound  absorption  characteristics  are  listed  in  Table  5.1.  A  square  bell 
mouth  was  constructed  out  of  split,  8  cm  PVC  pipe.  This  provides  a  smooth  inlet  into 
the  flow  management  section. 

Table  5. 1  Sound  Absorption  Coefficients  of  Foam  and  Wedges  Used  in  QAF 


Frequency  (Hz) 

SONEXsuper  Wedges 

Soundfoam  Embossed 

125 

0.18 

0.15 

250 

0.86 

0.45 

500 

1.22 

0.70 

1000 

1.24 

0.92 

2000 

1.14 

0.97 

4000 

1.20 

0.98 

5.1.3  Flow  Management 

The  preceding  plenum,  diffuser,  and  settling  chamber  with  flow  management 
screens  are  shown  in  Figure  5.10.  Ten  layers  of  screens  are  used  to  reduce  turbulence  and 
create  as  homogeneous  a  flow  as  possible.  The  entire  assembly  is  mounted  on  rubber 
pads  for  vibration  isolation.  Care  was  taken  to  make  all  comers  within  the  settling 
chamber  and  nozzle  as  smooth  as  possible  to  help  maintain  laminar  irrotational  flow. 
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A  new  nozzle  was  designed  with  a  greater  contraction  than  the  prior  square 
nozzle.  The  design  of  the  nozzle  followed  a  procedure  described  by  Morel  (1975).  The 
nozzle  is  shown  in  Figure  5.11.  The  design  calculation  is  described  in  more  detail  in  a 
Mathcad  document  in  Appendix  A  The  nozzle  pieces  were  rolled  into  the  correct 
curvature  using  templates,  welded  into  four  comer  sections  and  then  assembled. 

5.1.4  Test  Section 

The  test  section  consists  of  an  aluminum  test  plate  bottom  and  three  plexiglas 
sides  (Figure  5.12).  The  aluminum  plate  is  30.5  cm  wide  and  91.4  cm  long.  It  is  a  12.7 
mm  thick  sandwich  construction  with  damping  material  placed  in  the  middle.  A  series  of 
2.5  nun  holes  were  drilled  at  25.4  mm  intervals  along  the  centerline  of  the  plate  for 
insertion  of  a  hot  film  probe  and  turbulent  spot  trip  mechanisms  (see  Figure  5.13).  When 
the  tunnel  is  operating,  the  unused  holes  are  filled  with  wax. 

The  top  of  the  test  section  is  removable  and  can  be  adjusted  so  that  a  zero  pressure 
gradient  is  maintained  over  the  entire  test  plate.  The  top  is  shown  removed  in  Figure 
5.14.  The  opening  in  the  wedge  structure  for  air  recirculation  in  the  anechoic  chamber  is 
also  apparent.  The  top  of  the  test  section  has  a  61  cm  long  by  17.8  cm  wide  window  cut 
into  it  with  the  leading  edge  14  cm  downstream  of  the  leading  edge  of  the  test  section.  A 
0.0127  mm  (0.5  mils)  thick  piece  of  3M  mylar  film  was  stretched  over  the  entire  inner 
side  of  the  top.  A  spray  adhesive,  3M  Repositionable  75,  was  used  to  secure  the  mylar. 

A  hair  dryer  was  used  to  tighten  up  the  shrinkable  mylar.  One  of  the  sides  of  the  test 
section  has  been  modified  with  measurement  slots  so  that  hot  wire  probes,  stethoscope 
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burst  probes,  and  anemometer  probes  could  be  inserted.  The  forward  slots  were 
eventually  replaced  with  a  20.3  cm  long  by  10.2  high  window  cut  into  it  with  the  leading 
edge  5.1  cm  downstream  of  the  leading  edge  of  the  test  section.  It  was  covered  in  mylar 
in  a  manner  similar  to  above.  A  downstream  slot  was  located  at  85.75  cm  downstream  of 
the  leading  edge  of  the  test  section.  Typically,  the  inside  of  the  slot  would  be  covered 
with  clear  tape  and  a  small  hole  would  be  made  to  accommodate  the  desired  probe. 

By  surveying  over  the  plate  with  a  capillary  total  head  tube  on  a  stethoscope 
(Figure  5.15),  the  bursting  typical  of  transition  was  found  to  occur  near  the  trailing  edge 
of  the  test  and  working  forward  as  velocity  was  increased.  Side  wall  turbulent 
contamination  is  generally  thought  to  expand  at  an  angle  of  about  10  degrees  from  the 
comer  juncture  of  the  walls  of  the  test  section  (Motohashi  and  Blackwelder,  1983).  This 
contamination  should  join  across  the  test  section  nearly  at  the  end  of  the  test  section. 


5.1.5  Half-Plane 

A  half-plane  is  also  available  for  the  test  section.  The  splitter  plate  is  1.75  cm 
thick.  A  nose  portion  (aluminum)  is  the  leading  12.4  cm  of  the  plate.  It  is  attached  to 
another  71.3  cm  long  plate  (aluminum  anodized  blue)  as  shown  in  Figure  5.16.  The 
leading  edge  follows  design  guidance  from  Narasimha  and  Prasad  (1995)  using  a  nose 
length  of  3  times  the  plate  thickness.  The  super-elliptic  profile  is  given  by: 

=  1  ,  Q<,x<.a  n> 2  (5.1) 

where 
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a  length  of  the  nose 

b  half-thickness  of  the  plate 

n  was  2.5  for  this  application. 

A  diffuser  insert  can  be  attached  to  the  trailing  edge  of  the  plate  for  a  smooth  transition  to 
diffusion. 

5.1.6  Diffuser 

A  new  diffuser  needed  to  be  designed  to  work  with  the  rectangular  test  section. 
The  trick  in  diffuser  design  is  to  diffuse  the  flow,  slowing  it  down,  as  quickly  as  possible 
without  creating  flow  separation  on  the  diffuser  walls.  This  was  especially  necessary  in 
the  configuration  being  used  in  this  investigation  because  of  the  limited  distance  from  the 
test  section  to  the  wall  of  the  anechoic  chamber.  As  it  was,  the  discharge  was  purposely 
pointed  toward  the  air  outlet  wall  and  the  anechoic  wedges  were  removed  in  that  area. 
Soundcoat  material  was  placed  along  the  wall  of  the  chamber  at  a  slight  angle  to  reduce 
direct  refection  back  into  the  diffuser. 

The  diffuser  was  designed  following  a  method  described  by  Kline,  Abbott,  and 
Fox  (1959).  Details  can  be  found  in  section  B.5.  Even  though  the  design  methodology 
was  followed,  the  diffuser  walls  were  still  prone  to  a  transitory  stall  bubble  that  resulted 
in  unsteady  flow  discharging  from  the  diffuser.  The  final  configuration  created  an 
operating  point  just  on  the  edge  of  another  stall  avoidance  criterion  later  found  in  Barlow 
et  al.  (1999).  However,  that  unsteadiness  did  not  propagate  back  upstream  into  the  test 
section  based  on  the  hot  wire  measurements  near  the  trailing  edge  of  the  test  section. 
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The  diffuser  was  lined  with  2.5  cm  thick  embossed  Soundcoat  ®  to  absorb  sound 


propagating  down  the  diffuser  and  any  reflections  back  from  downstream.  The  lining  is 
shown  in  Figure  5.17.  The  absorption  properties  of  the  material  are  included  in  Table  5.1. 

With  the  mylar  window  in  the  top  of  the  test  section,  it  was  necessary  to  increase 
the  resistance  at  the  end  of  the  diffuser  to  increase  the  pressure  in  the  test  section  so  the 
mylar  was  not  drawn  down  into  the  test  section.  This  draw  down  introduced  significant 
asymmetry  in  the  flow  when  it  first  occurred.  The  test  requires  that  the  flow  on  the  top 
and  bottom  of  the  test  section  be  the  same  so  the  intermittency  measurements  could  be 
made  along  with  the  acoustic  measurements.  Four  1.2  cm  thick  fibrous  filters  were  used 
to  provide  this  resistance.  These  are  also  shown  in  Figure  5.17.  This  had  the  added 
advantage  of  eliminating  any  transitory  stall  in  the  diffuser. 

5.2  Instrumentation  and  Signal  Processing 

5.2.1  Hot  Wire  Anemometry 

5.2.1.1  Probe  Positioning 

Measurement  probes  are  positioned  in  the  test  section  with  a  three  dimensional 
Velmex  traversing  mechanism,  accurate  to  25.4  fi m,  through  two,  6.35  mm  wide  slots  cut 
in  the  plexiglass  tunnel  sidewall.  A  Mitutoyo  digital  caliper,  also  accurate  to  25.4  /tm, 
attached  to  the  vertical  stage  of  the  mechanism  to  facilitate  the  measurement  of  boundary 
layer  velocity  profiles.  Additionally,  the  vertical  stage  is  equipped  with  a  rotation  stage, 
thereby  allowing  a  measurement  probe  to  be  rotated  relative  to  the  test  section  axis  for 
added  ease  of  probe  positioning.  The  Velmex  traversing  mechanism  is  attached  to  a 
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precision  leveling  table  adjacent  to  the  wind  tunnel  test  section. 

5.2.1.2  Pitot-Static  and  Total  Head  Probe  Measurements 

The  mean  flow  velocity,  U,  is  determined  by  manipulating  the  Bernoulli  equation 
and  correcting  for  fluid  density  changes  due  to  temperature  variations.  The  relationship 
for  U  [m/s]  is 


where  P  atmospheric  pressure  in  psi 

T  temperature  ( degrees  Farenheit) 

V  voltage  relative  to  total  pressure 

V0  voltage  relative  to  static  pressure 

Although  a  pitot-static  probe  is  unsuitable  for  unsteady  velocity  measurements 
because  of  the  slow  response  of  the  fluid-filled  tubes  leading  to  the  differential  pressure 
sensor,  relatively  rapid  mean  pressure/velocity  measurements  were  made  with  the  use  of  a 
calibrated  Micro  Switch  differential  pressure  transducer  in  conjunction  with  a  DISA 
integrating  voltmeter.  The  maximum  error  due  to  voltmeter  uncertainty  is  estimated  to 
be  ±0.005  V. 

When  making  boundary  layer  velocity  measurements,  a  total  head  probe  with  an 
inner  diameter  of  0.0495  cm  was  used.  Additionally,  the  static  pressure  is  measured  for 
the  same  fluid  velocity  and  tunnel  streamwise  location  by  using  the  static  tap  of  a  pitot- 
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static  probe.  Pitot-Static  probes  may  be  used  for  velocity  measurements  only  when  at 
least  five  probe  diameters  separate  the  wall  and  the  probe.  This  avoids  the  creation  of  a 
venturi  passage  between  the  wall  and  the  probe.  Otherwise,  the  flow  would  be 
accelerated,  thereby  lowering  the  static  pressure  on  one  side  of  the  probe,  resulting  in 
inaccuracies  in  static  pressure;  consequently,  imprecise  fluid  velocities  would  be 
obtained.  When  conditions  permitted,  a  0.132  cm  diameter  pitot-static  probe  was  used. 

5.2.1.3  Hot  Wire  Anemometer  Measurements 

Both  mean  and  fluctuating  velocities  were  measured  with  a  two  channel  TSI IFA 
300  Research  Constant  Temperature  Anemometer.  The  anemometer  is  controlled  with 
TSI’s  ThermalPro  software  via  an  RS232  interface  with  a  Gateway  2000  486  DX2  66V 
personal  computer.  A  high  speed  A/D  converter,  mounted  internal  to  the  computer,  is 
used  to  convert  the  analog  data  output  from  the  anemometer  to  digital  data  for 
manipulation.  Hot-wire  data  are  sampled  at  a  frequency  of  10  kHz  with  a  low  pass  filter 
set  at  5  kHz.  Data  are  acquired  for  a  period  of  6.5536  seconds  resulting  in  64,000  points 
at  each  measurement  location.  Typically,  an  overheat  ratio  (the  ratio  of  operating 
resistance  to  actual  hot  wire  resistance)  of  1.72  was  used.  The  ThermalPro  software 
calculates  the  velocity  statistics  up  to  the  fourth  order  moment  (i.e.,  kurtosis).  All  data 
are  archived  on  Zip  disks  using  an  Iomega  Zip  100  drive  connected  to  the  host  computer. 

The  hot-wire  probes  are  visually  aligned  in  the  flow  field  and  positioned  as 
required  with  the  Velmex  traversing  mechanism.  A  TSI  Model  1218E  Boundary  Layer 
Probe  constructed  of  a  single  5  [im  diameter  platinum-coated  tungsten  wire  with  a 
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sensing  area  length  of  1.27  mm  was  used  for  streamwise  velocity  measurements. 

Single  component  wires  are  calibrated  in-situ  over  a  range  of  velocities  from  0.0 
m/s  to  25.0  m/s  with  pitot-static  probe  measurements  as  the  known  standard.  TSI’s 
ThermalPro  software  generates  a  calibration  curve  for  every  probe  calibration  and 
calculates  the  mean  square  error  which  is  usually  less  than  0.3.  Probe  velocity 
calibrations  are  checked  prior  to  every  profile  measurement  to  minimize  errors.  The 
experiment  temperature  is  checked  throughout  each  data  run  so  as  to  correct  for 
deviations  from  the  calibration  temperature.  Equations  (5.3)  and  (5.4)  illustrate  how  a 
small  change  in  temperature  can  cause  an  erroneous  fluid  velocity  to  be  determined  if  the 
temperature  departure  is  not  appropriately  accounted  for. 


E=Eh 

N 


(5.3) 


where 


Eb  wheatstone  bridge  voltage 

T,  hot-wire  sensor  operating  temperature 

Tc  hot-wire  sensor  calibration  temperature 

Te  hot-wire  sensor  experiment  temperature 


Veff=  K+AE+BE2+CE3+DE4  (5.4) 

where  Veff  =  effective  velocity  (velocity  component  normal  to  the  hot-wire) 

K,  A,  B,  C,  and  D  are  calibration  constants 
For  example,  if  the  experiment  temperature  where  to  unaccountably  increase,  the 
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effective  velocity  as  determined  from  Equation  (5.4)  would  be  lower  than  the  actual 
velocity. 

Because  hot-wire  probes  are  calibrated  against  pitot-static  probes  and  due  to  the 
sensitivity  of  hot-wire  measurements  to  temperature  changes,  pitot-static  probe  mean 
velocity  measurements  are  considered  more  accurate  than  hot-wire  measurements. 

5.2.1.4  Handheld  Anemometry 

Velocity  measurements  were  also  made  in  the  test  section  with  a  TSI  Model 
8345  VelociCalc  handheld  anemometer  (S/N  00030146).  It  has  an  accuracy  of  3%  of 
scale  (or  ±  0.015  m/sec,  whichever  is  greater)  although  the  calibration  was  less  than 
1.7%.  This  anemometer  reads  from  0  to  30  m/sec  with  a  resolution  of  0.01  m/sec.  It  was 
used  to  verify  flow  conditions  and  temperature  in  the  test  section. 

5.2.2  Sound  Pressure 

The  single  channel  sound  pressure  measurements  used  for  diagnostics  are  made 
with  a  Brilel  &  KjaerType  4135  Vi  inch  microphone  (serial  number  1946536)  with  Type 
2670  preamplifier  (serial  number  2003873).  A  Brilel  &  Kjaer  Type  5935L  power  supply 
and  amplifier  (serial  number  2016632)  is  used  with  linear  weighting  set.  The  amplifier  is 
capable  of  -5  to  +55  dB  gain.  It  is  battery  powered  for  low  noise.  If  necessary,  a  Krohn- 
Hite  3342  2  channel  filter  (serial  number  216)  is  used  as  a  high  pass  filter  to  remove  very 
low  frequency  facility  noise.  This  is  then  input  to  the  DSP  Technologies  SigLab  signal 
analyzer  described  in  section  5.2.5.  The  procedures  for  use  of  the  acoustic  intensity 
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methods  are  well  documented  in  Fahy  (1989). 


5.23  Acoustic  Intensity 

The  two  channel  sound  intensity  measurements  are  made  with  a  Bruel  &  Kjaer 
Type  4178  Intensity  Vi  inch  microphone  set  consisting  of  Type  4135  microphone  (serial 
number  1946506)  with  Type  2670  preamplifier  (serial  number  2003871)  for  channel  1 
and  Type  4135  microphone  (serial  number  1946464)  with  Type  2670  preamplifier  (serial 
number  2003873)  for  channel  2.  Otherwise,  the  signal  processing  equipment  is  the  same. 
Table  5.2  provides  a  full  list  of  equipment  used  for  the  time-averaged  intensity  method. 

The  microphones  are  generally  placed  outside  of  the  mylar  windows.  The 
acoustic  transparency  of  this  material  was  established  by  Kargus  and  Lauchle  (1997). 

The  insertion  loss  is  generally  less  than  3  dB  with  deviations  less  than  1  dB  from  300  Hz 
to  3  KHz. 

5.23.1  Real-Time  Analog  Intensity  Circuitry 

A  real-time  active  acoustic  intensity  circuit  has  been  developed  to  provide  a  signal 
that  can  be  processed  along  with  a  hot  film  shear  stress  sensor  to  yield  a  coherent  output 
intensity  similar  to  coherent  output  power  processing.  Figure  5.19  provides  a  process 
diagram. 

The  electronic  circuit  was  designed  and  fabricated  by  Todd  Batzel  of  ARL  Penn 
State.  It  is  described  in  detail  in  Appendix  C  including  circuit  diagrams.  The  gains  in  the 
circuitry  are  multiply  by  250  at  the  input  filter,  multiply  by  295  in  the  integrator, 


116 


and  divide  by  10  in  the  multiplier.  This  gives  a  total  gain  multiplier  of  7375 
or  38.7  dB.  A  high  pass  filter  has  been  implemented  with  a  cut-off  at  30  Hz.  The  filter  is 
a  fourth  order  Bessel  filter  with  linear  phase.  This  was  chosen  for  this  particular 
application  in  order  to  remove  some  low  frequency  tunnel  noise  that  otherwise  would 
have  used  up  about  40  dB  of  dynamic  range.  A  reset  button  is  installed  for  the  integrator 
to  allow  re-zeroing  in  case  of  a  built  up  bias.  The  circuit  can  be  powered  by  two  onboard 
9V  batteries  or  an  AC  adapter. 

5.23.2  Calibration 

Two  plane  wave  tubes  were  fabricated  for  calibration  of  the  two  microphone  intensity 
probe.  They  cover  frequency  ranges  800-5270  Hz  and  70-950  Hz,  respectively.  They  are  PVC 
pipe  with  aluminum  ends  for  probe  mounting.  They  are  shown  in  Figure  5.18.  Further 
details  are  provided  in  Table  5.3.  The  individual  microphones  are  also  calibrated  with  a 
BrUel  &  Kjaer  Type  4220  pistonphone  (serial  number  347321). 
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Table  5.2  Equipment  List  for  Time- Averaged  Acoustic  Intensity  Measurement. 


Component 

Channel  1 

Channel  2 

Microphone 

Brtiel  &  Kjaer  4178  pair 
(phase  matched  4135  mic) 
S/N  1946464 

3.67  mV/Pa 

Brtiel  &  Kjaer  4178  pair 
(phase  matched  4135  mic) 
S/N  1946506 

3.55  mV/Pa 

Preamplifier 

Brtiel  &  Kjaer  Type  2670 
Falcon®  Range  w  /  LEMO 
connectors 

S/N  2003871 

Brtiel  &  Kjaer  Type  2670 
Falcon®  Range  w/  LEMO 
connectors 

S/N  2003873 

Power  Supply  /  Amplifier 

Brtiel  &  Kjaer  Type  5935L 

2  channel 

S/N  2016632 

Same 

Digital  Signal  Processor 

DSP  Technologies 

SigLab  20-42 

S/N  11204 

Same 

Computer 

Dell  Latitude  LM  133 

S/N  6533L 

Microsoft  Windows  95 
SigLab  3.24 

Matlab  5.3 

Same 

Table  5.3  Calibration  Tube  Characteristics. 
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5.2.4  Coherent  Output  Intensity 

The  signal  processing  diagram  for  the  calculation  of  coherent  ouput  intensity  is 
provided  in  Figure  5.18.  The  output  of  the  real-time  active  intensity  processing  circuit  is 
input  to  channel  1  of  the  SigLab  analyzer.  The  output  of  the  TSI  anemometer  is  input  to 
channel  2  of  the  SigLab  analyzer.  The  coherence  is  computed  and  multiplied  along  with 
microphone  sensitivity  and  other  constants  with  channel  1  to  yield  the  desired  coherent 
output  intensity.  This  is  all  covered  in  more  detail  in  Appendix  C. 

5.2.4.1  Hot  Film  Sensor 

The  hot  film  signals  are  measured  with  a  TSI  Incorporated  model  1750  constant 
temperature  anemometer  (serial  number  2318).  A  TSI  probe  model  1237-10  (serial 
number  984071)  is  inserted  flush  into  holes  in  the  test  plate.  A  31  ohm  control  resistor  is 
necessary.  The  holes  were  spaced  25.4  mm  apart  longitudinally  along  the  centerline  of 
the  test  plate.  The  other  holes  were  waxed  over. 

5.2.4.2  Conditioning  with  a  Shear  Stress  Sensor 

It  is  desired  to  apply  the  power  of  coherent  output  processing  to  the  acoustic 
measurements.  In  this  case,  it  is  assumed  that  shear  stress  due  to  passage  of  a  turbulent 
spot  is  caused  by  the  same  physical  phenomenon  (i.e.  growth  of  the  spot)  that  results  in 
acoustic  radiation.  Therefore,  there  should  be  coherence  between  those  signals.  The 
noise  associated  with  each  sensor  is  assumed  to  be  independent.  Hot  films  have  been 
shown  to  be  insensitive  to  acoustic  noise.  The  coherence  is  calculated  using  Equation 
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(5.5)  using  the  shear  stress  signal  and  the  intensity  circuit  signal. 


ICV,WI; 

a,xV) 


(5.5) 


5.2 .5  Other  Processing  Hardware  and  Software 

The  analog  signals  are  processed  by  a  DSP  Technologies,  Inc.  SigLab  Model  20- 
42  dynamic  signal  analyzer.  This  has  four  channels  of  0-20  kHz  analysis.  Each  channel 
contains  a  fourth  order  quasi-elliptic  filter  with  liner  phase  for  anti-aliasing.  The  signal 
path  include  sigma  delta  analog-to-digital  converters,  an  ADSP-2105  decimating  digital 
filter,  and  a  Texas  Instruments  TMS320C31  digital  signal  processor.  The  dynamic  range 
is  greater  than  90  dB.  A  SCSI  interface  connects  it  to  a  Dell  Latitude  LM  133  notebook 
PC.  The  SigLab  has  a  hardware  performance  verification  procedure  that  calibrates 
channel  input  and  output  gain  and  offsets  which  are  written  to  the  EEPROM. 

The  SigLab  3.24  software  is  written  in  The  Mathworks,  Inc  MATLAB  5.3.  Plots 
can  be  made  on  a  Hewlett  Packard  DeskJet  870Cxi  inkjet  printer.  The  anemometer  and 
the  microphone  amplifier  are  shown  in  Figure  5.20.  The  rest  of  the  signal  processing 
instrumentation  is  shown  in  Figure  5.21. 
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Table  5.4  Equipment  List  for  Coherent  Output  Intensity. 


Component  Channel!^^j^happgj 2 

Microphone  Briiel  &  Kjaer  4178  pair  Briiel  &  Kjaer  4178  pair 

(phase  matched  4135  mic)  (phase  matched  4135  mic) 

S/N  1946464  S/N  1946506 

3.67  mV/Pa _ 3.55  mV/Pa _ 

Preamplifier  Briiel  &  Kjaer  Type  2670  Briiel  &  Kjaer  Type  2670 

Falcon®  Range  w/  LEMO  Falcon®  Range  w/  LEMO 

connectors  connectors 

S/N  2003871  S/N  2003873 _ 

Power  Supply  /  Amplifier  Briiel  &  Kjaer  Type  5935L  Same 

2  channel 

S/N  2016632 _ 

Analog  Filter  and  Krohn-Hite  3342  Same 

Amplifier  2  channel 

S/N  216 _ 

Real  -  Time  Processor  ARL  Penn  State  Same 

2  channel 

Real-Time  Active  Intensity 
Circuit 

Shear  Stress  Conditioning  TSI  Model  1237-10  Hot 

Signal  Film 

S/N 984071 

Model  1151  Probe  Holder 

Anemometer  TSI  Constant  Temperature 

Anemometer  Model  1750 
S/N 2318 _ 

Digital  Signal  Processor  DSP  Technologies  Same 

SigLab  20-42 

S/N 11204 _ 

Computer  Dell  Latitude  LM 133  Same 

S/N6533L 

Microsoft  Windows  95 
SigLab  3.24 

Matlab  5.3  _ _ _ _ 
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Figure  5. 1  Quiet  Airflow  Facility  configuration  for  present  experiment. 
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Figure  5.2  Flow-through  anechoic  chamber. 
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Figure  5.4  Blower  and  sound/vibration  isolation  box  . 
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Figure  5.9  In-flow  silencer  and  diffuser  into  acoustic  settling  chamber. 
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Figure  5.10  Flow  management  plenum,  diffuser,  screens,  and  settling  sections  prior  to 
nozzle. 
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Figure  5.13  Test  section  view  showing  drilled  holes  for  hot  film  probe  or  spot 
generation  excitation. 
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Figure  5. 14  View  of  nozzle,  plenum,  and  diffuser. 
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Figure  5.16  Half-plane  plate  desig: 


for  test  section. 
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Figure  5.17  Looking  upstream  through  diffuser  to  the  test  section  and  the  added 
diffuser  resistance. 
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Figure  5.19 


Processing  diagram  for  coherent  output  intensity. 


Figure  5.20  Anemometer  and  microphone  amplifier 


Figure  5.21  Test  instrumentation. 
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CHAPTER  6 


MEASUREMENT  RESULTS 


6.1  Test  Section  Flow 

The  test  section  flow  was  characterized  at  two  longitudinal  locations  in  the  test 
section  using  hot  wire  anemometry.  These  were  x  =  7.6  cm  and  77.5  cm  of  the  test 
section.  The  forward  measurement  was  used  to  establish  that  there  was  uniform  laminar 
entering  the  test  section.  This  is  shown  in  Figure  6.1  for  two  tunnel  velocities  and 
compared  to  the  Blasius  boundary  layer  profile.  The  uniformity  across  the  test  section 
area  is  seen  in  Figure  6.2.  The  hot  wire  measurements  were  made  traversing  across  the 
test  section  at  heights  of  25, 50,  and  75%  of  the  vertical  dimension  of  the  test  section. 

The  variation  is  less  than  1%.  It  should  be  noted  that  these  measurement  were  conducted 
prior  to  the  addition  of  the  resistance  on  the  end  of  the  diffuser  which  was  done  during 
the  acoustic  measurements.  Thus,  the  top  velocity  in  the  test  section  is  higher  than  could 
be  achieved  for  the  acoustic  measurements. 

The  downstream  measurement  location  was  used  to  establish  a  zero  pressure 
gradient  flow  condition.  Spacers  were  placed  under  the  support  arms  on  the  test  section 
top  until  the  proper  uniformity  of  flow  was  achieved.  A  similar  survey  uniformity  was 
completed  at  this  location.  Those  results  are  shown  in  Figure  6.3. 

A  15  cm  strip  of  6.3  mm  mesh  screen  placed  near  the  leading  edge  of  the  test 
section  was  used  to  trip  the  boundary  layer  in  order  to  check  the  fully  developed  turbulent 
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flow  at  the  end  of  the  test  section.  Figure  6.4  shows  a  comparison  of  the  boundary  layer 
profiles  and  turbulence  intensities  for  this  test  and  similar  data  taken  with  a  separate 
nozzle  and  test  section  in  the  Quiet  Airflow  Facility.  The  flow  behavior  appears  normal 
with  this  nozzle  and  test  section. 

6.2  Intensity  Measurement  Calibration 

A  calibration  sequence  was  performed  with  the  intensity  instrumentation.  The 
configuration  for  a  free-field  calibration  is  shown  in  Figure  6.5.  The  SigLab  signal 
generator  was  used  to  provide  a  2.5  Vrms,  10  kHz  bandwidth  white  noise  input  to  a  25 
cm  woofer  speaker  mounted  in  a  wooden  baffle.  The  intensity  probe  was  positioned 
above  the  speaker  81  cm.  Reflections  were  minimized  by  use  of  13  cm  thick  fiberglas 
pads. 

Before  the  speaker  test  was  done,  the  microphones  were  amplitude  calibrated 
using  the  pistonphone.  A  magnitude  and  phase  mismatch  test  was  performed  in  the  small 
calibration  tube  (Figure  6.6).  The  microphones  were  extended  into  the  tube  to  make  sure 
that  the  pressure  compensation  holes  on  the  preamplifiers  were  exposed.  This  is  the 
reason  for  the  notable  spikes  around  2  kHz.  This  corresponds  to  a  quarter  wavelength  for 
reflection  from  the  end  wall  of  the  tube  back  to  the  microphone.  The  phase  error  is 
generally  less  than  0.4  degrees  and  the  magnitude  error  is  less  than  0.5  dB.  For  this 
instrumentation  phase  error,  the  intensity  error  is  0  to  -1  dB  down  to  60  Hz.  The 
microphones  have  a  25.4  mm  spacing.  This  spacing  produces  an  accuracy  within  1.0  dB 
over  a  range  of  320  to  2500  Hz.  Above  this,  the  error  increases  to  5  dB  at  5000  Hz. 
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There  are  also  bias  errors  associated  with  the  type  of  source  field  and  the 
measurement  distance  (Fahy,  1989).  This  bias  error  is  shown  for  point  monopole  and 
dipole  fields,  and  a  measurement  distance  of  5  cm  in  Figure  6.7.  The  effect  is  very  small 
below  3  KHz.  All  of  the  mentioned  biases  apply  in  the  same  manner  to  each  of  the 
measurements  conditions,  so  relative  amplitude  differences  are  still  meaningful.  The 
random  error  was  reduced  by  taking  256  averages  for  each  measurement  which  yields  a 
±0.4  dB  random  error  for  the  cross-spectrum  used  in  the  analysis. 

The  results  of  the  speaker  calibration  verification  is  shown  in  Figure  6.8.  The  P-P 
intensity  is  compared  to  the  far  field  intensity  determined  by  P,2/pc.  The  agreement  is 
very  good  down  to  about  150  Hz  where  the  near-field  errors  increase  and  produce  the 
amplitude  fluctuations  seen  in  the  figure.  Above  3.5  kHz,  the  discrepancy  is  due  to  finite- 
difference  approximations.  The  test  section  top  with  the  mylar  window  was  then  placed 
between  the  speaker  and  the  intensity  probe.  The  results  in  Figure  6.8b  show  generally 
good  agreement  with  narrowband  deviations  being  more  pronounced  with  the  mylar. 
Figure  6.9  shows  the  same  data  processed  using  50  Hz  constant  bandwidth  and  1/3  octave 
bandwidths.  This  shows  agreement  within  1  dB  over  the  frequency  range  of  interest  from 
250  Hz  to  3  kHz. 

The  value  of  the  use  of  an  intensity  probe  (with  its  cosine  directivity)  over  a  single 
microphone  is  shown  in  Figure  6.10.  A  5-12  dB  improvement  in  signal-to-noise  ratio  is 
achieved.  The  equivalent  pressure  level  converted  from  the  intensity  level  is  compared  to 
the  pressure  level  from  only  microphone  1  of  the  intensity  probe  for  two  measurement 
locations. 
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6.3  Flat  Plate  Natural  Transition  Zone 


Acoustic  active  intensity  measurements  were  made  for  a  natural  transition  zone 
occurring  on  the  top  surface  of  the  test  section.  The  tunnel  was  run  at  its  highest  possible 
velocity  condition  of  21  m/sec.  In  this  case,  the  transition  began  just  downstream  of  the 
mylar  window  leading  edge.  The  downstream  edge  of  the  mylar  window  was  at  about 
65%  of  the  transition  zone.  It  is  within  this  range  that  the  predominant  portion  of  the 
radiated  acoustics  has  been  predicted  to  occur  according  to  Audet  et  al.(1989a,b)  as 
shown  in  Figure  2.2.  Any  further  extension  of  the  mylar  window  would  have  been 
unpractical  due  to  a  weakening  of  its  structure  and  higher  likelihood  of  surface  dynamics. 
A  longitudinal  survey  was  conducted  over  the  mylar  window  with  the  probe  orientation 
normal  to  the  window  as  shown  in  Figure  6.1 1. 

Hot  film  measurements  were  made  by  inserting  the  hot  film  through  predrilled, 
but  waxed  over,  holes  in  the  bottom  of  the  test  section.  It  is  presumed  from  the  uniform 
flow  measured  previously  that  the  boundary  layer  growth  and  transition  on  the  bottom  of 
the  test  section  matches  that  on  top  of  the  test  section.  A  comparison  is  made  in  Figure 
6.12  of  the  locations  and  amplitudes  of  the  measured  intermittency  to  the  Josserand 
intermittency  function. 

The  intermittency  was  determined  using  the  circuitry  used  by  Josserand  (1986) 
which  turns  turbulent  spot  shear  fluctuations  into  a  one  or  zero  volt  indicator  that  is 
integrated  to  show  an  intermittency.  The  match  with  the  Josserand  function  for  a  90  cm 
transition  zone  is  very  good.  The  intensity  measurement  locations  are  also  indicated  on 
the  graph.  The  Josserand  intermittency  values  at  those  locations  were  then  used  for  any 
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subsequent  intermittency  weighted  calculations. 

Recall  from  Figure  6.8  that  the  intensity  probe  showed  narrowband  fluctuations  at 
lower  frequencies  that  were  not  present  in  the  far  field  processed  output.  These 
fluctuations  are  due  to  tunnel  modes  with  their  characteristic  constructive  and  destructive 
interference.  Even  though  sources  on  the  bottom  of  the  test  section  should  be  at  least  12 
dB  lower  than  similar  sources  on  the  top  of  the  test  section  as  received  at  the  intensity 
probe,  due  to  a  certain  mode,  the  source  could  seem  to  be  located  near  the  top.  The  true 
spectrum  should  be  smoother  as  for  the  far  field  in  Figure  6.8.  So  further  processing  was 
done  to  smooth  the  spectra  by  using  50  Hz  constant  bandwidth  processing.  All  of  the 
presentations  of  measurements  will  therefore  be  shown  in  both  1.6  Hz  and  50  Hz 
processing  bandwidths.  Figure  6.13  shows  a  comparison  of  these  processing  bandwidths 
to  a  1/3  octave  processed  representative  measurement. 

The  active  intensity  measured  at  the  7  longitudinal  positions  are  shown  in  Figure 
6.14.  The  flow  condition  at  x=20  cm  is  laminar.  Selective  spectra  are  shown  in  Figure 
6.15  for  the  laminar  location  and  the  transition  locations  near  what  is  predicted  to  be  the 
peak  for  acoustic  pressure. 

A  second  experiment  was  conducted  with  a  boundary  layer  trip  consisting  of  a  6 
cm  long  strip  of  50  grit  sandpaper  located  at  x  =  2  cm  which  is  upstream  of  the  mylar 
window  (Figure  6.16).  The  measured  active  intensity  spectra  are  shown  in  Figures  6.17 
and  6.18.  The  intent  was  to  demonstrate  that  there  was  signal  to  noise  available  for 
measurement  of  laminar  and  turbulent  signatures  and  that  they  were  distinctive  from  the 
transitional  structure  radiation. 
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A  third  experiment  was  conducted  to  try  to  reduce  the  length  of  the  transition  zone 
by  stimulating  the  boundary  layer  to  transition  earlier.  Several  variations  of  sandpaper 
grits  and  sizes  were  tried  to  get  a  condition  that  was  determined  by  stethoscope  to  be 
intermittent  rather  than  turbulent.  Finally,  two,  10  mm  long  strips  of  150  grit  sandpaper 
were  located  at  x  12  and  15  cm.  The  results  arc  provided  in  Figures  6.19  and  6.20.  Even 
with  the  stimulation,  the  results  are  not  significantly  different  from  those  of  natural 
transition. 

Comparisons  are  made  for  each  measurement  locations  between  these  three  test 
cases  in  Figures  6.21  through  6.27.  The  spectra  for  natural  transition  tend  to  be  less  than 
the  turbulent  spectra  in  frequencies  above  1  kHz  by  2-3  dB  except  for  the  Z,  =  0.45 
position  where  the  transition  levels  exceed  the  fully  turbulent  measured  levels.  This 
correlates  well  with  the  wall  pressures  measured  by  Dufourq  (1984)  shown  in  Figure 
6.28.  Recall  from  Equation  2.17  in  the  description  of  the  two  fluids  model  that  the  wall 
pressure  is  modeled  as  an  intemittency  weighted  wall  pressure  spectrum  due  to  turbulence 
plus  an  acoustic  pressure  due  to  the  spot  dynamics  in  intermittency.  The  acoustic 
measurements  described  here  are  the  first  validation  of  that  hypothesis  since  the  active 
acoustic  intensity  is  the  radiating  component  and  not  unduly  influenced  by  the  near  field 
reactive  intensity.  At  the  next  Z,  positions,  the  turbulent  and  transitional  levels  are  more 
equivalent  similar  to  the  wall  pressure  measurement.  At  these  locations  and  beyond,  the 
quadrupole  turbulent  noise  production  is  dominating  the  higher  order  spot  related 
radiation  due  to  coalescence. 

The  transition  intensity  measurements  were  then  summed,  but  weighted  by  the 
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longitudinal  extent  of  the  area  covered  by  that  measurement  in  order  to  get  one  transition 
intensity  spectrum.  Thus  it  still  reflects  pW/m2  which  can  be  converted  to  the  pressure 
spectrum  using  the  far  field  assumption  for  comparison  to  the  computational  model 
predictions.  All  models  were  exercised  for  the  QAF  measurement  condition  of  U  =  21 
m/sec.  Ax  =  90  cm,  and  A8*=0.278  cm.  These  results  are  compared  to  the  mesured 
radiated  pressure  levels  in  Figure  6.29.  As  expected,  the  displacement  thickness  model 
with  the  dipole  source  and  the  Farabee  intermittency  and  burst  rate  functions  provides  the 
best  estimate  to  the  measured  from  300  Hz  to  3000  Hz  ( the  range  for  which  bias  and 
random  errors  are  within  2  dB). 

Both  the  transition  and  tripped,  fully  turbulent  measurements  were  summed  based 
on  an  area  weighting  of  power  to  get  a  representative  zone  acoustic  power  signature.  The 
turbulent  spectra  were  weighted  at  each  measurement  location  for  the  local  intermittency 
representative  of  that  area.  These  spectra  are  compared  in  Figure  6.30.  For  frequencies 
above  1250  Hz,  the  turbulent  spectra  is  about  3  dB  above  the  transition  levels.  This 
would  be  expected  for  the  frequencies  where  turbulent  boundary  layer  noise  would  tend 
to  dominate.  However,  for  lower  frequencies,  that  are  predicted  to  be  dominated  by 
transition  flow  structure  acoustic  sources,  the  transition  levels  exceed  the  turbulent  levels 
3  to  5  dB.  This  confirms  that,  locally,  the  transition  zone  noise  affects  low  frequency 
signature  and  has  a  source  order  higher  than  quadrupole. 
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Figure  6. 1  Boundary  layer  profile  for  U  =  23. 1  m/s,  x  =7.6  cm,  y=2.54  cm  from 

centerline. 
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Figure  6.3  Spanwise  uniformity  of  the  test  section  streamwise  velocity  for 

U=23.1  m/s,  x  =  77.5  cm,  at  25,  50,  and  75%  of  test  section  height. 
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Figure  6.4  Comparison  of  boundary  layer  profiles  and  turbulence  intensity  profiles 
tripped,  fully  turbulent  flow  for  nominal  U=23.9  m/s  near  end  of  test 
section. 
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Magnitude  and  phase  of  microphone  mismatch  as  measured  in  small 
calibration  tube. 
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degrees 


Intensity  Estimate  Error,  dB 


Point  Monopole  Field 
Point  Dipole  Field 


Figure  6.7  Intensity  estimate  bias  errors  due  to  point  monopole  and  dipole  source 
fields  for  microphone  separation  of  2.54  cm. 
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Speaker  Calibration 


Frequency  (Hz) 


Frequency  (Hz) 


Figure  6.8  Comparison  of  (a)  P-P  time-averaged  intensity  measurement  to  far  field 
intensity  at  0.8m  above  speaker,  and  (b)  P-P  time-averaged  intensity 
measurement  with  and  without  mylar  to  open  far  field  intensity  at  0.5m 
above  speaker  in  1.6  Hz  bands. 
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Figure  6.9  Comparison  of  P-P  time-averaged  intensity  measurement  with  and 
without  mylar  to  open  far  field  intensity  at  0.5m  above  speaker  in 
(a)  50  Hz  bands  and  (b)  1/3  octave  bands. 
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Figure  6.10  Intensity  probe  gain  over  a  single  microphone. 
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Figure  6. 12  Locations  and  levels  of  hot  film  measurements  for  intermittency 

compared  to  Josserand  intermittency  function  and  locations  of 
intensity  measurements. 
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Frequency  (Hz),  bands 


Figure  6.13  Examples  of  spectrum  smoothing  achieved  through  wider  band 
processing. 
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Active  Intensity  (dB  re  1  pW/m2),  1.6  Hz  band  Active  Intensity  (dB  re  1  pW/m2),  50  Hz  band 


21  m/sec,  natural  transition 


Frequency  (Hz) 

21  m/sec,  natural  transition 


Figure  6. 14  Active  intensity  in  1 .6  and  50  Hz  processing  bands  measured  for  natural 

transition  region  of  21  m/s  in  zero  pressure  gradient  above  mylar  wall 
for  -0.03<  Zx  <  0.60  . 
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21  m/sec,  natural  transition 


Frequency  (Hz) 


Figure  6.15  Active  intensity  in  1 .6  and  50  Hz  processing  bands  measured  for  natural 
transition  region  of  21  m/s  in  zero  pressure  gradient  above  mylar  wall 
for  Zx  = -0.03, 0.45  and  0.53. 
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Figure  6.16  Boundary  layer  trip  of  50  grit  sandpaper  in  6  cm  long  strip  at  x=  20  cm 
and  boundary  layer  stimulation  of  150  grit  sandpaper  in  two  10mm  long 
strips  at  x=  12  and  15  cm . 
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21  m/sec,  turbulent  tripped  flow 
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Figure  6. 17  Active  intensity  in  1 .6  and  50  Hz  processing  bands  measured  for  tripped 
turbulent  flow  at  21  m/s  in  zero  pressure  gradient  above  mylar  wall 
for  -0.03<  Z„  <  0.60 . 


21  m/sec,  2  strip  stimulated  transition 


Frequency  (Hz) 
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Figure  6. 19  Active  intensity  in  1 .6  and  50  Hz  processing  bands  measured  for 

stimulated  transition  region  at  21  m/s  in  zero  pressure  gradient  above 
mylar  wall  for  -0.03<  Zx  <  0.60  . 
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Figure  6.20  Active  intensity  in  1.6  and  50  Hz  processing  bands  measured  for 

stimulated  transition  region  at  21  m/s  in  zero  pressure  gradient  above 
mylar  wall  for  Zx  =  -0.03, 0.45  and  0.53. 
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Frequency  (Hz),  1.6  Hz  bands 


Figure  6.21  Comparison  of  active  intensity  in  1 .6  and  50  Hz  processing  bands 

measured  for  natural  and  stimulated  transition  and  tripped  turbulent  flow 
at  21  m/s  in  zero  pressure  gradient  above  mylar  wall  for  Zx  =  -0.03. 
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Figure  6.22  Comparison  of  active  intensity  in  1.6  and  50  Hz  processing  bands 

measured  for  natural  and  stimulated  transition  and  tripped  turbulent  flow 
at  21  m/s  in  zero  pressure  gradient  above  mylar  wall  for  Zx  =  0.09. 
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Frequency  (Hz),  1.6  Hz  bands 

Figure  6.23  Comparison  of  active  intensity  in  1.6  and  50  Hz  processing  bands 

measured  for  natural  and  stimulated  transition  and  tripped  turbulent  flow 
at  21  m/s  in  zero  pressure  gradient  above  mylar  wall  for  Zx  =  0.20. 


21  m/sec,  x  =  50  cm 


Frequency  (Hz) 


Frequency  (Hz),  1 .6  Hz  bands 


Figure  6.24  Comparison  of  active  intensity  in  1.6  and  50  Hz  processing  bands 

measured  for  natural  and  stimulated  transition  and  tripped  turbulent  flow 
at  21  m/s  in  zero  pressure  gradient  above  mylar  wall  for  Zx  =  0.31. 
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Figure  6.25  Comparison  of  active  intensity  in  1.6  and  50  Hz  processing  bands 

measured  for  natural  and  stimulated  transition  and  tripped  turbulent  flow 
at  21  m/s  in  zero  pressure  gradient  above  mylar  wall  for  Zx  =  0.45. 
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Figure  6.26  Comparison  of  active  intensity  in  1.6  and  50  Hz  processing  bands 

measured  for  natural  and  stimulated  transition  and  tripped  turbulent  flow 
at  21  m/s  in  zero  pressure  gradient  above  mylar  wall  for  Zx  =  0.53. 
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Figure  6.27  Comparison  of  active  intensity  in  1.6  and  50  Hz  processing  bands 

measured  for  natural  and  stimulated  transition  and  tripped  turbulent  flow 
at  21  m/s  in  zero  pressure  gradient  above  mylar  wall  for  Zx  =  0.60. 
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O  Normalized  pressure  due  to  radiated  acoustics  (70-5000  Hz) 
measured  in  current  investigation  in  air  (21  m/sec). 

Figure  6.28  Fluctuating  wall  pressure  measurements  made  in  air  by  Dufourq  (1984). 
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(Measured  data  presented  with  50  Hz  processing  but  corrected  to 
1  Hz  spectral  level) 

Figure  6.29  Comparison  of  transition  noise  model  predictions  with  measured 
directly  radiated  sound  pressure  levels. 
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21  m/sec,  transition  zone  acoustics 
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Figure  6.30  Comparison  of  intermittency  weighted  measured  turbulent  boundary 
layer  noise  with  measured  transition  zone  directly  radiated  sound 
pressure  levels. 


CHAPTER  7 


CONCLUSIONS  AND  FUTURE  WORK 


7.1  Conclusions 

The  research  described  in  this  dissertation  has  resulted  in  the  following  unique 

contributions  and  conclusions. 

1.  First  unambiguous  direct  measurement  of  the  acoustic  radiation  from  a  boundary 
layer  transition  zone. 

2.  Confirmed  that  the  wall  pressure  peak  measured  by  Dufourq  (1984)  in  the  range 
0.4  <  Z*  <  0.55  under  a  transition  zone  is  due  to  the  radiated  acoustic  pressure. 

3.  The  acoustic  theory  developed  by  Lauchle  using  the  Liepmann  fluctuating 
displacement  thickness  approach,  and  an  infinite  flat  plate  Green’s  function  with 
monopole  source  field  appears  to  be  incorrect  when  compared  against  the 
measurements.  Lauchle’ s  model  has  been  improved  by  using  a  dipole  source  field 
based  on  the  Krane  (1992)  water  channel  data,  and  on  qualitative  observations  of 
the  flow  structures  present  during  transition.  The  true  axis  angle  of  the  dipole 
relative  to  the  wall  has  yet  to  be  determined.  There  may  even  be  an  occurrence  of 
multiple  angles  in  a  natural  transition  zone. 

4.  Rather  than  using  a  constant  dipole  axis  (source  separation  distance),  which 
results  in  unrealistic  high-frequency  content  of  the  predicted  transition  noise 
(extensive  sinusoidal  null  variation  in  levels),  a  statistical  estimate  is  constructed 
by  treating  the  dipole  axis  as  a  random  variable.  This  variable  is  assumed 
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uniformly  distributed  over  physically  predictable  limits.  This  results  in  a 
smoothed  spectral  estimate  and  better  match  with  the  measured  data.  This 
prediction  model  provides  the  best  fit  to  the  measured  data  for  frequencies  above 
300  Hz.  Because  of  measurement  uncertainty,  a  proper  comparison  cannot  be 
made  below  that  frequency  at  this  time.  In  addition,  the  angle  of  the  dipole  may 
vary  depending  on  the  growth  dynamics  of  the  turbulent  spot.  A  treatment  similar 
to  the  one  described  for  dipole  source  separation  distance  in  section  3.1.2  could 
also  be  incorporated  for  this  dipole  axis  angle. 

3.  Demonstrated  that  the  intermittency- weighted,  fully  turbulent  boundary  layer 

radiated  sound  power  level  is  a  good  approximation  to  the  directly  radiated  sound 
power  level  from  the  transition  zone.  The  level  was  found  to  be  about  3  dB  below 
that  of  the  transition  zone  noise  for  f  <  1  kHz  .and  3  dB  high  for  f  >  1  kHz. 

6.  The  two-fluids  model  proposed  by  Lagier  and  Somette  offers  no  real  prediction 
benefit  over  the  Lauchle  model  once  the  spatial  dependence  is  brought  back  inside 
the  integral  and  given  the  simplifications  that  were  necessary  to  make  the  two 
fluids  model  integrable. 

7.2  General  Observations  on  Transition  Noise 

The  following  are  observations  that  have  resulted  from  the  cumulative  knowledge 
gained  over  the  course  of  this  research  including  the  background  materials  and 
understanding  of  the  previous  works  of  others. 

1.  Most  of  the  focus  in  boundary  layer  transition  acoustics  in  the  last  eight  years  has 
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been  on  developing  direct  numerical  simulation  (DNS)  approaches  to  modeling 
using  spatial  and  temporal  approaches  to  solving  the  complete  Navier-Stokes 
equations.  Based  on  the  results  of  Wang,  et  al.  (1996),  there  appears  to  be  more 
“noise”  in  the  calculation  than  is  produced  by  the  actual  fluid  dynamics.  However, 
the  value  of  DNS  may  really  be  to  provide  insight  on  the  necessary  modeling  of 
substructure  dynamics  that  influence  the  spatial  definition  of  the  normal  velocity 
used  in  the  fluctuating  displacement  thickness  model  and  contribute  in  their  own 
right  to  the  acoustic  pressure  particularly  influencing  the  near-field.  This  bridges 
the  gap  to  some  extent  between  the  two  fluids  model  using  the  Lighthill  analogy 
and  the  Liepmann  approach  used  for  the  fluctuating  displacement  thickness 
model. 

2.  The  predicted  low  frequency  contribution  of  transition  noise  seems  to  be 
significant,  particularly  for  f  <300  Hz  and  higher  mean  flow  velocities.  Audet  et 
al.  (1989a,  Figure  10)  show  this  as  a  domination  of  the  low  polar  order  growth 
mechanism  of  the  spot  (a  la  Lauchle)  over  the  quadrupole  source  due  to  the  TBL 
statistics  in  the  turbulent  spots.  Howe’s  model  (1999)  for  the  effect  of  a  transition 
zone  near  the  leading  edge  of  a  hydrofoil  clearly  show  an  enhancement  of  the 
same  low  frequency  regime  due  to  the  transition  zone. 

3.  Laboratory  measurements  performed  prior  to  this  dissertation  indicate  that 
transition  zone  noise  levels  may  be  contaminated  by  other  sources.  It  is  suspected 
that  one  contaminating  source  for  early  measurements  in  this  investigation  (with 
an  early  QAF  configuration)  was  trailing  edge  noise  when  a  jet  was  discharging 
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across  a  flat  plate.  However,  the  axisymmetric  body  tests  by  Perraud  (1989) 
indicated  propagation  of  acoustic  noise  from  the  transition  region  to  sensors  in  the 
laminar  flow  zone.  The  results  reported  here  qualitatively  supports  this  finding. 

4.  Transition  in  an  adverse  pressure  gradient  will  result  in  higher  radiated  and 
structural  excitation  levels  in  comparison  to  transition  in  a  zero  pressure  gradient. 
Huang  and  Shen  (1989)  found  wall  pressure  levels  2  to  3  times  higher  for 
transition  in  adverse  pressure  gradients  compared  to  TBL  wall  pressure  levels. 
Mautner  and  Van  Atta  (1982)  found  rms  wall  pressures  for  individual  turbulent 
spots  to  be  0.010-0.015  of  dynamic  pressure  (Vi  eU02)  for  zero  and  favorable 
pressure  gradients  and  0.015-0.030  of  dynamic  pressure  in  adverse  gradients. 
More  importantly,  in  the  adverse  gradient,  the  wall  pressure  spectrum  was  nearly 
constant  magnitude  with  increasing  frequency  rather  than  the  typical  decreasing 
slope.  This  is  enhanced  further  if  transition  is  induced  by  laminar  boundary  layer 
separation  and  re-attachment.  In  particular,  Huang  and  Shen  found  that  wall 
pressures  in  reattachment  zones  for  laminar  separation  were  10  to  15  times  higher. 
Hwang  and  Shang  (1989)  and  Hwang  and  Maidanik  (1990)  have  worked  on 
approaches  for  modeling  the  coupling  of  these  wall  pressures  within  transition 
zones  with  the  underlying  structure. 

5.  The  preceding  observation  is  perhaps  what  should  command  our  continued 
interest  because  of  the  engineering  relevance  associated  with  forebodies  of  bodies 
of  revolution  and  for  airfoil  leading  edges  where  acoustic  sources  may  be  located. 
Just  as  trailing  edge  flows  are  critical  to  the  direct  radiation,  excitation,  re- 
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radiation  and  scattering  of  noise  at  airfoil  trailing  edges  (see  Brooks  and  Hodgson, 
1981  and  Gershfeld,  1997),  so  may  transition  near  the  leading  edge  of  an  airfoil 
act  as  a  noise  source.  Of  course,  this  is  quite  dependent  upon  the  chordwise 
distribution  of  thickness  for  the  airfoil  and  the  local  radii  of  curvature  (Howe, 
1981). 

7.3  Recommended  Future  Work 

The  recommendations  for  future  work  are  as  follows. 

•  CFD  Validation  -  More  experimental  data  will  be  needed  at  increasing  Reynolds 
numbers  to  validate  the  computational  fluid  dynamics  solution  approaches  using 
Reynolds  Averaged  Navier-Stokes  (RANS)  or  Large  Eddy  Simulation  (LES) 

•  Transition  Noise  in  Adverse  Pressure  Gradients  -  Radiated  noise  and  fluid 
mechanics  measurements  need  to  be  made  for  a  systematic  variation  of  adverse 
pressure  gradient.  Pursuit  of  further  modeling  or  measurement  of  zero  pressure 
gradient  flows  does  not  seem  warranted  based  on  the  small  contribution  to  the 
overall  flow-induced  radiated  noise  for  a  flat  plate.  However,  the  indications  are 
that  adverse  pressure  gradients  should  shorten  the  transition  zone  and,  thus,  result 
in  a  higher  radiated  noise  component. 

•  Transition  on  Blades  -  Radiated  noise  due  to  boundary  layer  transition  on  air-  and 
hydrofoils  is  a  topic  of  continued  interest.  The  usefulness  of  the  models  described 
in  this  thesis  when  coupled  with  the  proper  Green’s  function  should  be  evaluated 
against  measurements.  But,  it  is  noted  that  the  mechanism  for  transition  on  blades 
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is  often  a  bypass  phenomenon  or  a  result  of  separation. 

Experimental  Validation  -  Continue  to  utilize  the  current  Quiet  Airflow  Facility  to 
measure  the  radiation  from  the  leading  edge  of  the  existing  a  half-plane  foil  in  the 
presence  of  various  pressure  gradients,  particularly  adverse.  Also,  perform  similar 
experiments  on  axisymmetric  headforms. 

Coherent  Output  Intensity  -  This  procedure  should  continue  to  be  pursued  with  a 
better  implementation  of  the  real-time  intensity  circuitry.  The  concept  of  using  a 
sensor  unaffected  by  the  acoustic  field  to  condition  the  intensity  measurement  is  very 
powerful. 
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APPENDIX  A 


MATHCAD  MODEL  FILES 


This  appendix  contains  the  worksheets  for  several  of  the  boundary  layer  transition 
noise  models  and  experiment  diagnostic  routines  as  developed  in  the  Mathcad  2000 
Professional  computer  environment.  Mathcad  ®  is  a  product  of  Mathsoft,  Inc.  of 
Cambridge,  MA.  This  product  is  documented  in  a  user’s  guide  and  reference  manual 
(1999). 
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A.1  Fluctuating  Shear  Stress  Using  a  Lighthill  Analogy 
Mathcad  file:  GLPROG80.MCD 
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RADIATED  NOISE  DUE  TO  FLUCTUATING  SHEAR  STRESS  FROM 
BOUNDARY  LAYER  TRANSITION 


[GLPROG80.MCD,  Mathcad  2000  Professional] 

This  program  is  the  implementation  of  a  boundary  layer  transition  noise  model 
described  in: 

Lauchle,  G.  C.,  1980,  ’On  the  radiated  noise  due  to  boundary- layer  transition,!/. 

Acous.  Soc.  ofAmer.,  67,  pp.  158-168. 

ASSUMPTIONS 

-  zero  pressure  gradient  flow  and  natural  transition  flow  structures,  i.e.  no  separation). 

-  flow  is  over  an  infinite,  planar,  rigid  surface 

-  based  on  very  low  Mach  number  -  use  only  dipole  components  due  to  fluctuating 
shear  stress 

-  flow  is  statistically  homogeneous  in  spanwise  direction,  X3,  non-homogeneous  in 
streamwise  direction,  xi,  stationary  but  Poisson  distributed  in  time  so  we  can  use  a 
classical  correlation  function  for  a  random  sequence  of  bursts,  exp(-2N|t|) 

-  neglect  quadruples  and  octupoles  due  to  fluctuating  Reynolds  stresses  due  to  lower 
radiation  efficiencies  -  however,  their  amplitudes  relative  to  shear  stresses  may 
overcome  this  lower  efficiency. 

-  the  space-time  correlation  function  for  the  indicator  function  <H’>is  direction 
separable 

-  the  lateral  term  integration  derived  from  this  separation  is  equal  to  2L3(n  1)  =  2(1.5 

t| ,  tan  a)  which  is  about  r\  ,/4  for  a=9.6  degrees  (angle  of  spreading  for  turbulent  spot) 

-  Uc/Uo  is  independent  of  frequency 
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DEFINITIONS 


Ax  is  the  length  of  the  transition  zone  in  feet  or  meters 

A5*  is  is  the  change  of  the  boundary  layer  displacement  thickness  from  laminar  to 
turbulent  states  in  feet  or  meters  (A5*  =  5*turbulent  -  5*laminar)  (A5star  in 
worksheet) 

Uo  is  the  freestream  velocity  in  feet/sec  or  meters/sec 

Uc  is  the  convection  velocity  in  feet/sec  or  meters/sec  (nominally  Uc  =  0.7*Uo) 
a*  is  the  flare  parameter  that  accounts  for  moving  axis  decorrelation 
c  is  the  speed  of  sound  in  feet/sec  or  meters/sec 
Ret  is  the  transition  Reynolds  number  based  on  Ax 

Rexo  is  the  Reynolds  number  based  on  the  distance  to  the  start  of  the  transition  zone 
W  is  the  spanwise  width  of  the  transition  zone  in  feet  or  meters 

v  is  the  kinematic  viscosity  in  fee^/sec 

p  is  the  fluid  density  in  slugs/feet3 
r  is  the  far-field  observation  distance  in  feet  or  meters 
0  is  the  angle  from  the  plane  of  the  wall  to  the  observation  point 
o  is  the  estimate  of  the  fluctuating  wall  shear  stress 
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SPECIFY  VALUES  FOR  FLOW  AND  GEOMETRY  VARIABLES 


[Note:  Mathcad  is  operating  in  SI  units  mode  for 
all  calculations.  The  user  may  enter  values  for 
variables  in  any  appropriate  measurement  system 
so  long  as  the  units  are  specified.  The  program 
will  automatically  make  the  conversion  to  SI  units 
during  calculation.] 

Define  1/3  Octave  band  center  frequencies 
a  :=  0.. 9 

Freqa  :=  Fra-1.0-Hz 
b  :=  10..  19 

Freqb  :=  Fr(b-io)10.0-Hz 
d  :=  20..  29 

Freqj  :=  Fr^O)- 100.0- Hz 
e  :=  30..  39 

Freqe  :=  Fr(e_30)- 1000.0- Hz 
e  :=  40.. 49 

Freqe  :=  Fr(e_40)- 10000.0- Hz 


Constant  Air  Fresh  Water 

Constants  for  user c 

reference.  v  (m2/s)  1.51  x  10~5  1.01  x  10'6 

p  (kg/m3)  1.21  999 

Pref(pPa)  20  1 


Salt  Water 

1487 

1.01  x  10‘6 

1027 

1 


Ax  :=  90.0- cm 
Ads  tar  :=  0.278- cm 


Uo  :=  21- 
sec 


W  :=  .305- m 


r  :=  1-m 


0  :=  0-deg 


m 

c  :=  343 - 

sec 


v  :=  1.51-10  5  — 


Uc  :=  0.7- Uo 
Ret  :=  0.8- 106 


Pref  :=  20-10  6-Pa 
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POWER  SPECTRAL  DENSITY 


This  section  will  calculate  some  factors  which  will  be  used  in  the  leading  terms  for  the 
final  spectrum  calculation. 


Xo  :=  v -  calculated  transition  point 

Uo 


Rexo(Xo)  :=  Uo-— 


transition  point  Reynolds  number 


-_1_ 

o(Xo)  :=  0.0288-p-Uo2-Rexo(Xo)  5 


(  — 1 

vl-  11.53- Rexo(Xo)  10  } 


fluctuating  wall 
shear  stress 


This  section  will  define  the  intermittency  function,  y(Z),  and  the  turbulent  spot 
burst  rate,  N(Z),  where  Z  =  qx/Ax.  To  lessen  computational  burden  later,  the 
integrations  will  use  another  function  B(Z)  =  N(Z)(Ax/Uo).  Obviously,  this  just 
removes  the  factor  (Uo/Ax)  from  the  definition  of  N(Z)  but  this  factor  is  common 
to  both  models  of  interest  and  is  accounted  for  in  the  overall  equation.. 


INTERMITTENCY  FACTOR 

This  is  the  time-averaged  value  of  the  indicator  function  and  represents  the  percentage 
of  time  that  the  flow  is  in  a  turbulent  state.  Either  of  two  models  for  this  factor  can  be 
chosen  by  the  user  for  the  computation  of  the  spectrum.  The  empirical  relation  from 
Farabee,  et  al.  (1974)  was  used  in  Lauchle  (1980). 

yFarabee(Z)  :=  1  -  exp(-4.185-Z2)  Empirical  relations  from  Farabee,  et  al.  (1974). 

yJL(Z)  :=  1  -  exp[-(l  +  3.4-Z)-Z2~J  Empirical  relations  from  Josserand  and 

Lauchle  (1986)  which  was  based  on  the  work 
of  Gedney  (1979)  asessing  both  a  Dirac  line 
source  hypothesis  and  a  constant  source 
hypothesis  for  intermitency. 

♦♦CHOOSE  INTERMITTENCY  FACTOR:  y  (Z)  :=  yFarabee(Z) 
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TURBULENT  SPOT  BURST  RATE 


This  is  the  time-averaged  number  of  turbulent  bursts  sensed  per  second.  Either  of  two 
models  for  this  factor  can  be  chosen  by  the  user  for  the  computation  of  the  spectrum. 
The  empirical  relation  from  Farabee,  et  al.  (1974)  was  used  in  Lauchle  (1980). 


BFarabee(Z)  :=  1.272  Z  exp(-4.185  Z2) 
BMarboe(Z)  :=  1.8  Z15  exp(-4  Z2) 


BJL(Z)  :=  2.38- 


(1 


_Y(Z)).ln  _ 


_1 _ 

-Y(Z)) 


**CHOOSE  SPOT  BURST  RATE  FUNCTION:  B  (Z)  :=  BFarabee(Z) 

FLARE  PARAMETER 

This  parameteris  used  in  the  longitudinal  space-time  correlation  function  to  account 
for  moving  axis  decorrelation  due  to  c  aburst  not  remaining  perfectly  correlated  as 
the  distance  increases  beyond  some  characterisitc  scale. 

Q 

astarAx(Z)  :=  1  +  83.35- Z  Normalized  flare  parameter  used  in  the 

longitudinal  space-time  correlation  function 
developed  by  Lauchle  (1980)  using  data  fron 
Shivitz  (1975). 
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F  FUNCTION 


This  section  will  define  the  F  function  given  by  Lauchle(1980)  equation  37  which 
is  integrated  over  the  transition  zone.  It  is  assumed  tha  Uc/Uo  is  independent  of 
frequency.  The  spectrum  is  then  calculated  for  single  frequencies  which  represent 
1/3  octave  center  frequencies  from  10  Hz  to  80  KHz. 


T1(KD,Z) 


2-astarAx(Z) 

__  _ 

(ostarAx(Z)2  +  KD2) 2 


T2(KD,Z)  :=  exp(-astarAx(Z)-z)-sin 


( 


KDZ  -  atan 


V 


^astarAx(Z)^ 
V  135  „ 


T3(KD,Z)  :=  exp[-astarAx(Z)-(l  -  Z)]  sin 


KD-(l-Z)-atan 


^  ocstarAx(Z)  ^ 
KD 


Func  (KD ,  VEL ,  Z) 


B(Z)-y(Z)»Z-(T1(KD,Z)  -T2(KD,Z)  +  T3(KD,Z)) 

\_ 

[[4-(B(Z))2-VEL'2]  +  KD2](  txstarAx  (Z)2  +  KD2) 2 


F(A,C) 


Func(A,C,Z)dZ 


Equation  implemented  below 
with  trapezoidal  integration. 
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SPECTRAL  CALCULATIONS 


Define  frequency  range 

fmax  :=  5000- Hz 
points  :=  1000 
nn  :=  1 ..  points  +  1 


ORIGIN  :=  1 


KL  :=  1000  k  :=  I ..  KL  fk:=fmax-10 


^KL  , 


Znn  :_ 


nn  -  1 
points 


intk,nn  :=  Func 
1 


(  Ax  Uc  ^ 

2-TC-fk - , —  ,Znn 

Uc  Uo  J 


dZ  := 


points 


FFk  •=  Qntk,nn+1  +  intk,nn)  ^ 


nn  =  1 


Definition  of  function  representing  power  spectral  density: 


Fk  :=  I  2-n-ftf— 
1  Uc 


\2 


FFk 


Gk  := 


[w-(o(Xo))2-(cos(e))2Uo-(Ax)2]-Fk 


»  2  2  2 
2-r  -n  c 


Calculation  of  power  spectral  density: 


SPECGk  :=  10- log 


Gk 


(Pref) 


Hz 
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RADIATED  NOISE  DUE  TO  SHEAR  STRESS  FROM  BOUNDARY  LAYER 
TRANSITION  USING  LIGHTHILL  ANALOGY 

GLPROG80.MCD 


Ax  =  0.9  m 
A5star  =  0.00278  m 

Uo  =  21msec’1 

Uc  =  14.7  msec’1 

Ret  =  800000 


W  =  0.305  m  c  =  343  msec  1 

r  =  lm  v  =  0.0000151m2  sec'1 

9  =  0  P  =  1-21  kgm’3 

Pref  =  0.00002kgm  1  sec  2 


Y(Z)  =  yFarabee(Z)  B(Z)  =  BFarabee(Z) 
SPECTRAL  PLOTS 
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MODIFICATION  HISTORY 

Original  7  May  1998  Conversion  of  elements  of  GLPROG  VAX 
program 

Mod  1  14  May  1998  Addition  of  choices  for  intermittency  and  burst  rate 

functions 

Mod  2  4  July  2000  Substitution  of  trapezoidal  integration. 

Spectrum^  :=  f-sec 
Spectrum^  :=  SPECG 


Q 

Zi\80QAFFF21.TXT 

Spectrum 
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A.2  Fluctuating  Boundary  Layer  Displacement  Thickness  Using  a  Liepmann 
Analogy  -  Theoretical  Space-Time  Correlation  Function 

Mathcad  file:  GLPROG81.MCD 
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DEFINITIONS 


Ax  is  the  length  of  the  transition  zone  in  feet  or  meters 

AS*  is  is  the  change  of  the  boundary  layer  displacement  thickness  from  laminar  to 
turbulent  states  in  feet  or  meters  (AS*  =  8*turbulent  -  5*laminar)  (ASstar  in 
worksheet) 

Uo  is  the  freestream  velocity  in  feet/sec  or  meters/sec 

Uc  is  the  convection  velocity  in  feet/sec  or  meters/sec  (nominally  Uc  =  0.7*Uo) 
a*  is  is  the  flare  parameter  accounts  for  moving  axis  decorrelation 
c  is  the  speed  of  sound  in  feet/sec  or  meters/sec 
Ti  is  the  time  it  takes  S*iaininarto  transition  to  dS^^j  in  sec 

rise  is  the  non-dimensional  rise  time  for  the  AS*  (rise  =  Uc*Ti/Ax) 

Ret  is  the  transition  Reynolds  number  based  on  Ax 

Rexo  is  the  Reynolds  number  based  on  the  distance  to  the  start  of  the  transition  zone 
W  is  the  spanwise  width  of  the  transition  zone  in  feet  or  meters 

v  is  the  kinematic  viscosity  in  fee^/sec 

p  is  the  fluid  density  in  slugs/feet3 

r  is  the  far-field  observation  distance  in  feet  or  meters 

0  is  the  angle  from  the  plane  of  the  wall  to  the  observation  point 
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SPECIFY  VALUES  FOR  FLOW  AND  GEOMETRY  VARIABLES 


(  1.0  > 
1.25 
1.6 
2.0 
2.5 
3.15 
4.0 
5.0 
6.3 


[Note:  Mathcad  is  operating  in  SI  units  mode  for 
all  calculations.  The  user  may  enter  values  for 
variables  in  any  appropriate  measurement  system 
so  long  as  the  units  are  specified.  The  program 
will  automatically  make  the  conversion  to  SI  units 
during  calculation.] 

Define  1/3  Octave  band  center  frequencies 
a  :=  0.. 9 

Freqa  :=  Fra-1.0-Hz 


b  :=  10..  19 

Freqb  :=  Fr(b_10)-10.0-Hz 
d  :=  20.. 29 

Freqd  :=  Fr(d-20)*  100.0- Hz 
e  :=  30..  39 

Freqe  :=  Fr(e_30)- 1000.0- Hz 
f  :=  40.. 49 

Freqf  :=  Fr(f_4Q)- 10000.0- Hz 


Constant- 
s  for  user 
reference. 


Constant  Air  Freshwater  Saltwater 

c  (m/s)  343  1525  1487 

v  (m2/s)  1.51  x  10'5  1.01  x  10’6  1.01  x  10’6 

p  (kg/m3)  1.21  999  1027 

Pref  (pPa)  20  1  1 


Ax  :=  2.82-cm  Ti  :=  0.45-sec 

ASstar  :=  0.021-cm  rise  :=  0.15 


Uo 


21.0- 


m 


sec 


W  :=  .305- m 


c 

V 


1525  - 
sec 


1.01-10"6 


2 

m 

sec 


Uc  :=  0.7-Uo 
Ret  :=  1.5-106 


r  :=  1.0-m 


m 


Pref  :=  1-10  6-Pa 
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POWER  SPECTRAL  DENSITY 


This  section  will  define  the  intermittency  function,  y(Z),  and  the  turbulent  spot 
burst  rate,  N(Z),  where  Z  =  rix/Ax.  To  lessen  computational  burden  later,  the 
integrations  will  use  another  function  B(Z)  =  N(Z)(Ax/Uo).  Obviously,  this  just 
removes  the  factor  (Uo/Ax)  from  the  definition  of  N(Z)  but  this  factor  is  common 
to  both  models  of  interest  and  is  accounted  for  in  the  overall  equation.. 


INTERMITTENCY  FACTOR 

This  is  the  time-averaged  value  of  the  indicator  function  and  represents  the  percentage 
of  time  that  the  flow  is  in  a  turbulent  state.  Either  of  two  models  for  this  factor  can  be 
chosen  by  the  user  for  the  computation  of  the  spectrum.  The  empirical  relation  from 
Farabee,  et  al.  (1974)  was  used  in  Lauchle  (1980). 


yFarabee(Z)  :=  1  -  exp(-4.185-Z2)  Empirical  relations  from  Farabee,  et  al.  (1974). 

yJL(Z)  :=  1  -  exp[-(l  +  3.4-Z)-Z2]  Empirical  relations  from  Josserand  and 

Lauchle  (1986)  which  was  based  on  the  work 
of  Gedney  (1979)  asessing  both  a  Dirac  line 
source  hypothesis  and  a  constant  source 
hypothesis  for  intermitency. 


♦♦CHOOSE  INTERMITTENCY  FACTOR:  y  (Z)  :=  yFarabee(Z) 
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TURBULENT  SPOT  BURST  RATE 


This  is  the  time-averaged  number  of  turbulent  bursts  sensed  per  second.  Either  of  two 
models  for  this  factor  can  be  chosen  by  the  user  for  the  computation  of  the  spectrum. 
The  empirical  relation  from  Farabee,  et  al.  (1974)  was  used  in  Lauchle  (1980). 


BFarabee(Z)  :=  1.272-Z-exp(-4.185-Z2) 
BMarboe(Z)  :=  1.8-ZL5exp(-4-Z2) 


BJL(Z)  :=  2.38- 


♦♦CHOOSE  SPOT  BURST  RATE  FUNCTION:  B  (Z)  :=  BFarabee(Z) 

FLARE  PARAMETER 

This  parameter  is  used  in  the  longitudinal  space-time  correlation  function  to  account 
for  moving  axis  decorrelation  due  to  a  burst  not  remaining  perfectly  correlated  as  the 
distance  increases  beyond  some  characteristic  scale. 

cts  tar  Ax  (Z)  :=  1  +  83.35-Z^  Normalized  flare  parameter  used  in  the 

longitudinal  space-time  correlation  function 
developed  by  Lauchle  (1980)  using  data  fron 
Shivitz  (1975). 


F*  FUNCTION 


This  section  will  define  the  F  function  given  by  Lauchle(1980)  equation  37  which 
is  integrated  over  the  transition  zone.  It  is  assumed  tha  Uc/Uo  is  independent  of 
frequency.  The  spectrum  is  then  calculated  for  single  frequencies  which  represent 
1/3  octave  center  frequencies  from  10  Hz  to  80  KHz. 


T1(KD,Z)  := 


2-astarAx(Z) 


(astarAx(Z)2  +  KD2) 


/  \  (  ( cxstarAx(Z)^ 

T2(KD,Z)  :=  exp(-astarAx(Z)-Z)-sin  KDZ-atan  - 

V  135  ; 

T3(KD,Z)  :=  exp[-astarAx(Z)*(l  -  Z)]  sin  KD(1  -  Z)  -  atanf— 


astarAx(Z) 

KD 


Func  (KD ,  VEL ,  Z)  := 


B(Z)-y(Z)-Z-(T1(KD,Z)-T2(KD,Z)  +  T3(KD,Z)) 


[[4-(B(Z))2-VEL  2]  +  KD2]( astarAx (Z)2  +  KD2) 


F(A ,  C)  :=  A*-  Func  (A ,  C ,  Z)  dZ 

Jo 


FSTAR(A,C,T)  := 


(i.o  +  a2-t2) 


•F(A,C) 
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SPECTRAL  CALCULATIONS 


Define  frequency  range 


ORIGIN  :=  1 


fmax  :=  5000- Hz  KL  :=  1000  k  :=  1 ..  KL  :=  fmax- 10 


(  k 

3-  — - 
KL 


points  :=  1000 
nn  :=  1 ..  points  +  1 


Znn  •“ 


nn  -  1 
points 


I  Ax  Uc 

in*k,nn  :=  Punc  2*7C*f —  ,2.m 

Uc  Uo 


dZ  := 


points 


ppk  •=  Qntk,nn+1 +»ntk,nn)  ^ 


nn  =  1 


Definition  of  function  representing  power  spectral  density: 


Fstark  := 


2-rc-ffc- 


,  f  c  Ax  )  .  2 

1  +  2-Jt-ffc - -rise 

Uc 


Gk  := 


[( W-p2- ASstar2) -Uo-Uc^-Fstark 


8-r2 -it2 


Calculation  of  power  spectral  density: 


SPECGk  :=  101o 


(Pref) 
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dB  //  ( lor  20)  uPaA2/Hz  at  r  m 


RADIATED  NOISE  DUE  TO  FLUCTUATING  DISPLACEMENT  THICKNESS 
FROM  BOUNDARY  LAYER  TRANSITION  USING  LIEPMANN  ANALOGY 

GLPROG81.MCD 


Ax  =  0.028  m 
A5star  =  0.00021m 

Uo  =  21msec"1 

Uc  =  14.7  msec"1 

Ret  =  1500000 


W  =  0.305  m  c=  1525  msec"1 

r  =  lm  v  =  0.00000101m2  sec  1 

rise  =  0.15  p  =  999  kg  m’3 

Uc  _ 

Pref  =  0.000001  kg  m  1  sec  ~ 


SPECTRAL  PLOTS  y(Z)  =  yFarabee(Z)  B(Z)  =  BFarabee(Z) 

120 
110 
100 
90 
80 
70 
60 
50 
40 
30 
20 

1  10  100  1  103  1  104 

Frequency,  Hz 

-  rise  =  1.00 


Power  Spectral  Density 
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MODIFICATION  HISTORY 

Original  7  June  1998  Conversion  of  elements  of  GLPROG  VAX 

program 

Mod  1  7  June  1998  Addition  of  choices  for  intermittency  and  burst  rate 

functions 

Mod  2  4  July  2000  Substituted  trapezoidal  integration 

Spectrum^  :=  f-sec 

Spectrum^  :=  SPECG 

a 

Z:\-.\81WATERFF21.TXT 

Spectrum 
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A .3  Fluctuating  Boundary  Layer  Displacement  Thickness  Using  a  Liepmann 

Analogy  -  Theoretical  Space-Time  Correlation  Function  with  Krane  Dipole 
Source  Model 


Mathcad  file:  KMTRANS00.MCD 
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RADIATED  NOISE  DUE  TO  FLUCTUATING  DISPLACEMENT  THICKNESS 
FROM  BOUNDARY  LAYER  TRANSITION  USING  A  LIEPMANN  ANALOGY 
WITH  KRANE  DIPOLE  MODEL 


[KMTRANSOO.MCD,  Mathcad  2000  Professional] 

This  program  is  the  implementation  of  a  boundary  layer  transition  noise  model 
described  in: 

Lauchle,  G.  C„  1981,  ’Transition  noise  -  The  role  of  fluctuating  displacement 
thickness,’/.  Acous.  Soc.  ofAmer.,  69,  pp.665-671. 

However,  the  source  term  has  benn  modified  to  replace  the  monopole  with  a 
dipole  derived  from  the  measurements  cited  in: 

Krane,  M.  H.,  1992,  "Estimation  of  the  Direct  Acoustic  Radiation  from  Turbulent 
Spots,  Using  Boundary  Layer  Velocity  Measurements,"  Ph.D.  Thesis  in 
Aerospace  Engineering,  The  Pennsylvania  State  University. 


ASSUMPTIONS 

-  builds  on  fluctuating  shear  stress  model  [Lauchle  1981] 

-  zero  pressure  gradient  flow  and  natural  transition  flow  structures,  i.e.  no 
separation). 

-  flow  is  over  an  infinite,  planar,  rigid  surface 

-  based  on  very  low  Mach  number  -  use  only  dipole  components  due  to 
fluctuating  displacement  thickness 

-  flow  is  statistically  homogeneous  in  spanwise  direction,  X3,  non-homogeneous 
in  streamwise  direction,  xi,  stationary  but  Poisson  distributed  in  time  so  we  can 
use  a  classical  correlation  function  for  a  random  sequence  of  bursts,  exp(-2N|t|) 

-  neglect  quadruples  and  octupoles  due  to  fluctuating  Reynolds  stresses  due  to 
lower  radiation  efficiencies  -  however,  their  amplitudes  relative  to  shear  stresses 
may  overcome  this  lower  efficiency. 

-  the  space-time  correlation  function  for  the  indicator  function  <H’>is  direction 
separable 

-  the  lateral  term  integration  derived  from  this  separation  is  equal  to  2Lj(Hi)  = 

2(1.5  Hi  tan  a)  which  is  about  Hi/4  for  a=9.6  degrees  (angle  of  spreading  for 
turbulent  spot) 

-  Uc/Uo  is  independent  of  frequency 

-  the  normalized  rise  time  is  described  as  a  function  of  t|i/Ax 
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DEFINITIONS 


Ax  is  the  length  of  the  transition  zone  in  feet  or  meters 

AS*  is  is  the  change  of  the  boundary  layer  displacement  thickness  from  laminar  to 
turbulent  states  in  feet  or  meters  (AS*  =  ^^turbuient  -  ^*iaminar)  (ASstar  in 
worksheet) 

Uo  is  the  freestream  velocity  in  feet/sec  or  meters/sec 

Uc  is  the  convection  velocity  in  feet/sec  or  meters/sec  (nominally  Uc  =  0.7*Uo) 
a*  is  is  the  flare  parameter  accounts  for  moving  axis  decorrelation 
c  is  the  speed  of  sound  in  feet/sec  or  meters/sec 
Ti  is  the  time  it  takes  S*laminarto  transition  to  dS^^  in  sec 

rise  is  the  non-dimensional  rise  time  for  the  AS*  (rise  =  Uc*Ti/Ax) 

Ret  is  the  transition  Reynolds  number  based  on  Ax 

Rexo  is  the  Reynolds  number  based  on  the  distance  to  the  start  of  the  transition  zone 
W  is  the  spanwise  width  of  the  transition  zone  in  feet  or  meters 

v  is  the  kinematic  viscosity  in  feetVsec 

p  is  the  fluid  density  in  slugs/feet3 

r  is  the  far-field  observation  distance  in  feet  or  meters 

0  is  the  angle  from  the  plane  of  the  wall  to  the  observation  point 
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SPECIFY  VALUES  FOR  FLOW  AND  GEOMETRY  VARIABLES 


l.0\ 

1.25 

1.6 

2.0 

2.5 

3.15 

4.0 

5.0 

6.3 


[Note:  Mathcad  is  operating  in  SI  units  mode  for 
all  calculations.  The  user  may  enter  values  for 
variables  in  any  appropriate  measurement  system 
so  long  as  the  units  are  specified.  The  program 
will  automatically  make  the  conversion  to  SI  units 
during  calculation.] 

Define  1/3  Octave  band  center  frequencies 
a  :=  0..9 

Freqa  :=  Fra-1.0-Hz 


b  :=  10..  19 

Freqb  :=  Fr(b_10)-10.0-Hz 
d  :=  20..  29 

Freqa  :=  Fr(d_2o>- 100.0- Hz 
e  :=  30.. 39 

Freqc  :=  Fr(c_30)- 1000.0- Hz 
f  :=  40..  49 

Freqf  :=  Fr(f^o)- 10000.0- Hz 


Constant- 
s  for  user 
reference. 


Constant  Air  Fresh  Water  Salt  Water 

c  (m/s)  343  1525  1487 

v  (m2/s)  1.51  x  10‘5  1.01  x  10'6  1.01  x  10'6 

p  (kg/m3)  1.21  999  1027 

Pref  OiPa)  20  1  1 


Ax  :=  90-cm 


ASstar  :=  0.278- cm 


A5  :=  1.6-cm 


m 

v  :=  1.51- 105-— 

Uo  :=  21,0 - 

W  :=  .305- m 

sec 

sec 

Uc  :=  0.7-Uo 

r  :=  1.0-m 

Pref  :=  20- 10"  6- Pa 

Ret  :=  1.5- 106 

0  :=  45-deg 

p  :=  1.21-— 

3 

m 
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POWER  SPECTRAL  DENSITY 


This  section  will  define  the  intermittency  function,  y(Z),  and  the  turbulent  spot 
burst  rate,  N(Z),  where  Z  =  r|x/Ax.  To  lessen  computational  burden  later,  the 
integrations  will  use  another  function  B(Z)  =  N(Z)(Ax/Uo).  Obviously,  this  just 
removes  the  factor  (Uo/Ax)  from  the  definition  of  N(Z)  but  this  factor  is  common 
to  both  models  of  interest  and  is  accounted  for  in  the  overall  equation.. 

INTERMITTENCY  FACTOR 

This  is  the  time-averaged  value  of  the  indicator  function  and  represents  the  percentage 
of  time  that  the  flow  is  in  a  turbulent  state.  Either  of  two  models  for  this  factor  can  be 
chosen  by  the  user  for  the  computation  of  the  spectrum.  The  empirical  relation  from 
Farabee,  et  al.  (1974)  was  used  in  Lauchle  (1980). 

yFarabee(Z)  :=  1  -  exp(-4.185-Z2)  Empirical  relations  from  Farabee,  et  al.  (1974). 

yJL(Z)  :=  1  -  exp[-(l  +  3.4  Z)-Z2]  Empirical  relations  from  Josserand  and 

Lauchle  (1986)  which  was  based  on  the  work 
of  Gedney  (1979)  asessing  both  a  Dirac  line 
source  hypothesis  and  a  constant  source 
hypothesis  for  intermitency. 

**CHOOSE  INTERMITTENCY  FACTOR:  y  (Z)  :=  yFarabee(Z) 
NORMALIZED  RISE  TIME 

This  value  has  previously  always  been  held  as  a  constant  throughout  the  transition  zone. 
It  is  often  used  as  a  parametric  variable  in  the  use  of  GLPROG81.  However,  Krane 
measured  a  variation  of  this  time  with  lonitudinal  position  in  the  transition  zone.  The 
user  may  choose  to  use  the  Krane  function  or  just  set  a  constant.  Parametric  values  of 
0.05, 0.15, 0.45,  and  1.0  are  common. 


TK(Z)  :=  0.30  -  «p(-1.45  -  3  Z2)  E"I?irica*  r±?°n  from  flgures  415  an<l  4  20 
rv  '  m  Krane  (1992) 


**CHOOSE  RISE  TIME  FUNCTION: 


Use  constant  for 
higher  frequency 
calculations. 
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TURBULENT  SPOT  BURST  RATE 


This  is  the  time-averaged  number  of  turbulent  bursts  sensed  per  second.  Either  of  two 
models  for  this  factor  can  be  chosen  by  the  user  for  the  computation  of  the  spectrum. 
The  empirical  relation  from  Farabee,  et  al.  (1974)  was  used  in  Lauchle  (1980). 


BFarabee(Z)  :=  1.272Z-exp(-4.185-Z2) 


BJL(Z)  :=  2.38- 


-  Y(Z))-ln 


_L _ 

-Y(Z>) 


♦♦CHOOSE  SPOT  BURST  RATE  FUNCTION:  B  (Z)  :=  BFarabee(Z) 
FLARE  PARAMETER 


This  parameter  is  used  in  the  longitudinal  space-time  correlation  function  to  account 
for  moving  axis  decorrelation  due  to  a  burst  not  remaining  perfectly  correlated  as  the 
distance  increases  beyond  some  characteristic  scale. 

astarAx(Z)  :=  1  +  83.35-Z*  Normalized  flare  parameter  used  in  the 

longitudinal  space-time  correlation  function 
developed  by  Lauchle  (1980)  using  data  fron 
Shivitz  (1975). 
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F*  FUNCTION 


This  section  will  define  the  F  function  given  by  Lauchle(1980)  equation  37  which 
is  integrated  over  the  transition  zone.  It  is  assumed  tha  Uc/Uo  is  independent  of 
frequency.  The  spectrum  is  then  calculated  for  single  frequencies  which  represent 
1/3  octave  center  frequencies  from  10  Hz  to  80  KHz. 


T1(KD,Z) 


2-exstarAx(Z) 


(astarAx(Z)2 


+  KD2) 


2 


T2(KD,Z)  :=  exp(-exstarAx(Z)*z)sin 


f 


KDZ-atan 


V 


W 


exstar Ax(Z) 


T3(KD,Z)  :=  exp[-astarAx(Z)-(l  -  Z)]-sin 


KD-(1  -Z)  -  atan 


exstarAx(Z) 
KD 


M 


Func  (KD ,  VEL ,  Z) 


B(Z)-y(Z)Z(Tl(KD,Z)  -  T2(KD,Z)  +  T3(KD,Z)) 

\_ 

[[4-(B(Z))2-VEL'2]  +  KD2](  exstarAx(Z)2  +  KD2) 2 


F(A,C) 


Func(A,C,Z)dZ 


This  is  the 
integration  that  is 
performed  below 
using  a  trapezoidal 
method. 


FSTAR(A,C,T) 


(l.0  +  A2-T2) 


•F(A,C) 
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SPECTRAL  CALCULATIONS 


Define  frequency  range  ORIGIN  :=  1 

fmax  :=  5000- Hz  KL  :=  1250  k  :=  1..KL 
points  :=  1000 

nn  :=  1 ..  points  +1  ^  nn  -  1 

Am  : - : — 

points 


fk 


3- 

fmax- 10 


1“-  J 


f 


intk.nn  :=  Func 


2-n-fjc*— 

Uc  Uo 


\ 

>^nn 

/ 


dZ  := 


points 


FPk  (intk,nn+l -t- intk.nn)  ^ 

nn  =  1 
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Definition  of  function  representing  power  spectral  density: 


FFSTARk  :=  F 


'Y  e  Ax 
2-n-fk-— 
Uc 


Y» 


f 


1.0  + 


2-tt-fk 


Ax 

Uc 


=  FFk 


T2 


2-71  fk  /  \ 
Argk  :=  0.1 - Axcos(0J 


Dipolek  :=  4- 


sin 


L  L 


2-rc-fk 


-AS-cos(o) 


Y 


From  Kinsler  &  Frey  - 
constant  rise  time  (and 
therefore,  d) 

T50  :=  0.1 


Dipole2k  :=  4 


sm 


\2J 


A52  +  (Ax-T50)2]  -cos(o) 


n2 


DipoleE2k 


;=  2  +  "S17-Arg-  (l  -  2-cos(Argk)) 


Argk 


Gk  := 


(wp2A8star2)-Uo-Uc2-FFSTARk 


8-r2-7t2 


Modified  by  calculating  the 
expected  value  for  a  uniform 
distribution  of  the  source 
separation,  d  which  is  a 
function  of  rise  time  which  is 
itself  a  function  of  Z 


ODk  :=  DipolekGk  0D2k  _  Dipoll.2k.Gk 

GDE2k  :=  (DipoleE2k)  -Gk 


Calculation  of  power  spectral  density: 


SPECGk  10-lod 


Gk 


(Pref) 


Hz 


SPECGDk  :=  10-lod 


GDk 

(Pref)2 


Hz 


SPECGD2k  :=  10-lod 


GD2k 

(Pref)2 
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SPECGDE2k  :=  10- log) 


GDE2k 

(Pref)2 


Hz 


Y(Z)  =  YFarabee(Z)  B(Z)  =  BFarabee(Z) 
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RADIATED  NOISE  DUE  TO  FLUCTUATING  DISPLACEMENT  THICKNESS 
FROM  BOUNDARY  LAYER  TRANSITION  USING  LLEPMANN  ANALOGY 
BUT  DIPOLE  SOURCE  FROM  KRANE 

KMTRANSOO.MCD 


Ax  =  0.9  m 
ASstar  =  0.00278  m 

Uo  =  21msec  1 

Uc  =  14.7  msec'1 
Ret  =  1500000 
SPECTRAL  PLOTS 


W  =  0.305  m  c  =  343  msec  1 

A5  =  0.016  m  v  =  15 1000  m2  sec'1 

r  =  lm  p  =  1.21  kg  m'3 

0  =  45  deg 

Pref  =  0.00002  kg  m'1  sec*2  —  =  0.7 

Uo 

Y(Z)  =  vFarabee(Z)  B(Z)  =  BFarabee(Z) 
T(Z)=0.15  Integration  points  =  1000 


10  100  l  io3  l  104 

Frequency,  Hz 


-  Dipole  d=delta  del* 

-  Dipole,  d=2deltax  T 

-  Dipole-  d  varies  with  uniform  prob 
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MODIFICATION  HISTORY 


Original  8  June  2000 

time 

Mod  1  14  June  2000 


Conversion  of  GLPROG81  to  replace  monopole 
with  dipole  source  and  put  in  normailized  rise 

as  a  function  of  Z 

Upgraded  dipole  source  term  to  eliminate  kd«l 
low  frequency  assumption. 


Spectrum^  :=  f-sec 

Spectrum^  :=  SPECGD2 
Spectrum^  :=  SPECGDE2 

a 

ZUKMTRWB  OAF  F  21  NBW.TXT 


Spectrum 
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A.4 


Fluctuating  Boundary  Layer  Displacement  Thickness  Using  a  Liepmann 
Analogy  -  Empirical  Space-Time  Correlation  Function 


Mathcad  file:  JLTRANS99.MCD 
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BOUNDARY  LAYER  TRANSITION  ZONE  SOUND  RADIATION  USING 
JOSSERAND’S  EMPIRICAL  SPACE-TIME  CORRELATION 


[JLTRANS99.MCD,  Mathcad  2000  Professional] 

This  program  is  the  implementation  ot  a  boundary  layer  transition  noise 
model  described  in: 

Lauchle,  G.  C.,  1981.  Transition  noise  -  The  role  of  fluctuating 
displacement 

thickness,’/  Aeons.  Soc.  ofAmer..  69,  pp.665-671. 

Josserand,  M.  A.  and  Lauchle.  G.  C.,  1989.  "Cross-Spectral  Density  of  the 
Wall  Pressure  Fluctuations  Under  a  Turbulent  Boundary  Layer," 

/  Sound  Vib.  128.  pp.  519-522. 


Josserand,  M.  A.  and  Lauchle.  G.  C..  1990.  "Modeling  the  Wavevector- 
Frequency  Spectrum  of  Boundary  Layer  Wall  Pressure  During 
Transition  on  a  Flat  Plate."  Tran  s.  ASME,  J.  Vib.  and  Acoustics 
112,  pp.  523-534. 


ASSUMPTIONS 

-  builds  on  fluctuating  displacemtn  thickness  model  [Lauchle  1981] 

-  zero  pressure  gradient  flow  and  natural  transition  flow  structures,  i.e.  no 

separation ). 

-  flow  is  over  an  infinite,  planar,  rigid  surface 

-  based  on  very  low  Mach  number 

-  uses  enmpitrical  description  for  Cross  correlation  function  developed  by 

Josserand 

-  Uc/Uo  is  independent  of  frequency 


DEFINITIONS 


Define  transition  zone  length  (Ax)  and  free-stream  velocity  (U) 

Ax  :=  0.9*m  U  :=  21- 
sec 

Turbulent  spot  convection  velocity  Uc  :=  0.8- U 


Intermittancy  factor.  Josserand  &  Lauchle  (1990)  eqn.  1 1 
c3:=  1.0  c4  :=  3.4 


yM  :=  1  -  exp) 


-fd  +  c4JE\f’E')2 

l  a*A“IJ 
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r(tix,Sx,5z,t)  :=  yu(t\x,^)-[yc1(tix,4x)  +  (l  -  Yd(rix,4x))-efl(^t,T)-ef2(5x’t)-ef3<^‘,5z’t)]-ef4(5x,T> 
Discretize  the  separation  distance  range 
Az(rix)  :=  3-iix  tan(9.4  deg) 

IL  :=  31  i  :=  O..IL 

JL:=15  j:=0..JL  dx,,:=Ax/— ' )  di  s  :=  -  - 

J  I^ilJ  JL  2 

4x(tix)  :=  dx  -  t]x 
^z(tix)  :=  Az(tix)-dz 


236 


Discretize  the  transition  zone  range  and  time  delay 
KL  :=  31  N  :=  1000 

k :=  0..KL  n  :*  0..N-  1 


T  :=  12-— 
U 


^  :=  Ax - 

KL 


tn:=  2---1  - 

n  ,  : 


ffk.n  :=  ^x(tlXk) ,  ^(tixj ,  x  j 


Integrate  over  the  separation  range 
Ax-Az^xJ  tl1  %!' 

2  2  [(*«..), j + *  (*<•■),.», + <M. ,»,] 

i  =0  j  =0 


Integrate  over  the  transition  zone 
Ax  KLrl 

n*"  “TIE'S  (intk.»+intk+'.n) 

k  =0 


237 


Discrete  Fourier  transform  of  spatially  integrated  correlation  function 


G:=  — -CFFT(RR) 
fs 


k:=  1..0.5-N 


Farfield  radiated  sound,  eqn.  21 

p  :=  1.21 ~  W  :=  0.305*  m  ASstar  :=  0.278  cm  ti  :=  0.2  — 

m3  Uc 

pref  :=  2010'6  Pa 


4>v:= 


p2*A8star2*(2-n-fk)4 
4-n2-[l  +  (2-«-fk-ti)2] 


dBk:= 


/ 

10-log 

V 


»k-Hz 

_ -f  2 

prct 


\ 

/ 


Lauchle  Model  w/Josserand  Correlation 


r:=  l-m 
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Spectrum^  :=  f-sec 
Spectrum^  :=  dB 


B 

ZIU.TRANSW211.TXT 


Spectrum 


AJ 


Noise  Generated  by  Boundary  Layer  Transition  on  Symmetric  Hydrofoils  • 
Howe  Formulation 


Mathcad  file:  MSHTRANS99.MCD 
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BOUNDARY  LAYER  TRANSITION  NOISE 
FROM  A  SYMMETRIC  HYDROFOIL 


[MSHTRANS99.MCD,  Mathcad  2000  Professional] 

This  program  is  the  implementation  of  a  boundary  layer  transition  noise  model 
described  in: 

Howe,  M.  S.,  1999,  ’Noise  Generated  at  Low  Mach  Numbers  by  Boundary  Layer 
Transition  on  Symmetric  Hydrofoils  and  Axisymmetric  Headforms,’  Applied 

Research  Laboratory  /  Penn  State  University  Technical  Memorandum,  99-017 
February  1999. 

ASSUMPTIONS 

-  very  low  Mach  number.  2D  flow; 

symmetric,  finite  thickness  plane  (better  termed  a  hydrofoil); 

-  hydrofoil  chord  is  acoustically  non-compact; 

•  hydrofoil  span  is  much  greater  than  the  boundary  layer  thickness; 

-  hydrofoil  thickness  is  acoustically  compact,  and  that  the  leading  edge  is 

within  the  nearfield  of  the  transition  zone; 

-  zero  angle  of  attack  and  boundary  layer  fluctuations  on  opposite  sides  of 

the  foil  are  equal  and  independent  so  an  estimate  of  one  side  can 
just  be  increased  by  3  dB; 

-  the  integration  of  the  vorticity  is  restricted  to  the  non-linear  region  of  the 

boundary  layer  above  the  viscous  sublayer,  so  bound  vorticity  on  the 
hydrofoil  is  ignored;  and 

;  h,§h  Reynolds  number  so  the  contribution  of  tangential  shear  stresses  is 
ignored. 


SPECIFY  VALUES  FOR  FLOW  AND  GEOMETRY  VARIABLES 


ORIGIN*  | 


Physical  constants  for  water  at  10°C 

paiooo-ii  V  a  1.284- JO'6  — 

m3  sec 

Define  hydrofoil  parameters 

a  s  0.0013  p  a  0.02  h  3  0.856  m 


C  *  1447  —  Prof  :=  110"6  Pa 
sec 


Ctrans  *  1  Ntrans  *  151 
241 


Define  free-stream  velocity,  calculate  Mach  number  and  Reynolds  number 

Uh 


U«  15-2- 

sec 


S-OOl 

c 


=  1  x  10' 


Define  convection  velocity  using  assumption  stated  on  page  12  Uc  :=  0.7-U 
Define  frequency  range 

fmax  a  2000  Hz  KL  a  200  k  :=  1 ..  KL  fk  :=  fmax- 10 1  ') 
Specify  frequency-dependent  window  length  for  integration 


Uc 


v0.7 


(Ok  :=  2-refk 


OnaXjj  :=  2-  — —  I  +  ^trails  Nwindow  :=  151 

v*h) 

Generate  hydrofoil  profile  up  to  transition  point  via  conformal  mapping  (eqns.  2.11  and  2.12) 
n  :=  1 ..  Ntrans 


Con  :=  Orans - - 

Ntrans  -  1 

fc(0  :=  — -  1  -  ln($  +  yjc,2  -  l))  +  — 

*  2 

fcfoO  + 

:=  .  a  X0» M20")  yOn  -  WzOn) 
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Calculate  arc  length  up  to  transition  point 


n  :=  1 ..  (Ntrans  -  1 ) 

dso,  :=  ^(x on+l  -  xon)2  +  (y0n+1  -  yon)2 

a 

Calculate  arc  length  for  points  beyond  transition  point  (different  grange  for  each  frequency) 
n  :=  1 ..  Nwindow 


Cn.k  •“  (Cn^k  ~  Ctrans) - - - +  £tran$ 

Nwindow  -  i 


Vk 


fc(Co.k)  +  «-fc 

h - - - 1 


P  ; 


1  +  a 


*«.k  :=  Refa.k)  yn.k  :=  Im(*n.k) 


nn  :=  1 ..  (Nwindow  -  1) 


dSna-k  :=  >/(* on+l.k  ~  *nn.k)2  +  (ynn+l.k  ~  Ynn.k)2 
sl,k:=  SO 


an 

snwl,k  :=  SO  +  £  dsjik 
j  =1 
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Plot  profile  up  to  transition  point  (red)  and  windowed  integration  range  (blue) 
for  a  specified  frequency 


kk:=50  — =  5.761 

Uc 


nl  :=  1 ..  Ntrans  n2  :=  1 ..  Nwindow 


SPECTRAL  CALCULATIONS 


Calculate  transition  zone  intermittency  function  (eqns  3.14  and  3.15) 
f  U-so  V0'2 


A  :=  17so 


V  v  ) 


cl  :=  4.18  O'* 


Y.,k  :■  1  -  50 j 


Calculate  boundary  layer  displacement  thickness  and  friction  velocity  (eqn.  3.17  and  page  14) 

Vuv 

V 


Sstarn>k:=!  0.046-h-f  — 
h 


V 


V 


vstar  :=  0.035- U 


Calculate  point  pressure  spectrum,  based  on  Chase  data  (eqn.  3.16) 
ctp  :=  0.12 


^PPn.k :»  (p 


&tarn.k  (  &urn.kV 

vstar  )  --.Yo>k - 2 - 1 - “ - L_ 

u  f/-  .  s->  Tl-5 


'  &tarn.kV  : 

»k - —\  +  ctp 


r  f  a)th^l  ®k 

Calculate  integrand  of  eqn.  3.19  im,  k  :=  +  d  Uc  ij-e" 


nn  :=  l ..  (Nwindow  -  1 ) 
^Cmi.k  :=  ?na+l.k  ~  Cnn.k 


Apply  window  function  (n  C„  k  -  Grans  V 

wn,k  :=  cod  — ■  *  1 


2  £maxk  -  £trans 


*ntwn,k  :=  wn.k-*nt„,k 
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Plot  the  integrand  for  a  specified  frequency 


kJc  :=  50 


(0]£*h 

Uc 


=  5.761 


Evaluate  windowed  integral  using  trapezoidal  integration 
nn  :=  1 ..  (Nwindow  -  1 ) 

c  V't  (intwnn+l.k  +  >ntwnn.k) 

Fk’2v - 2 - d?ra,-k 


nn 
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Integral  (magnitude  squared)  for  low  frequencies,  eqn  3.20. 


F2.owk:=0,4^V06(p.vstarfA\f^a 


V  v  ) 


1” 


i  +  iL 

v  32, 


Uc 


Integral  (magnitude  squared)  for  high  frequencies,  eqn  3.21. 
^2,k  ~  Cl.k 


d^dsov  := 

^.k  ~  sl,k 

F2hlgh|C  :=  cl-(p-vstar')  /“l 


Sound  radiation  from  eqn.  3.19 


XX  Uc 

Mc:= —  L  :=  10- h 
c 


'h*  * 

— d^dsok 

V*  J 


_U 

Uc 


^oyhY5 

l  Uc  J 


r:=  l  m  y  :=  90-deg  9  := 


*  OJSMcfhlA 

®‘s — ; — r 

*  r ) 


(  a  N 
1  +  — 


v  i  +  ay 


sin(v) 


•sir|fj  (lFk|)2 


<t>lowk  := 


0.35 -Me  ( h-L^) 


f  a  ^ 

1  +  — 


1  +a 


)s# 


sin(v)-sinj  —  I  F21owk 


‘fthighfc  := 


0.35-Mc  ( h-L') 


f  a  >k 

1  +  — 


v  1  +  a/ 


•sin(v)-sin^-jj  F2highk 


90- deg 
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Calculate  and  plot  dimensional  spectrum 


°k  := 


10- log 


(Pref)2 


+  11 


f 


k;: 


(Oy_ 

2-71 


The  1 1  dB  accounts  for  conversion 
to  G(f),  the  one  sided  spectrum. 


Calculate  and  plot  Chase  spectrum  for  TBL  noise 


CT  :=  0.0046 


b  :=  0.75 


hbar:=  3 


Gchase*  :=  10- loo 


‘WZhasCk-Hz 

(Pref)2 


+  11 


Calculate  combined  spectrum  for  transition  noise  and  TBL 


Gtotalk  :=  10- log 


4>v 


+  <t>Chase)JJ-Hz 
(Pref)2 


+  11 
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dB  //  t  uPi»A2/Hz  9rm 


Transition  Noise 


*“•  Quadrupole  Noise 
•••••  Total  Noise 


MODIFICATION  HISTORY 


Original 
Mod  1 
Mod  2 


2  June  2000 

5  June  2000 

6  June  2000 


Full  implementation 

Addition  of  windowing  function  for  integration 
Addition  of  dimensional  spectra  and  quadrupole 
spectra 
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dB //  1  uP#A2/Hz  #rm 


BL  Transition  Noise  For  Symmetric  Foil 


— *  Quadmpoie  Noise 
Total  Noise 


MODIFICATION  HISTORY 


Original 
Mod  1 
Mod  2 


2  June  2000 

5  June  2000 

6  June  2000 


Full  implementation 

Addition  of  windowing  function  for  integration 
Addition  of  dimensional  spectra  and  quadrupole 
spectra  and  Josserand  intermittency  function 
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A.6  Calculation  of  Boundary  Layer  Descriptors  -  AIR 

Mathcad  file:  TRANSITION  AIR.MCD 


253 


TRANSITIONAL  FLOW  HYDRODYNAMIC  DESCRIPTORS 
-AIRFLOW 


[TRANSITION  AIR.MCD,  Mathcad  2000  Professional] 

This  worksheet  is  to  be  used  to  calculate  the  transition  zone  length,  transition  zone 
Reynolds  numbers,  convection  velocity,  and  displacement  thickness  peak  amplitude. 

From  Josserand,  I  am  assuming  a  initial  transition  at  a  Reynolds  number  of  8  x  10A5. 
The  transition  Reynolds  number  (where  intermittancy  is  50  percent)  is  equal  to  1.5  x 
10A6.  The  Reynolds  number  at  the  end  of  the  transition  zone  is  an  average  2.2  x  10A6. 
These  are  based  on  experiments  in  air  for  zero  pressure  gradient  flat  plate  flows. 


^9^ 


Uo  := 


12 

15 

18 

21 

24 


m 

sec 


DEFINITIONS 


Uc  :=  Uo-0.7 


Convection  velocity  ratio 


x  :=  0.90  m 


V 


:=  1.5- 10' 


5  m_ 

sec 


j  :=  0..6 


Measurement  position 

Kinematic  viscosity, 
air  =  1.50  x  10'5  m2/sec 
water  =  1.14  x  10-6  n^/sec 


Rexj  := 


8xlamj  :=  5.0-x-(Rexj)  2 


Reynolds  number  based  on 
measurement  position  for  a 
flat  plate 

Laminar  99%  boundary  layer 
thickness 
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-1 

Sxturbj  :=  0.37 -x- (Rex  j)  5 

Turbulent  boundary  layer 
thickness 

800000 

XO;  := - V 

J  Uoj 

Effective  leading  edge  point  of  transition 
zone 

XO; 

RexOj  :=  Uo; - 

V 

Transition  LE  Reynolds  number 
(by  definition) 

Length  of  transition  zone 

Chen  and  Thyson  use  60, 

Josserand  found  40 

«  xoi 

Sstarlam;  :=  1.743 - - — 

(RexOj)0-5 

Laminar  displacement  thickness 
at  LE  of  transition  zone 

2200000 

xt,  := - v 

J  Uoj 

Effective  distance  at  beginning 
of  fully  turbulent  region 

e.  Xti 

Sstarturb:  :=  0.046- - - - 

J  \ 

(2200000) 5 

Displacement  thickness  at 
beginnning  of  fully  turbulent 
region 

ASstarj  :=  Sstarturbj  -  Sstarlamj 

AS*  calculation 

6 

Refistarturbj  :=  0.018-(Rcxj)7 
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CALCULATED  TRANSITIONAL  FLOW  VARIABLES 


Uo  = 


(  9^ 

12 

15 

18 

21 

24 

v27y 


m 

sec 


5x1am  = 


xo  = 


Uc  = 


f  6.3  ^ 

8.4 

10.5 

12.6 

14.7 

16.8 
18.9 


m 

sec 


—  Rex  = 


,5  \ 


5.4  x  Iff 


7.2  x  ur 
9x  105 
1.08  x  106 
1.26  x  106 
1.44  x  106 
1.62  x  106 


Ax  = 


(  86.569 
64.927 
51.941 
43.284 
37.101 
32.463 
28.856  j 


cm 


'0.612^ 

'  2.377  > 

r0.909> 

f  0.26  "l 

0.53 

2.244 

0.682 

0.195 

0.474 

2.146 

0.545 

0.156 

0.433 

cm  Sxturb  = 

2.069 

cm  Sstarturb  = 

0.454 

cm  Sstarlam  = 

0.13 

0.401 

2.006 

0.39 

0.111 

0.375 

1.953 

0.341 

0.097 

^ 0.354 j 

1.908  j 

^0.303; 

^0.087y 

cm 


f  133.333  > 

1.475  x  103 

'0.649  > 

100 

1.887  x  103 

0.487 

80 

2.285  x  103 

0.389 

66.667 

cm  Re&tarturb  = 

2.672  x  103 

Afistar  = 

0.325 

57.143 

0.278 

50 

3.049  x  103 

0.243 

44.444  j 

3.419  x  103 

^  3.782  x  103y 

cm 


xt  -  xo  = 


^233.333  ^ 

'  86.569  N 

175 

64.927 

140 

51.941 

116.667 

cm  Ax  = 

43.284 

100 

37.101 

87.5 

32.463 

^  77.778  j 

^  28.856  , 

cm 
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Non-dimensional  boundary  layer  axial  velocity  calculation  for  a 
laminar  boundary  layer. 
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Calculation  of  Boundary  Layer  Descriptors  ■  WATER 


Mathcad  file:  TRANSITION  WATER.MCD 


Sxturbj  :=  0.37  x- (Rex j)  5 

Turbulent  boundary  layer 
thickness 

800000 

XO;  := - V 

J  Uoj 

Effective  leading  edge  point  of  transition 
zone 

XO; 

Rexoj  :=  Uoj — - 

Transition  LE  Reynolds  number 
(by  definition) 

luv 

Length  of  transition  zone 

Chen  and  Thyson  use  60, 

Josserand  found  40 

XOi 

Sstarlam;  :=  1.743 - - — 

(Rexoj)05 

Laminar  displacement  thickness 
at  LE  of  transition  zone 

2200000 

xti  := - v 

J  Uoj 

Effective  distance  at  beginning 
of  fully  turbulent  region 

Xtj 

Sstarturb;  :=  0.046 - - - 

J 

(2200000) 5 

Displacement  thickness  at 
beginnning  of  fully  turbulent 
region 

ASstaij  :=  Sstarturbj  -  Sstarlamj 

AS*  calculation 

6 

ReSstarturbj  :=  0.018*(Rcxj)7 
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Non-dimensional  boundary  layer  axial  velocity  calculation  for  a 
laminar  boundary  layer. 


'  0.0  > 

^0.0> 

0.1 

0.21 

0.2 

042 

0.3 

0.64 

0.4 

0.86 

0.5 

1.10 

unondim  := 

0.6 

n  “ 

1.33 

0.7 

1.61 

0.8 

1.94 

0.9 

2.4 

0.95 

2.79 

^0.9999^ 

,6.0; 

h  :=  0..  11 


( 


zh :=  nh 


2-v* 


V 


UOj 


For  Uo;  =  15  msec 


-l 


unondim  = 


0  ^ 

f  0 

0.1 

7.767  x  10" 

0.2 

1.553 

0.3 

0.024 

0.4 

0.032 

0.5 

0.041 

0.6 

Z  = 

0.049 

0.7 

0.06 

0.8 

0.072 

0.9 

0.089 

0.95 

0.103 

1  ; 

^  0.222 

^-3 


cm 
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QAF  NOZZLE  CONTOUR  DESIGN 


[NOZZLE  CONTOUR.MCD,  Mathcad  2000  Professional] 

Program  to  plot  shape  of  a  nozzle  contour.  It  uses  a  seventh  order  polynomial.  The 
coefficients  were  found  using  a  Mathematica  program  by  setting  the  following  boundary 
conditions: 

-  the  inlet  and  outlet  dimensions  were  set  at  semi-diameter  dimensions  of 

Inlet=  29  inches  and  outlet=6  and  3  inchesdepending  on  side 

-  the  first  and  second  derivatives  at  x=0  and  L  were  set  to  zero 

-  the  first  derivative  at  x=L/2  was  set  to  -1.0 

-  the  coefficient  for  x3  was  set  to  zero  in  order  to  get  a  determinate  set  of 

equations 

-  L  was  chosen  to  be  inlet  diameter  times  0.80 


Coefficients  for  top  and  bottom  Coefficients  for  sides  of 
nozzle  contour  nozzle  contour 


at:=  -0.000969861 
bt:=  0.00476745 
ct  :=  0.0166601 
dt  :=  -0.089322 
et:=  0.000 
ft  :=  0.000 

gt  :=  0.000 
ht  :=  2.41667 


as  :=  0.00135161 
bs  :=  -0.0214136 
cs  :=  0.113844 
ds  :=  -0.206753 
es  :=  0.000 
fs:=  0.000 
gs  :=  0.000 
hs  :=  2.41667 


ytop{x)  :=  at-x7  +  btx6  +  ct-x5  +  dt-x4  +  etx3  +  ft-x2  +  gt-x  +  ht 

yside(x)  :=  as-x7  +  bs-x6  +  cs-x5  +  ds-x4  +  es-x3  +  fs-x2  +  gs-x  +  hs 

dimensions  are  in  feet  L=3.867  ft 

m:=  0..  186 

yclm  :=  O.OOdistn, 

ySm  :=  ysidc{distm) 


distg,  :=  m-0.0208333 
ytm  :=  ytof^dis^) 
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Y  height  above  enterline  in  feet 


Contour  for  side  walls  of  nozzle 
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m:=  0..33 


Longitudinal  axis 
point  in  feet  from 
LE  of  Nozzle  (X) 

distjn  = 


0 

0.0208 

0.0417 

0.0625 

0.0833 

0.1042 

0.125 

0.1458 

0.1667 

0.1875 

0.2083 

0.2292 

0.25 

0.2708 

0.2917 

0.3125 

0.3333 

0.3542 

0.375 

0.3958 

0.4167 

0.4375 

0.4583 

0.4792 

0.5 

0.5208 

0.5417 

0.5625 

0.5833 

0.6042 

0.625 

0.6458 

0.6667 

0.6875 

Distance  from  Nozzle  Distance  from  Nozzle 
centerline  in  feet  for  centerline  in  feet  for 
top  and  bottom  both  sides 


ytm  = 


2.4167 

2.4167 

2.4167 

2.4167 

2.4167 

2.4167 

2.4166 

2.4166 

2.4166 

2.4166 

2.4165 

2.4164 

2.4163 

2.4162 

2.4161 

2.4159 

2.4156 

2.4154 

2.415 

2.4147 

2.4142 

2.4137 

2.4131 

2.4124 

2.4117 

2.4108 

2.4099 

2.4088 

2.4076 

2.4063 

2.4049 

2.4033 

2.4016 

2.3997 

ySm  = 


2.4167 

2.4167 

2.4167 

2.4167 

2.4167 

2.4166 

2.4166 

2.4166 

2.4165 

2.4164 

2.4163 

2.4162 

2.416 

2.4157 

2.4154 

2.415 

2.4146 

2.414 

2.4134 

2.4126 

2.4118 

2.4108 

2.4097 

2.4084 

2.407 

2.4054 

2.4037 

2.4017 

2.3996 

2.3973 

2.3948 

2.392 

2.389 

2.3858 

Example  of 
coordinates 
passed  to  CAD 
system  for 
manufacture. 
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A.9  Calculation  of  Measurement  Errors  from  a  Two  Microphone  Intensity  Probe 
Mathcad  file:  INTEN_ERROR.MCD 
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MEASUREMENT  ERRORS  FROM  A  TWO  MICROPHONE 
INTENSITY  PROBE 


[INTEN_ERROR.MCD,  Mathcad  2000  Professional] 

This  workspace  is  for  calculation  of  measurement  errors  which 
are  associated  with  two  microphone  P-P  intensity  probes. 


Co  :=  343--— 

Speed  of  sound ,  m/sec 

sec 

f  :=  5000  Hz 

Acoustic  frequency,  Hz 

k  :=  2-n-— 

Co 

Acoustic  wavenumber,  1/m 

s  :=  .0254- m 

Microphone  separation  distance, 

Finite  difference  approximation  error.  This  is  a  high 
frequency  limit. 


FDAE  =  -5.047 


FDAE  :=  10.0-log^  -°(k'S)j 


To  stay  within  a  1  dB  error  limit,  the  upper  frequency  is: 

2470  Hz 

An  alternative  upper  threshold  is  given  by  the  -3  dB  point: 

_  , .  Co 

Fu  :=  1.9 - 

2-Tt-s  Fu  =  4083.515  sec*1  Hz 

A  low  frequency  threshold  can  be  developed  for  a  source  such 
as  the  TBL  which  has  a  correlation  length  associated  with  it. 

In  Lauchle’s  INCE  paper,  he  provides  an  equation: 

u  :=  21  m/sec 

Fl  =  795.768  m'1  Iz 


for  flow  velocity, 


0.7- 


dB 
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Thumbrule  for  maximum  inaccuracy  of  +/- 1.0  dB  is 
0.15<ks<1.16. 

Co 


fl  :=  0.15- 


f2  :=  1.16 


2-n-s 

Co 


fl  =  322.383  msec 


i 


2  ns  f2  =  2493.093  msec'1 

Bias  errors  for  different  source  fields,  Fahy,  102-105 

i  :=  0..  1250  freq.  :=  4-i-Hz 

Point  Monopole  field 
r  :=  5-cm  h  :=  1.27-cm 

EIMj  :=  10log(l  +  elMj) 

Point  dipole  field  intensity  bias  error 

/  c—~  _  \2 


elD-,  := 


-2 


EID;  :=  101og(l  +  eH>,) 


Point  Dipole  Field 
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APPENDIX  B 


QUIET  AIRFLOW  FACILITY  DEVELOPMENT 

This  appendix  provides  some  historical  perspective  of  development  efforts  for  the 
Quiet  Airflow  Facility.  The  facility  has  progressed  through  many  stages  during  the 
several  years  covered  by  this  investigation.  All  of  the  efforts  have  been  focused  on 
reducing  the  background  noise  levels.  Section  5.1  describes  the  current  facility  while  this 
appendix  describes  interim  configurations  and  data  acquisition  tools. 

This  facility  is  basically  a  quieted  air  source  connected  to  a  flow  management 
nozzle.  The  original,  symmetric  nozzle  was  fitted  with  an  internal  contoured  ramp  to 
create  a  semi-circular  orifice  rather  than  the  circular  one  that  previously  existed.  This 
nozzle  formed  an  open  jet  with  an  opening  axis  perpendicular  to  a  planar  baffle.  A  plate 
mounted  in  the  semi-circular  orifice  formed  a  wall  jet  over  the  test  plate.  This 
configuration  is  shown  in  Figure  B.l. 

B.l  Power  Source  Improvements 

The  Buffalo  Forge  centrifugal  blower  was  originally  powered  by  a  20  HP  AC 
motor  driven  by  a  Borg-Wamer  inverter.  This  has  been  replaced  by  an  ABB  inverter . 
This  was  necessary  because  of  the  lack  of  spare  parts  for  the  Borg-Wamer  system.  This 
old  inverter  had  been  used  a  long  as  possible  because  it  produced  a  “sloppy”  excitation 
that  was  a  much  smoother  sinusoid  than  the  modem,  six-step,  digitally  generated 


271 


sinusoidal  excitation.  This  tended  to  provide  a  less  tonal  rich  background  noise  level  than 
the  more  modern  inverters.  A  series  of  measurements  were  made  to  characterize  the 
impact  of  this  change.  One  example  of  the  reduction  in  the  tonal  levels  due  to 
replacement  of  the  Borg-Wamer  inverter  with  the  ABB  inverter  is  shown  in  Figure  B.2. 

In  some  cases,  the  levels  decreased  by  30  dB  in  the  frequency  range  of  interest. 

B.2  Blower  Noise  Improvements 

The  initial  blower  isolation  box  did  not  have  a  muffler  on  the  inlet.  A  sheet  metal 
transition  section  was  added  from  the  plenum  box  intake  hole  to  the  blower  suction  port. 
A  rubber  bellows  attaches  the  blower  discharge  to  the  attenuation  plenum.  A  Conic-Flow 
Silencer  muffler  has  been  attached  to  the  inlet  to  the  blower.  A  throttle  can  be  mounted 
on  the  muffler  intake  to  provide  the  coarse  adjustment  of  the  suction  strength.  It  is  not 
known  if  the  presence  of  the  branches  off  the  inlet  tube  between  the  muffler  and  the 
blower  intake  causes  any  three  dimensional  distortions  to  be  ingested  into  the  blower. 
(These  branches  are  for  the  suction  manifold  which  has  been  used  in  previous  tests.) 
However,  some  of  the  low  frequency  rumble  may  be  attributable  to  this  but  the  turning  at 
the  elbows  and  probable  separated  flow  is  likely  a  more  important  source  of  noise.  The 
air  is  discharged  from  the  source  plenum  vertically  through  a  40  cm  diameter  flexible 
vinyl  hose  and  turned  with  a  damped  sheet-metal  elbow  into  another  piece  of  flexible 
hose  for  entry  in  the  flow  management  section  and  nozzle  of  the  wall  jet  facility.  This 
still  didn’t  prove  to  be  enough  quieting,  so  an  additional  plenum  was  added  before  the 
flow  management  section  in  the  final  configuration. 


272 


B3  Noise  Reduction  at  the  Test  Plate 


The  entire  flow  management  assembly  was  mounted  on  pneumatic  tires  for 
vibration  isolation  and  ease  of  movement  into  and  out  of  the  chamber  wall  opening.  In  the 
final  configuration,  the  wheels  were  not  used  because  the  facility  did  not  have  to  be 
moved  when  located  on  the  floor  of  the  anechoic  chamber.  However,  some  rubber 
isolation  material  is  used  between  the  floor  and  the  flow  management  assembly. 

The  original  configuration  for  the  test  plate  called  for  a  false  bottom  to  be 
installed  in  the  existing  nozzle  so  that  the  discharge  orifice  could  be  changed  from 
circular  to  semi-circular.  The  plate  was  to  be  mounted  such  that  it  extended  the  false 
wall.  This  is  shown  in  Figure  B.3. 

In  order  to  provide  for  realistic  boundary  layer  growth  over  the  flat  plate,  as  well 
as  to  prevent  the  formation  of  a  free  shear  layer  at  the  outer  edges  of  the  jet  orifice,  an 
acoustically  transparent,  but  flow  impermeable  semi-circular  mylar  tube  was  placed  over 
the  test  plate  (Figure  B.4).  This  was  engineered  by  Walter  A.  Kargus  (1997).  It  is 
attached  to  aluminum  semi-circular  supports  at  each  end  and  is  lightly  restrained  by  a  few 
thin  diameter  wire  hoops  in  between.  During  early  testing  of  the  acoustically  transparent 
tunnel,  it  collapsed  due  to  its  lack  of  inherent  structural  rigidity  when  subjected  to  the  low 
pressure  region  created  from  the  flow.  In  an  attempt  to  prevent  this,  cheese-cloth  was 
added  across  the  tunnel  termination  as  a  resistance  to  create  a  positive  gage  pressure 
within  the  tunnel.  This  was  aerodynamically  effective  but  noisy.  A  wind  sock-like 
termination  was  created  next,  and  was  found  to  be  much  quieter  and  tended  to  move  exit 
sources  of  noise  farther  away  from  the  sources  of  interest  on  the  test  plate. 
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The  sidewall  contamination  is  reduced  by  incorporation  of  suction  along  the  test 
plate  following  Motohashi  and  Blackwelder  (1983).  A  narrow  (6  mm)  slot  along  the  side 
of  the  test  plate  provides  a  gap  in  which  some  of  the  over-pressurized  air  within  the  mylar 
tube  is  allowed  to  discharge  to  eliminate  secondary  flows  at  the  plate-tube  junction. 
Therefore,  a  much  greater  percentage  of  the  flow  over  the  test  plate  remains  laminar 
before  wall  contamination  extends  to  the  center  of  the  plate. 

Care  was  taken  to  make  all  comers  within  the  settling  chamber  and  nozzle  as 
smooth  as  possible  to  help  maintain  laminar  irrotational  flow.  Initial  measurements 
indicated  that  a  maximum  velocity  of  31.7  m/sec  at  the  exit  could  be  achieved.  By 
surveying  over  the  plate  with  a  capillary  total  head  tube  on  a  stethoscope  (Figure  5.14), 
the  bursting  typical  of  transition  was  found  to  occur  about  half  way  down  the  plate  which 
indicated  the  beginning  of  natural  transition.  Subsequent  measurements  with  a  tuft  wand 
indicated  that  small  comer  vortices  were  being  created  inside  the  nozzle  at  the  junction  of 
the  nozzle  wall  and  the  flat  bottom.  To  eliminate  this,  the  test  plates  were  redesigned  to 
include  an  extension  of  the  plate  15  cm  into  the  nozzle.  The  entire  plate  was  raised  about 
2.5  cm  from  its  previous  position  to  create  a  1.2  cm  high  by  30.5  cm  wide  duct  below  the 
plate  through  which  wall  flow  (including  the  vortices)  could  be  removed  via  a  suction 
manifold.  Suction  for  flow  management  on  the  test  plate  was  provided  by  the  blower 
intaking  air  from  a  10  cm  hose  attached  to  the  duct  between  the  muffler  and  the  rubber 
isolation  hose  into  the  blower.  The  manifold  was  developed  to  provide  a  more  even 
suction  distribution.  This  is  shown  in  Figure  B.4.  The  manifold  itself  is  shown  in  Figure 
B.5.  The  end  plate  is  removed  in  one  of  the  photographs  to  show  the  rod  mechanism  for 
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opening  and  closing  the  orifices  with  wooden  cones. 

When  designing  the  extension  of  the  test  plate  into  the  nozzle  care  was  taken  to 
avoid  leading  edge  separation  and  premature  transition  by  incorporating  a  4:1  elliptical 
leading  edge  (as  suggested  by  Gad-el-Hak  et  al.,  1981)  and  the  windsock  resistance  at  the 
end  of  the  plate.  This  assured  that  the  stagnation  line  near  the  leading  edge  was  on  the 
upper,  working  surface. 

Provision  was  made  for  either  mechanical  or  wind  puff  turbulent  structure  trips  to 
be  attached  below  the  test  plate  and  reside  in  the  laminar  region. 

Despite  all  of  this  effort,  it  was  determined  that  this  configuration  did  not  meet  the 
low  background  noise  levels  required  to  make  transition  zone  noise  measurements. 

Figure  B.6  shows  a  sound  pressure  level  measurement  made  above  a  natural  transition 
zone  with  a  flow  velocity  of  1 1.8  m/s.  An  early  version  of  GLPROG81  was  used  to 
calculate  predicted  levels  for  transition  noise.  This  required  the  rather  unlikely  values  of 
normalized  rise  time  to  be  between  0.05  and  0.15  in  order  to  match  the  measurement. 
Given  the  shape  of  the  spectrum,  it  is  likely  that  other  contaminating  noise  sources  where 
measured.  Figures  B.7  provides  the  prediction  for  the  same  flow  conditions  using  the 
latest  version  of  GLPROG81. 

B.4  Relocation  of  Test  Section  and  New  Nozzle  Design 

After  the  success  of  Brungart  (1997)  operating  with  the  flow  management  section 
and  square  nozzle  located  in  the  anechoic  chamber  (Figure  B.8),  it  was  decided  to  keep 
this  configuration  intact  and  add  an  additional  plenum.  The  intention  was  to  further 
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reduce  the  low  frequency  background  noise. 

However,  the  46  cm  x  46  cm  outlet  nozzle  needed  to  be  replaced  with  one  that 
provided  a  15.2  cm  x  30.5  cm  outlet.  The  design  of  this  nozzle  is  described  in  section 
5.1.3. 

B.5  Diffuser  Design 

A  new  diffuser  needed  to  be  designed  since  the  Brungart  configuration  did  not 
utilize  one.  The  intent  was  to  slow  the  flow  and  move  any  possible  trailing  edge  noise 
sources  as  far  as  possible  from  the  test  section  measurements. 

Kline  (1959)  and  Kline,  et  al.  (1959)  were  used  to  design  the  diffuser.  The 
Mathcad  Hie  for  calculating  and  specifying  the  nozzle  contours  is  provided  in  section  A.6. 
In  the  design  of  a  straight  walled  diffuser,  one  needs  to  be  aware  of  the  potential  for  a 
stalled  flow  to  develop  in  the  diffuser.  Kline  (1959)  provided  a  criterion  for  determining 
if  stall  might  exist.  This  is  shown  in  Figure  B.9  as  the  line  on  the  plot.  The  diffuser  went 
through  several  deign  iterations  which  are  listed  in  Figure  B.9.  The  open  circles  refer  to 
the  design  condition  for  the  side  walls  and  the  dark  circles  refer  to  the  top  and  bottom 
walls.  A  natural  diffusion  angle,  0,  of  7  degrees  was  designed  first;  however,  a  steady 
stall  condition  resulted.  The  diffuser  length  was  decreased  by  30%  but  the  stall  was  still 
present  The  diffuser  was  then  reconstructed  with  a  4.5  degree  diffuser  angle.  Transitory 
stall  was  seen  depending  on  the  measurement  openings  in  the  test  section.  When  the 
Soundcoat  ®  foam  was  installed,  this  reduced  the  effective  angle  to  3.5  degrees. 
Transitory  stall  still  occurs  at  times  but  it  is  a  very  weak  circulation.  Barlow  et  al.  (1999) 
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recommend  an  area  ratio  of  the  outlet  to  the  inlet  of  5  or  6:1  to  avoid  separations.  This 
diffuser  has  an  area  ratio  of  7.3  which  may  be  why  it  is  just  on  the  edge  for  this  transitory 
flow.  However,  the  desire  was  to  keep  any  potential  trailing  edge  noise  sources  as  far 
from  the  test  section  as  possible. 

The  possible  consequences  of  a  design  point  being  near  the  stall  criterion  line  are 
shown  in  Figure  B.10.  The  steady  flow  condition  should  have  a  uniform  flow  diffusing 
following  the  open  area  ratio  determined  by  the  actual  geometry  of  the  diffuser.  The 
transitory  stall  condition  may  occur  due  to  any  fluctuating  disturbance.  This  is  the  case 
when  one  of  the  measurement  slots  in  the  test  section  is  opened.  The  steady  stall 
condition  effectively  reduces  the  flow  area  of  the  diffuser  resulting  in  higher  discharge 
velocities  (and  mote  noise).  However,  the  addition  of  a  resistance  to  the  end  of  the 
diffuser  (in  the  form  of  four  1.2  cm  thick  furnace  filters)  results  in  a  sufficient  pressure 
rise  that  eliminates  the  possibility  of  this  stall  under  the  run  conditions  for  measurements 
made  in  this  investigation. 
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Figure  B.l  Original  orientation  of  Quiet  Airflow  Facility  to  the  flow-through  anechoic 
chamber. 
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Figure  B.2  Comparison  of  narrowband  signature  of  Borg-Wamer  inverter  and 
replacement  ABB  inverter  for  (a)  upstream  plenum  and  (b)  airborne 
noise  in  control  room. 
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FLOW  STRAIGUTNER  PLENUM  WEATHER 


Figure  B.3  Original  semi-circular  wall  jet  over  flat  plate  configuration. 
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Figure  B.4  Addition  of  mylar  tunnel  above  and  suction  manifold  below  test  plate. 
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Figure  B.5  Views  of  suction  manifold  ports  and  rods  with  cones. 
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RADIATED  NOISE  DUE  TO  FLUCTUATING  DISPLACEMENT  THICKNESS 
FROM  BOUNDARY  LAYER  TRANSITION  USING  LIEPMANN  ANALOGY 
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Figure  B.7  GLPROG81  sound  pressure  level  predictions  for  U=  11.8  m/s  and 
several  \ictJAx  values  otherwise  matching  Figure  B.6. 
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Figure  B.8 


Relocation  of  facility  and  addition  of  moving  wall  test  plate  [from 
Brungart  1997]. 


285 


L/Wt 


Condition  1  7  degree,  2.03  m  diffuser  walls 

Condition  2  7  degree,  1.42  m  diffuser  walls 

Condition  3  4.5  degree,  1 .42  m  diffuser  walls 

Condition  4  3.5  degree,  1 .42  m  diffuser  walls  (after  foam  panel 

insertion) 

Open  circles  Diffuser  side  walls  (W=0.305  m) 

Dark  circles  Diffuser  top  and  bottom  walls  (W=0. 153  m) 


Figure  B.9  Diffuser  design  iterations  versus  stall  design  criterion  [stall  curve  from 
Kline  1959]. 
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(a)  Steady  flow  condition 


(b)  Transitory  stall  condition 


(c)  Steady  stall  condition 


Figure  B.10  Possible  flow  conditions  in  diffuser. 
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APPENDIX  C 


REAL-TIME  ACOUSTIC  INTENSITY 
SIGNAL  PROCESSING 

This  appendix  describes  the  real-time  acoustic  active  intensity  signal  processing 
performed  during  the  air  flow  experiment.  The  processing  uses  inputs  from  two 
microphones  using  a  pressure  -  pressure  (P-P)  configuration.  The  two  pressures  are 
averaged  to  represent  the  pressure  at  a  location  intermediate  to  the  two  microphones  and 
acoustic  particle  velocity  is  derived  by  a  finite  difference  integration  of  the  pressures. 

This  is  predominately  performed  in  a  time-averaged  sense.  Procedures,  governing 
equations,  and  uncertainty  analysis  are  well  documented  in  Fahy  (1989). 

The  general  concept  for  making  measurements  of  the  acoustic  signature  of  the 
transition  zone  or  individual  turbulent  spots  is  shown  in  Figure  C.l.  Transition  is  allowed 
to  occur  naturally  or  a  laminar  flow  condition  is  set  and  a  stimulation  is  used  to  trigger  a 
turbulent  spot.  This  stimulation  can  be  an  acoustic  impulse  from  a  speaker  or  a 
mechanically  generated  air  puff.  A  flush  mounted  hot  film  is  used  to  provide  a  spectral 
conditioning  signal.  The  microphones  may  be  flush  mounted  as  shown  in  the  figure  or 
mounted  outside  of  the  tunnel.  The  idea  is  to  have  the  microphone  separation  large 
enough  so  that  there  will  be  little  correlated  flow  noise  between  them.  Also  by  orienting 
the  two  microphones  in  a  line  normal  to  the  tunnel  axis,  the  null  of  the  cosine  directivity 
pattern  associated  with  the  intensity  probe  will  reject  plane  wave  tunnel  noise  propagating 
along  the  tunnel  axis. 
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C.l  Time-Averaged  Intensity  Processing 

The  traditional  method  of  time-averaged  processing  a  P-P  probe  is  shown  in 
Figure  C.2.  The  active  intensity  for  a  time-stationary  sound  field  is  obtained  from: 

'»  =  ■(^b),m(G^(“))  ecu 

An  alternative  method  often  used  for  handheld  1/3  octave  intensity  measurement 
devices  is  shown  in  Figure  C.3.  Here  sum  and  differencing  is  done  to  the  two  analog 
microphone  signals  before  they  are  passed  through  1/3  octave  filters  and  the  difference 
signal  is  analog  integrated.  Analog-to-digital  conversion  is  then  performed  before 
multiplication  and  division  by  p0d.  Averaging  is  often  then  performed. 


C.2  Real-Time  Intensity  Processing 

An  approximation  for  the  instantaneous  particle  velocity  taken  from  the 
component  of  the  pressure  gradient  in  a  direction  n  is  given  by: 


(C.2) 


With  In(t)  =  p(t)  u„(t),  the  instantaneous  active  intensity  component  is  approximated  by: 


)(Pi (f)+p2(t)\  | \px(i)-p2{x)]di 


(CJ) 


An  example  of  the  analog  integrator  needed  to  implement  Equation  C.3  is  given  in 
Figure  C.4. 
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The  portions  of  the  processing  diagram  prior  to  the  A/D  conversion  have  been 
implemented  in  an  electronic  circuit.  The  design  of  this  circuit  and  its  fabrication  was 
performed  by  Todd  D.  Batzel  of  the  Applied  Research  Laboratory  at  Penn  State.  The 
circuit  diagram  is  shown  in  Figure  C.5.  The  component  layout  is  provided  in  Figure  C.6. 

This  circuit  has  two  analog  inputs  which  can  be  up  to  5  volts.  A  high  pass  filter 
has  been  implemented  with  a  cut-off  at  30  Hz.  The  filter  is  a  fourth  order  Bessel  filter 
with  linear  phase  (Burr-Brown  Universal  Active  Filter,  UAF42G).  This  was  chosen  for 
this  particular  application  in  order  to  remove  some  low  frequency  tunnel  noise  that 
otherwise  would  have  used  up  about  40  dB  of  dynamic  range.  The  gains  in  the  circuit  are 
as  follows: 

multiply  by  250  in  input 

multiply  by  295  in  integrator 

divide  by  10  in  multiplier 

total  gain  =  multiply  by  7375  or  38.7  dB 

A  reset  button  is  installed  for  the  integrator  to  allow  re-zeroing  in  case  of  a  built  up  bias. 
The  circuit  can  be  powered  by  two  onboard  9V  batteries  or  an  AC  adapter. 

The  remainder  of  the  processing  is  performed  in  the  DSP  Technologies,  Inc 
SigLab  ®  20-42  dynamic  signal  processor  which  uses  The  Mathworks,  Inc  MATLAB  ® 
software. 

It  is  important  for  the  signal  processing  of  a  P-P  probe  to  have  two  channels  well 
matched  in  amplitude  and  phase.  The  output  of  the  real-time  intensity  circuit  should  be 
zero  for  the  same  signal  being  input  to  both  channels.  Any  residual  intensity  is  an 
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indication  of  an  error  which  consists  of  random  error  and  bias.  The  “residual  intensity” 
as  compared  to  the  same  signal  input  to  channels  1  and  2  of  the  intensity  circuit  is  shown 
in  Figure  C.7.  The  upper  plot  shows  the  phase  and  amplitude  matching  of  the  input 
signals  even  with  any  channel  mismatch  on  the  SigLab  analyzer.  The  wide  variation  in 
phase  and  amplitude  with  frequency  indicate  problems  in  the  circuitry  that  were  not  able 
to  be  overcome  without  a  complete  redesign.  The  impact  can  be  seen  in  Figure  C.8  with 
a  comparison  of  the  same  measurement  for  the  case  of  a  fully  turbulent  flow  condition 
between  the  time-averaged  active  intensity  method  and  the  real-time  intensity  circuit.  A 
contributing  problem  may  be  the  closeness  of  the  probe  to  the  source  in  these 
measurements.  The  highly  reactive  acoustic  field  may  not  be  conducive  to  the  real-time 
measurement. 

C3  Conditioning  with  a  Shear  Stress  Sensor 

It  is  desired  to  apply  the  power  of  coherent  output  processing  to  the  acoustic 
measurements.  In  this  case,  it  is  assumed  that  shear  stress  due  to  passage  of  a  turbulent 
spot  is  caused  by  the  same  physical  phenomenon  (i.e.  growth  of  the  turbulent  spot)  that 
results  in  acoustic  radiation  (represented  by  the  transfer  function  h(x).  Therefore,  there 
should  be  coherence  between  thiose  signals.  The  noise  associated  with  each  sensor  is 
assumed  to  be  uncorrelated.  These  are  identified  in  Figure  C.9.  It  has  been  demonstrated 
in  other  experiments  that  the  shear  stress  sensor  (hot  film)  is  insensitive  to  background 
acoustic  noise. 

This  represents  a  single-input  /  single-output  problem.  The  measured  coherence 
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function  for  the  shear  stress  and  the  active  acoustic  intensity  is  (Bendat  and  Piersol, 
1993): 


V  0 = 


(C.4) 


where  the  hat  represents  a  measured  value.  Note  that  this  same  approach  can  be  used 
when  a  single  pressure  microphone  is  substituted  for  the  intensity  probe  microphone  pair. 
It  is  assumed  that  the  cross  spectral  terms 


*  G,„(/>  *  OJi)  -  0 


(C.5) 


It  is  also  known  that 


I <v»l2  -  |Gv,tf)|!  =  \H(f>\%(f)  -  av<f) 


(C.6) 


When  the  measured  quantities  are  resolved  into  true  value  and  added  noise,  the 
measured  coherence  becomes: 


v,» ' 


cv y> a »sf> 


l G,,V)+Gmnm  iav(f)*aj/a 


(C.7) 


When  the  measured  intensity  is  multiplied  by  the  coherence  function,  one  gets  the 
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equation: 


v,0)  Gt;if)  -  G,j(j)  [ 


GrJfl 


GtJSf)+GJ[f) 


] 


(C.8) 


Therefore,  a  good  estimate  of  the  true  spectrum  of  the  acoustic  intensity  is  obtained  when 
Gmm(f)«Gr0T0(f)  regardless  of  the  amplitude  of  G^f).  Josserand’s  prior  measurements 
(Josserand  1986)  and  intermittency  measurements  made  during  this  investigation  have 
shown  that  the  shear  stress  measurements  have  a  very  high  signal-  to-noise  ratio  and 
should  meet  the  criterion  well. 
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MEASUREMENT  OF  ACOUSTIC  SIGNATURE 
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Figure  C.  1  General  processing  diagram. 
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Figure  C.2  Traditional  method  of  FFT  processing  of  a  P-P  intensity  probe. 


Figure  C.3  Processing  diagram  for  analog  and  digital  calculation  of  intensity  from  a 
P-P  probe. 
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IF  Ri  =  R2  =  Rj  =  R4 

vp2-p\  =  vp2  ~  vp\ 


Figure  C.4  Processing  diagram  for  integrator. 
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[Circuit  designed  and  constructed  by  Todd  D.  Batzel,  Applied  Research 
Laboratory,  Penn  State  University.] 


Figure  C.5  Circuit  diagram  for  real-time  processing  of  a  P-P  intensity  probe. 
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Figure  C.6  Component  layout  for  real-time  active  intensity  circuit. 


299 


Figure  C.7 


Real-Time  Intensity  Circuit 
Input  J1  to  Input  J2  Transfer  Function 
For  Residual  Intensity  Check 


Output  of  Real-Time  Intensity  Circuit 
to  Same  Input  into  Each  Channel 
Transfer  Function 
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Residual  intensity  measured  for  same  signal  input  to  channels  1  and  2 
compared,  to  input  signal. 
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Figure  C.8  Comparison  of  real-time  intensity  circuit  output  to  time-averaged  intensity 
method  for  a  fully  turbulent  flow  condition. 
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Figure  C.9  Processing  diagram  for  conditioning  a  real-time  processing  of  a  P-P 
intensity  probe  with  a  hot  film  probe. 
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